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IMPACT WAVES IN SAND: THEORY COMPARED WITH 
EXPERIMENTS ON SAND COLUMNS 


By Blaine R. Parkin! 


SYNOPSIS 


A phenomenological theory is developed in order to study the propagation 
of unidimensional compression waves in columns of sand. The medium of the 
theory is treated as an elastic-plastic continuum. It is assumed that each ele- 
ment of the substance exhibits astrain-rate effect such that, at a given strain, 
the plastic strain-rate is proportional to the difference between the compres- 
sive stress on the particle and that stress which would act on the element 
under static conditions. Published experimental results on the propagation of 
stress waves in sand are used as a basis of comparison between experiment 
and theory. The present theory gives satisfactory agreement with experiments 
on stress propagation in two dry sands. 


INTRODUCTION 


From July, 1951, to August, 1954, the United States Army Corps of Engi- 
neers sponsored a program of research into the behavior of soils under 
transient loadings. The work was carried out at the Massachusetts Institute of 
Technology (MIT), Cambridge, Mass., under the direction of the late Donald 
W. Taylor. In 1957 R. I. Whitman,2 M. ASCE, published a brief summary of 


Note.—Discussion open until November, 1961. Separate discussions should be sub- 


mitted for the individual papers in this symposium. To extend the closing date one 
month, a written request must be filed with the Executive Secretary, ASCE. This paper 
is part of the copyrighted Journal of the Soil Mechanics Division, Proceedings of the 
American Society of Civil Engineers, Vol. 87, No. SM 3, June, 1961. 

1 Physics Dept., The RAND Corp., Santa Monica, Calif. 

2 “The Behavior of Soils under Transient Loadings,” by R. V. Whitman, Proceedings, 
Fourth Internatl. Conf. on Soil Mechanics and Foundation Engrg., Butterworths Scien- 
tific Publications, London, 1957. 1 


< 
poe 


2 June, 1961 SM 3 


the principal laboratory techniques and significant results of this work. 
Although Whitman’s paper gives a clear discussion of the scope and main re- 
sults of the MIT experiments, it is still necessary to consult the reports 
originally issued under the contract3,4 in order to view the results in detail. 
From that portion of this work which was concerned exclusively with wave 
propagation, experimental results on Ottawa and Fort Peck sands are re- 
ported, with major emphasis given to wave transmission through dry Ottawa 
sand. Physical descriptions of the two sands found are available.4a 

In these experiments on wave propagation, a typical test assembly con- 
sisted of a cylinder of sand 2 in. in diameter, having a pressure gage in con- 
tact with the sand at each end of the column. In order to give some rigidity to 
this assembly, the sides of the sand column together with portions of the 
pressure cells were encased in a rubber membrane so that most of the air 
could be removed from the sample. Fig. 1 shows schematically a pressure 
gage in contact with one end of a sand column. (Fig. 1 is derived from data 
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FIG. 1.—SCHEMATIC REPRESENTATION OF GAGE 


given by Whitman and Taylor3a,4b.) The test assembly was mounted in a 
special machine which permitted one end of the assembly to be struck by a 
50-lb weight which drove the impact-end pressure gage against the sand. The 
waves, produced in the sand by the motion of the impact gage, traveled the 
length of the sample and were reflected by the stress gage at the reaction 
end of the column. Pressure indications from the gages gave time records of 
longitudinal stress in the sand at the impact and at the reaction ends of the 


3 “The Behavior of Soils under Dynamic Loadings.3. Final Report on Laboratory 
Studies,” by R. V. Whitman, et al., MIT, Dept. of Civ. and Sanitary Engrg., Soil Mechan- 
ics Lab., August, 1954. AFSWP-118, Contract DA-49-129 Eng-227, Office of the Chf. 
of Engrs. ASTIA No. AD-64032. (a) p. 159. (b) Sec. 11.4. (c) p. 184. (d) p. 188, Fig. 
11.9. 

4 “The Behavior of Soils under Dynamic Loadings.2. Interim Report on Wave Pro- 
pagation and Strain-Rate Effect,” by D. W. Taylor and R. V. Whitman, MIT Dept. of Civ. 
and Sanitary Engrg., Soil Mechanics Lab., July, 1953. AFSWP-117, Contract DA-49-129 
Eng-227, Office of the Chf. of Engrs. ASTIA No. AD-28-650. (a) p. 106. (b) p. 189. 
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sample. In addition, position-time records for the impact gage were obtained 
by means of a Schaevitz transducer. The load gage data are as follows: 


Load gage Diaphragm Thickness, T, Computed Capacity, 
in inches in pounds 
No..1 0.045 250 
No. 2 0.071 600 
No. 3 0.119 1,800 


In addition to reporting the data resulting from their experiments, the 
authors also review their laboratory techniques and the characteristics of 
their apparatus. They point out certain experimental difficulties which may 
have limited the validity of some of their results. Even though these difficul- 
ties may have lead to certain inconsistencies in their data, it seems only fair 
to note that experiments of the type undertaken by the group at MIT are 
notoriously difficult.5 These data on the propagation of waves in sand appear 
to constitute the best body of laboratory results now available. 

The purpose of the present paper is to summarize the results of a 1959 
study§ comparing the results of one phenomenological theory for wave propa- 
gation in a slender column of dry sand with experiment. Such an analysis is 
certainly not the first of its kind, and will probably not be the last theoretical 
investigation to use the MIT data as a basis for comparison. For example, 
R. H. Smith and N. M. Newmark, F. ASCE,7 have modeled dry sand by means 
of a theory based upon a viscoelastic-plastic constitutive relationship which 
gives fair agreement with experiment. However, one aspect of the MIT wave 
propagation data which Smith and Newmark excluded from their considera- 
tions consists of a pronounced initial peak or spike in the records of compres- 
sive stress versus time measured at the impact end of the sand columns. 
Such a peak was not found in stress-time records from the reaction end. 

In agreement with the investigators at MIT, Smith and Newmark ascribe 
the occurrence of this stress peak to the action of important lateral inertia 
effects. They contend that such effects are specifically excluded from ele- 
mentary one-dimensional theories of wave propagation, and therefore should 
be ignored in any analysis of the data which uses such one-dimensional theo- 
ries of wave transmission. This contention is a reasonable one. In fact, some 
investigators justify it in a qualitative sense by citing the well-known results 
of Pochhammer and Chree® or other elastic theories which account approxi- 
mately for the effects of lateral inertia.9 Of course, it is generally recognized 


5 “Wave Propagation in Solids,” by E. Pinney, et al., Tech. Report No. 4, Office of 
Naval Research, Contract Nonr-222 (Or) (NR340 040), Dept. of Mathematics, Univ. of 
California, Berkeley, January, 1955, pp. 1-2. 

6 “Impact Wave Propagation in Columns of Sand,” by B. R. Parkin, The RAND Corp., 
Research Memorandum RM-2486, ASTIA No. AD 234961, Santa Monica, Calif., Novem- 
ber, 1959. 

7 “Numerical Integrations for One-Dimensional Stress Waves,” by R. H. Smith and 
N. M. Newmark, Civ. Engrg. Studies Struct. Research Series, No. 162, Tech. Report, 
Office of Naval Research, Dept. of the Navy, Contract Nonr 1834 (03) Project NR064 183, 
Dept. of Civ. Engrg., Univ. of Illinois, Urbana, August, 1958. 

8 “A Treatise on the Mathematical Theory of Elasticity,” by A. E. H. Love, Dover 
Publications, New York, 1944, Art. 201, p. 289. (a) Sec. 309 and Sec. 310c. 

9 “Stress Wave Propagation in Rods and Beams,” by H. N. Abramson, H. J. Plass, 
and E. A. Ripperger, “Advances in Applied Mechanics,” Vol. V, Academic Press, Inc., 
New York, 1958, pp. 111-188. (a) p. 144. 
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that the Pochhammer-Chree or similar results hold only for cylinders of an 
ideal elastic substance which has a well-defined constitutive relationship, 
whereas the constitutive relationship appropriate to any sand is by no means 
well-defined. Therefore when one seeks to analyze the influence of lateral 
inertia as a possible cause for the impact stress peaks observed in cylinders 
of sand, he is forced to rely upon rather crude qualitative arguments. Such 
arguments seem best suited for determining whether or not a given set of 
assumptions concerning the observations is plausible. (An example of such an 
argument is available.3b) Even in those instances where they lend credence to 
a particular theory, such rough arguments certainly donot exclude other sets 
of assumptions which one might wish to make in connection with the interpre- 
tation of a particular body of data. 

Although it cannot be denied, in principle, that effects due to lateral inertia 
may exist, the purpose of the present study is to test an alternative theory by 
comparing numerical results derived from it with experiment. In essence we 
ask if it is possible to exclude all appeals to the effects of lateral inertia and 
to choose a specific constitutive relationship for the medium which, in con- 
junction with a sufficiently detailed one-dimensional representation of the 
experimental arrangement, is able to account for the observed behavior. A 
general study of soil mechanics theories which might apply to the present 
situation will not be undertaken. Such a program is far beyond the scope of 
this paper. It appears that there are still many difficulties which stand in the 
way of an early success with such abroad approach.5,10 Instead a particular 
one-dimensional theory previously used by L. E. Malvern!1,12 in an investi- 
gation of delayed yielding in metals which have been subjected to impact 
loads will be tested. The specific feature of Malvern’s analysis which inter- 
ested this writer was its prediction that immediately after an impact a sharp 
stress peak should appear at the impact end of rods which were made of 
materials obeying the constitutive relationship of the theory. This striking 
qualitative resemblence to the similar stress peak appearing in the MIT 
data suggested that such behavior might be worthy of closer scrutiny. 

Notation.—The letter symbols adopted for use in this paper are defined 
where they first appear, in the illustrations or in the text, and are arranged 
alphabetically, for convenience of reference in the Appendix. 


THE IDEAL MEDIUM 


In an elastic rod, if c denotes the speed of a longitudinal wave it is well- 
known that 


10 “Some Notes on Rheological Properties and Experimental Data for Soils Relevant 
to Theoretical Studies in Plasticity,” by A. D. Cox and H. G. Hopkins, Ministry of Sup- 
ply, Armament Research and Development Establishment, ARDC Memorandum (B) 2/58, 
Fort Halsted, Kent, January, 1958. 

11 “Plastic Wave Propagation in a Bar of Material Exhibiting a Strain-Rate Effect, ” 
by L. E. Malvern, Quarterly of Applied Mathematics, Vol. 8, 1951, p. 405. 

12 “The Propagation of Longitudinal Waves of Plastic Deformation in a Bar of Ma- 
terial Exhibiting a Strain-Rate Effect,” by L. E. Malvern, Transactions, Amer. Soc. of 


Mech. Engrs., Journal of Applied Mechanics, Vol. 18, 1951, p. 203. 
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in which E is the modulus of elasticity andp is the density of the medium. 
Of course sand, under a slowly increasing load, does not exhibit a linear 
elastic behavior at any stress. Even so, the quantities entered in Eq. 1 may 
be given a useful interpretation. It is convenient to definep as the initial dens- 
ity of the medium, and the modulus E will be taken as the slope of a static 
stress-strain curve at zero stress. (Throughout this paper the word stress 
will designate the difference between the total axialstress and the hydrostatic 
pressure exerted by the atmosphere upon the evacuated sample.) This static 
stress-strain curve gives the compressive axial strain resulting from a single 
triaxial test in which the axial load is slowly increased. For the moment, 
stress-strain curves corresponding to unloading or to subsequent loading of 
the specimen are not considered, Eq. 1 then gives the speed of a longitudinal 
stress wave of vanishingly small amplitude. For Ottawa sand appropriate 
values for E and p appear to be 


E = 40,000 psi 


and p = .00016 lb sec? per in.4 


With these quantities c is found, from Eq. 1, to be c = 15,600 in. per sec. The 
preceding values fall within the range given by the MIT group) on the basis 
of their own measurements and also from their study of the static test results 
of L. S. Chen.13 From the MIT results it also appears that the static stress- 
strain curve for Ottawa sand with an initial void ratio of about 0.53 should 
show a maximum strength at 4% strain of about 28 psi and that the stress- 
strain curve should be essentially nonlinear up to this value. Although Ottawa 
sand having the void ratio noted previously was used in most of the wave 
propagation experiments, no static stress-strain data for this void ratio is 
given by Whitman.3 However, in his Fig. 11.1, a static stress-strain curve 
is given for a sand with a void ratio of 0.49 which has a compressive strength 
of 32 psi. In order to derive a stress-strain curve more nearly like that of 
the sand actually used inthe wave propagation tests, the ordinates of the graph 
from Whitman3 were reduced in the ratio of 28 to 32 for axial strains equal to 
or exceeding 0.13%. For values of the compressive strain « in the range 


0 < € <.0013, 
the static compressive stress f(€) was taken as a segment of the ellipse 
(£+3.73)2 (€-.00213)? _ 


4.73x107° 

Eq. 2 provides that f(0) = 0, thatf'(0) = E, and that at « = 0.0013 the slope and 
the ordinate of the elliptic approximation agree with the values obtained from 
the stress-strain curve derived from experiment by the method indicated pre- 
viously. The resulting graph of static stress f versus strain € is shown in 
Fig. 2. This yield curve together with the values of p, E, and c mentioned 
previously constitute the basic parameters defining the medium. 

For purposes of comparison a line through the origin having the slope 
E is shown as the dashed “elastic line” in Fig. 2. Also shown in Fig. 2 is an 
unloading line having a slope parallel to the elastic line which illustrates 


13 “An Investigation of Stress-Strain and Strength Characteristics of Cohesionless 


Soils by Triaxial Compression Tests,” by L. S. Chen, Proceedings, Second Internatl. 
Conf. on Soil Mechanics and Foundation Engrg., Rotterdam, Holland, Vol. 5, 1948, p. 35. 


! 


6 June, 1961 SM 3 


how the material is assumed to behave elastically if the load is removed 
after a static stress state on f(e€) has been established. It will be supposed 
that subsequent reloading up to the yield curve f(¢) will also progress 
elastically. The assumption that points between the static yield curve f(«) 
and the € axis correspond to elastic states of stress in the medium, is a 
drastic simplification. Even at best such ideal behavioris only approximately 
realized in a real sand. Nevertheless, the simplification to be obtained from 
the use of this assumption causes its adoption. If states of compressive stress 
not necessarily on the static yield curve f(€) are denoted by p, then at stress 


Ottawa sand 
Void ratio e= 0.53 


f(e) static compressive stress (psi) 


1?) 0.5 1.0 5 2.0 2.5 3.0 
€ longitudinal compressive strain (percent) 


FIG. 2.—STATIC STRESS-STRAIN CURVE FOR OTTAWA SAND 


points for which p(¢) < f(€) we may characterize the medium by the differential 
form of Hooke’s law, 


Next we may decompose the longitudinal strain € into a plastic strain ¢«P and 
an elastic strain «© so that the total strain is given by 


The elastic strain corresponds to stress points onthe elastic line through the 
origin in Fig. 2 so that 
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whereas the plastic strain corresponds to that part of the strain which is not 
recovered upon release of the load. Consider next a stress point (p,¢) located 
in the region of Fig. 2 between the elastic line and the static yield curve f(e). 
Such a stress, being greater than the static yield stress p = f(€), will tend to 
create a plastic flow of the medium. If inertia effects are excluded it seems 
reasonable to suppose that the time rate of change of plastic straining at an 
instantaneous value of the strain should depend upon the difference between 
the stress p andthe static value f(€) at the instantaneous strain € . If we denote 
time derivatives by dots then, following Malvern,!1,12 we might expect a 
simple equation for the plastic strain rate to be 


when 
p >f(e) 


We note that when p = f(€) the plastic strain rate is zero in accordance with 
our hypothesis that a static stress-strain curve exists. We note also that, 
according to Eq. 6, p is a monotonically increasing function of the plastic 
strain rate when the constant k is positive. Indeed if one wished he might con- 
sider Eq. 6 as thefirsttermina Taylor series for a more general strain rate 
function G(p,¢€), the expansion being made at fixed € about any point on the 
static stress-strain curve. (Some other possible forms of G(p,¢) which have 
received study are discussed by Abramson, Plass, and Ripperger.94) However, 
for the present investigation Eq. 6 will be taken as the appropriate strain- 
rate function for dry sand. We may now combine Eqs. 6, 5, and 4 and the 
approximations leading to Eq. 3 to write the constitutive relationahip for an 
idealized sand as 


k [p-f(e)], p2 fle) 

0 » p< fle) 


In Eq. 7 we have used the partial derivatives with respect to time instead 
of the dot notation in order to emphasize the fact that now the quantities « and 
p are to be regarded as functions of position, x, as well as time, t. In particu- 
lar we will suppose that x is a Lagrange coordinate denoting some initial 
position of the particle under examination. If the longitudinal displacement in 
the positive x direction of this particle at some subsequent time is s(x,t) the 
strain in the medium at that point is given by 


Compressive strains are counted positive. Moreover, the medium at this 
point, being in a state of motion, has a longitudinal velocity given by 


with the result that 
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also holds. Finally it is observed that if p is the initial density of the sand as 
before, the Lagrangian equation of motion for the medium is14 


Thus we see that motions in the medium may be characterized by three 
functions, p(x,t), € (x,t), and v(x,t). Whenp2f(€) these quantities are governed 
by the system of equations, 


and 


These equations constitute a totally hyperbolic system14a having three 
characteristic directions, in a Lagrange or t-x plane, which are given by 
dx = cdt, dx = 0, and dx = - cdt. The quantity c is the elastic wave speed of 
Eq. 1. The following ordinary differential equations hold along such lines: 


dp tpcdv = -k [p - f(c)] dt(on dx = cdt) ....... (18a) 
Ede - dp=k [p-f(c)]dt (on dx=0)........... (13b) 
- dp - pedv = -k [Pp - f(e) | dt (on dx = - cdt) ........ (13c) 


The numerical methods used in the present study are based upon the char- 
acteristic system Eqs. 13 together with appropriate initial and boundary con- 
ditions which must be satisfied if a state of motion in the column is to be de- 
termined without ambiguity. 


INITIAL CONDITIONS AND BOUNDARY CONDITIONS 


In the calculations which follow, interest willbe centered upon reproducing 
the conditions of the MIT experiments. In those tests all motions in the sand 
were started from rest, and the velocity of the impact-end stress gage V(t) 
was given some initial value as a result of being hit by the 50-lb weight. It 
was found that when the 50-lb weight struck the 6-lb stress gage at the impact 
end of the column, the two bodies separated immediately after the collision. 
Thus the gage, which acted as a free body, was pushed against the sand with 
an initial velocity V,>. After impact the larger weight was brought to rest by 
means of suitable stops. In this way further collisions between the gage and 
the larger mass were prevented. 


14 “Supersonic Flow and Shock Waves,” by R. Courant and K. O. Friedricis, Inter- 
science Publishers, Inc., New York, 1945, pp. 30-31, 235-238. (a) pp. 70-75. 
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We may summarize the initial conditions as follows: 


p(x,0) = e(x,0) = v(x,0) =0 ......... 


with 


In order to formulate appropriate boundary conditions it is recalled that 
a test assembly consisted of a column of sand of length 1 witha pressure 
gage on each end. Moreover, all pressure-time records at the impact end 
of the samples showed a continuous increase of stress with time, even though 
the gage at this end of the sample had acquired discontinuously a velocity Vo 
as a result of being struck by the 50-lb weight. If the gage had acted like a 


Fixed support 
Xz 


v(0,t) v(2,t) 


FIG. 3.—IDEALIZED EXPERIMENTAL ARRANGEMENT 


rigid body one would expect the pressure-time records from it to start dis- 
continuously from some nonzero pressure pp, corresponding to the impulsive 
velocity Vo. In the experiments no such records were obtained. The gages do 
not act like rigid bodies. They seem to be compliant bodies and their elastic 
characteristics in the present analysis must be accounted for. Simple cal- 
culations show that speeds of elastic wave propagation in the metal and the oil 
used in these gages are very large compared to the wave speed in the sand, 
and that the main source of flexibility in the gage assembly is probably con- 
tained in the circular plate to which the strain gages are glued (Fig. 1). More- 
over, from the fact that the natural frequency of these gages was found to be 
at least 8,000 cycles per sec, it seems safe to assume here that the gages 
indicate the applied stresses instantaneously. Accordingly, the gages will be 
represented as lumped parameter systems consisting of a massless spring 
connected to a 6-lb weight. The spring of stiffness S is placed between the 
mass and the end of the sand columnas indicated in Fig. 3. It is convenient to 
take the stiffness S in units of pressure rather than in units of force so that 
this quantity will be given in pounds per cubic inch. Similarly the mass m of 
the gage weight will be taken in gravitational units of mass per unit area so 
that the 6-lb weight spread over acircular area 2 in. in diameter gives a value 
for m of 


m = .00494 lb per 
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Of course, it is unnecessary to account for the mass m at the reaction end of 
the test assembly because at this point it is fixed to the ground as shown in 
Fig. 3. In Fig. 3, m denotes the mass of the gage, S refers to the stiffness of 
the gage spring, V(t), is the gate velocity at impact end, v(0,t) describes the 
sand particle velocity at the impact end, and v(1,t) refers to the sand particle 
velocity at the reaction end. 

At the impact end of the test assembly the stress in the spring is p(0,t) and 
the following equations hold at x = 0: 


dv 


and from Eq. 13b, 


Thus, at the impact end of the sand there are four functions to be determined 
by Eqs. 16, 17, and 18 and Eq. 13c. Eq. 13c holds along the backward running 
characteristic between (0,t) and some other point, (x,t,), in the medium for 
which Vo,€, and po are known. 

Similarly at the reaction end x = 1, the equations to be satisfied are: 


and Eq. 13, which holds along the forward running characteristic between 
(1,t) and some other point (x1,t,) inthe medium, At the point (x;,t,) the quanti- 
ties vj, €1, and p; are known. 

Evidently, Eqs. 13 through 20 permit the computation of the functions p, v, 
and € whenever (p-f) = 0. When this inequality is not satisfied Eq. 7 shows 
that the medium behaves elastically, although the possibility that different 
elements in the medium may have undergone different amounts of plastic 
straining requires that the particle paths retain their role as characteristic 
curves in this “elastic” regime. 

These calculations by the method of characteristics were eventually pro- 
grammed for execution by a digital computer. It app-ared that the chosen 
integration scheme leads to a stable numerical procedure which is accurate 
enough for present purposes. For example, when numerical results for a com- 
plete history of the stresses, strains, and particie velocities after an impact 
had been computed, the initial momentum of the impact gage was compared 
with the impulse under the curve for the impact-end stress history from the 
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time t = 0 until the instant at which the impact gage velocity became zero. 
It was found that these quantities differed by an amount which was less than 
1% of the initial momentum of the impact gage. 

In some of the computations which follow only the very first part of the 
impact stress response is required. In such cases it may be seen from con- 
siderations of impulse and momentum that the impact gage velocity V(t) will 
change by only a few percentage points. In those instances where simplicity 
has greater value than precision, the function V(t) will be replaced by its 
initial value V,. In this case Eqs. 16 and 17 may be replaced by the single 
condition, 


= s [vq - v(0t)] .. 


Finally it is observed that two undetermined constants, k and S, are con- 
tained in the theory. The values of these parameters must be determined from 
experiment. In its essentials this determination is carried out by adjusting 
the values of k and S sothat calculated results match the experimental obser- 
vations of the initial impact-stress rise shortly after an impact at a specified 
value of Vj). Since there are two constants to be found, two pieces of experi- 
mental data are needed to determine them. For these data we may take the 
peak stress in the stress spike and the time required for the stress to reach 
this peak. In order to match these datafor a given value of V,, various values 
of k and S were used in successive calculations based on the characteristic 
system, Eqs. 13, and the boundary conditions, Eqs. 16, 17, and 18. Details of 
the computation methods have been given by the author.6 


OTTAWA SAND 


The considerations of the two preceding Sections enable results from the 
theory to be compared with the dynamic stress histories obtained from the 
experiments on Ottawa sand, It is only necessary to select one representative 
impact gage record in order to determine the values of k and S by means of 
the procedure outlined previously. However, examination of the impact stress 
histories reported by Whitman? reveals a marked inconsistency in these data 
which has already been noted by Smith and Newmark.” It happens that the initial 
momenta reported for the impact gage are always significantly less than cor- 
responding values of the impulse obtained from the impact gage stress-time 
records. In fact the ratio of impact impulse to reported values of initial mo- 
mentum has an average value of 1.8, even when that portion of the stress-time 
records corresponding to unloading of the impact gage is excluded from the 
calculations. Although this discrepancy canbe ascribedtoa number of causes, 
there is strong evidence in support of the view that a major contribution to this 
inconsistency arises from the fact that all measurements of impact velocity 
Vo are probably in error.6 For the sake of simplicity in the following discus- 
sion we will assume that all of this discrepancy is due to difficulties en- 
countered in the measurement of the impact velocities. 

As noted previously, the impact gage, after being struck by the 50-lb 
weight, was observed to act as a free body with an initial velocity V,. The 
initial momentum of this gage is expended against the impulse of the reaction 
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force in the sand. Of course some elastic energy is stored in the impact gage. 
This energy together with the elastic energy of the reaction gage and some 
residual elasticity in the sand acts to impart a rebound velocity to the impact 
gage. Therefore the total change in the magnitude of momentum suffered by 
the impact gage and the total impulse under the impact stress curve undoubt- 
edly exceed the initial gage momentum. However, the difference between the 
initial momentum and the total change in momentum is not large, for one can 
estimate the amount of recoverable energy noted previously to show that the 
rebound velocity of the impact gage is generally less than 10 in. per sec. The 
start of the rebound phase of the gage motion is probably near the beginning 
of the unloading portion of the experimental curves of impact stress versus 
time given by Whitman,3 and the gage velocity probably passes through zero 
near the start of this unloading phase. Hence, the area under the impact stress 
versus time curves, up to the start of the unloading phase, should give an im- 
pulse which is nearly equal to the initial momentum of the impact gage. The 
stress-time curves presented by Whitman? were not intended to be precise 


TABLE 1 


7 in, samples 14 in, samples 20 in. samples 32 in, samples 
Yo | Pmax | Yo | Pmax| Yo | Pmax | | Pmax 


25 39 22.5 | 35,7 25.0 33.6 
63.4 33, 63.4 46. 
116, 51 94,7 54 92.4 
155. 


Vo = uncorrected impact velocity, inches per second 


Pues * peak impact stress, pounds per square inch 


% = corrected impact velocity, inches per second 


reproductions of the experimental results. Therefore exact equality between 
the measured impulse and the correct value of the initial momentum can not 
be expected to hold even if our estimate of the time at which the impact gage 
starts to rebound should prove to be perfectly accurate. 

With this understanding, the impulse, measured in the manner outlined 
previously may be used, as the basis for computing corrected values of the 
impact velocity. These corrected values are denoted by the symbol Vo in 
Table 1. In this table values of peak stress and corrected impact velocity 
are given for each sample length andfor each nominal value of impact velocity 
reported for the tests on Ottawa sand. 

In order to find values ofk andS appropriate to these data, a corrected im- 
pact velocity V3 of 80 in. per sec, a peak stress of 45 psi, and a rise time of 
0.49 millisecond as basic parameters in calculations of the sort outlined at 
the end of the preceding Section were chosen. The resulting values of spring 
stiffness and strain rate constant are 


k = 1.68 x 10° per sec 
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S = 2.2 x 103 lb per cu in. 


respectively. The fit to the corrected data which results from the use of these 
particular values of k and S incalculations which permit the corrected impact 
velocity to vary over the range of the experiments, is shown in Fig. 4 as the 
solid curve. In this illustration the curve labeled “Uncorrected” shows the 
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FIG. 4.—PEAK STRESS VERSUS CORRECTED IMPACT VELOCITY 


theoretical variation of peak pressure with impact velocity when the theory is 
fit to the data without the use of any correction procedures. Detailed discussion 
of the “uncorrected” results will be given subsequently. For the present the 
analysis will be confined to calculated results which are based upon the cor- 
rected data of Table 1. In Fig. 4 the curve through the origin which is labeled 
with the equation 
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shows the peak stresses predicted by this result from elementary elastic wave 
theory when values of p andc appropriate to Ottawa sand are used in the cal- 
culations. Evidently this result from elastic wave theory overestimates the 
maximum impact stress by an order of magnitude for most impact velocities. 
Moreover, it can be argued that this result from the elementary elastic theory 
leads to an estimate of the maximum compressive impact stress which is 
greater than or equal to that predicted by the present anelastic theory for all 
non-negative impact velocities. 

The satisfactory agreement between experiment and theory shown by Fig. 4, 
although encouraging, is not entirely sufficient to give credibility to the present 
theory. Since those theoretical calculations were designed with an eye upon 
events at the impact end of the sand column during only the very first part of 
the motion, it may be argued that additional comparisons between experiment 
and theory over longer time intervals and at other points in the sand must be 
carried out. In particular, experimental and theoretical stress histories at 
both the impact and the reaction ends of a typical sample during the entire 
motion must certainly be compared. 

In order to carry out the necessary computations, the gage stiffness S at 
the reaction end of the sample must be determined. Indirect evidence from 
Whitman” leads to the belief® that a gage conforming to the specification for 
gage No. 1 of Fig. 1 was probably used at the impact end ofall samples. (Un- 
fortunately no details were given2,3,4 concerning gage calibration curves, 
various gage stiffnesses, or which gages were used in the tests. In fact the 
total number of gages which were fabricated is not disclosed.) It is possible 
that any one of the three gages of Fig. 1 could have been used at the reaction 
end of these samples. Now if the stiffness of gage No. 1 in combination with 
Ottawa sand has the value of S calculated above, then the stiffness of the other 
two gages may be estimated from the fact that the stiffness of all gages should 
be at least roughly proportional to the cube of the diaphragm thickness, T. 
Some preliminary calculations using these three stiffnesses as well as a rigid 
constraint at the reaction end of a14-in, sample indicate that the use of gages 
No. 1 or No. 2 at this positionis equally likely. We have based our subsequent 
computations upon the use of gage No. 2 at the reaction end. 

Results from a complete set of computations for the response of a 14-in. 
sample of Ottawa sand to an uncorrected impact velocity of 60 in. per sec are 
shown in Fig. 5. The corrected impact velocity used in the computations may 
be identified from Table 1 as Vg = 94.7 in. per sec. The overall agreement 
between experiment and theory is thought tobe quite good, although the theory 
shows no stress plateau near the startofthe stress-time curves from the re- 
action gage. Moreover the duration of the theoretical stress-time curves ex- 
ceed those obtained experimentally by 1.2 milliseconds or about 5%. 

This difference in duration appears to be due in large part to the fact that 
the calculations give values of p(t) at the impact gage which are less than those 
measured experimentally after the peak impact stress has been attained. but 
the computations require that the initial gage momentum must equal the im- 
pulse under the impact p(t) curve from t = 0 until the gage velocity V(t) = 0. 
Therefore if, as is supposed, the initial gage momentum in the computations is 
equal to that of the experiments, the computed response time must exceed that 
measured in the experiments. Of course, the correction method supposes that 
the experimental data concerning the magnitude of p(t) is correct and that the 
estimate of the time corresponding to the condition V(t) = 0 from the experi- 
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mental p(t) records have been correctly made. Corrected impact velocity 
Vo = 92.72 in. per sec. The strain-rate parameter k = 1.68 x 105 per sec. 
Although errors due to both of these causes are undoubtedly contained in the 
computations of Fig. 5 it seems safe to say that they are probably of minor 
importance compared to the gross errors contained in the values of impact 
velocity reported by Whitman.3 
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FIG. 5.—EXPERIMENTAL AND THEORETICAL STRESS HISTORIES FOR A 60- 
IN. PER SEC. IMPACT ON A 14-IN. COLUMN OF OTTAWA SAND 


Aside from the selection of the reaction stress-gage stiffness (which has 
been found to constitute a choice of minor importance), the adjustment of the 
theory to fit the experiments is made to depend mainly upon events measured 
at the impact end of the sample. Therefore the fidelity with which the theory 
reproduces stress histories measured at the reaction end of the sand column 
gives some indication of the theory’s efficacy. For this reason it is of value to 
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compare experiment and theory for impacts on sand columns having various 
lengths, and to this end the early portions of impact-end and reaction-end 
stress histories for a 20-in. column of Ottawa sand are represented in Fig. 6. 
Except for the lengthened column, the values of all parameters which enter 
into the computations of Fig. 6 correspond exactly to those used in the com- 
putations leading to Fig. 5. Again the computed and the experimental results 
show satisfactory agreement. 

Whitman®¢ has presented asummary of data upon longitudinal strain versus 
time for Ottawa sand, These compressive strain measurements were obtained 
by imbedding SR-4 strain gages in the sand at distances from the impact end 
equal to x = 6 in. and x = 9 in. Transient measurements were obtained at the 
6-in. location on 14-in. samples only. The data for the 9-in. location were ob- 
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FIG. 7.—LONGITUDINAL STRAIN HISTORIES IN THE SAND AT 9 IN. FROM IMPACT 
END OF SAMPLE 


tained from 20-in, or 32-in. columns. For both locations data are reported at 
impact velocities of 15, 45, 60, and 90 in. per sec. In all cases data were re- 
corded for a 3 millisecond interval. 

The measurements were obtained by means of strain gage assemblies which 
were approximately 3 in. in length. No corrections to the data which allow for 
the effects of this length are reported. The gages were tested statically by 
embedding them in short samples of sand and subjecting these samples to 
standard triaxial tests. It was found that the gages responded linearly with 
respect to the axial strain of the triaxial samples up to strains of 0.4%, and 
that the static calibration was not linear in the region of compressive strains 
greater than 0.4%. This lack of linearity was attributed to the occurrence of 
slippage between the gages and the sand. 
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The only complete curve of compressive strain versus time reported by 
Whitman? pertains to a 32-in sample which suffered an impact of 15 in. per 
sec (V5 = 42 in. per sec). The measurements were taken at x = 9 in. This 
curve is reproduced in Fig. 7. For these conditions a theoretical curve of 
strain versus time is also shown in this graph. The computations were made 
for corrected values of impact velocity, k and S in accordance with the data 
of Table 1. The dashed part of this curve indicates an extrapolation of these 
results which is based upon similar computations for values of k and S which 
differ slightly from those of Table 1. This dashed portion of the curve is added 
in order to give at least qualitative indication of the theoretical trend as the 
time increases. Within the time intervai shown in Fig. 7, no reflected waves 
have reached the station at x = 9 in. from the impact end of the column. 

Comparisons between experimental and theoretical strains for other sam- 
ple lengths and at other impact velocities have also been made. These com- 
parisons give results which show both better and poorer agreement between 
experiment and theory than is indicated in Fig. 7. These comparisons have 
shown three essential points of difference between experiment and theory. 


1. The measurements indicate the existence of states of constant strain 
shortly after an impact, whereas the theory predicts that a maximum strain 


TABLE 2,—PEAK STRAIN RATES.COMPUTED FOR OTTAWA SAND 


Corrected impact velocity V', 
inches per second 


Maximum strain rate, €| percentage 
max 


per second 


will be reached at a later time and at a higher value of € . However, because 
of the static calibration results noted above we are prepared to accept this 
difference between experiment and theory as a possible indication of slippage 
between gages and the sand. 

2. Once the process of deformation is well started the theory predicts a 
rate of rise which is generally higher than the rates observed. That such a 
discrepancy might be caused at least in part by the effects of the nonzero length 
of the strain gages is supported by some results summarized by R. M. 
Davies. 15 However, no analysis was carried out to estimate the magnitude of 
this effect for gages which are embedded inthe present idealized medium. 

3. The experimental findings suggest that the increase of strain with time 
should be linear whereas results from the present theory show that the function 
€ (x,t) is not linear in time at any of the stations where this behavior was 
studied. In view of the fact that the MIT investigators seem to regard their 
transient strain measurements as being mainly of qualitative value, it may be 


15 “Stress Waves in Solids,” by R. M. Davies, “Surveys in Mechanics,” Batchelor 


and Davies, Eds., Cambridge Univ. Press, 1956. pp. 87-93. 
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that the indifferent agreement between experiment and theory obtained in the 
best case examined is as good as one should expect. The agreement between 
experiment and theory in the case illustrated by Fig. 7 is certainly of a rather 
low order. 


There are no data concerning the strains in the sand at the impact end of 
the samples. Since the strains are applied most rapidly at this location some 
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FIG. 8.—PEAK STRESS VERSUS IMPACT VELOCITY FOR OTTAWA 
SAND (UNCORRECTED DATA) 


information concerning their histories would certainly be of interest. Of 
course, we can obtain theoretical «- t curves fromthe present computations and 
from these curves we can find the maximum theoretical strain rate experienced 
by the sand. Such results are given in Table 2 below where values of the peak 
strain rate are given for some of the corrected impact velocities shown in 
Fig. 4. The entry at a V4 of 94.7 in. per sec corresponds to an uncorrected im- 
pact velocity Vo of 60 in. per sec. The effect of the present correction pro- 
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cedures upon the theoretical peak strain rates may be estimated from the fact 
that at this impact velocity the peak strain rate is 2,400% per sec when no 
corrections are employed. 

Next experimental and theoretical results for Ottawa sand are compared 
when the theory is fitted to the Whitman’ with no corrections being applied to 
them. In order to make this comparison new values of the parameters k and S 
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FIG. 9.—EXPERIMENTAL AND THEORETICAL IMPACT STRESS HISTO- 
RIES FOR A 14-IN. COLUMN OF OTTAWA SAND 
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and 
S = 4.6x 103 lb per cuin. ....... 


In order to test further these values of k and S computations at other im- 
pact velocities were carried out by the method of characteristics. In Fig. 8 
these theoretical results for the variation of peak impact pressure with impact 
velocity are compared with the uncorrected experimental data for Ottawa 
sand. With the range of velocities considered the agreement between experi- 
ment and theory again seems satisfactory. Since nocorrections are applied to 
the data of Fig. 8 the experimental points show somewhat less scatter than do 
the corrected values plotted in Fig. 4. The curve labeled “Uncorrected” in 
Fig. 4 is simply another plot of the theoretical curve of Fig. 8. 

If the stiffness of gage No. 2 is adjusted to account for the fact that gage 
No. 1 must now have the preceding value of S the computations which lead to 
Fig. 5 may be repeated. Results from such computations for a 14-in column 
of Ottawa sand responding to an impact at 60 in. per sec are shown in Fig. 9. 
Uncorrected impact velocity, V_ = 60 in. per sec. Strain-rate parameter, k = 
2.62 x 104 per sec. It willbe seen that the computed peak stress at the impact 
end of the cylinder now exceeds the experimental value by 8 psi. This differ- 
ence is close to the magnitude of the experimental scatter shown in Fig. 8. 
In fact, the particular experimental situation chosen for comparison may be 
identified as one of the plotted points in Fig. 8. Thus, in Fig. 8 a group of 
three experimental points near an impact velocity of 60 fps are seen. The 
lowest of these points results from the experiment ona 14-in, sample. 

There are two essential differences between experiment and theory which 
are shown by the curves of Fig. 9. The first difference concerns the total 
duration of the impact process, Clearly, the theory significantly under- 
estimates this duration and the present example is presented chiefly to em- 
phasize this difference. The cause for this underestimate appears to reside 
with the inconsistencies in the data which have been previously examined.3 
The second difference between experiment and theory is shown by the reaction 
gage curves in Figs. 5, 6, and 9. The experiments show an interval of nearly 
constant stress at the reaction end of the sample shortly after the start of its 
motion. Such stress plateaus are not accounted for by the present form of the 
theory although similar jogs in the reaction end stress-time curves are ex- 
hibited by all experimental curves reported for Ottawa sand. No satisfactory 
explanation for these observations has been advanced. However, the present 
study has led to a tentative hypothesis concerning this behavior which can at 
least be given theoretical justification. 

That it might be possible to alter the formulation of the present theory 
so as to account for this feature ofthe reaction gage response, was the result 
of the study of theoretical calculations in which the smooth nonlinear static 
stress-strain function f(€) was replaced by broken-line approximations. It was 
found that such approximate computations lead to stress-time curves at the 
impact end and in the interior of the sand which differ only slightly from 
similar results obtained with the smooth stress-strain curve of Fig. 2. How- 
ever, for the broken-line theory it was found that the stress-time curves 
computed at the reaction end contained jogs which bore some resemblence 
to those observed experimentally. Moreover it was found that these jogs in 
the reaction gage response were definitely associated with the corners in the 
static stress-strain curve which are contained in any broken-line approxima- 
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tion of the f(€) curve shown in Fig. 2. Physically, the presence of a corner in 
a static stress-strain curve represents a static stress at which the rate of 
slippage in the idealized medium suddenly increases. Therefore it becomes 
necessary to examine the following extreme case of sudden slippage in Ottawa 
sand. 

Suppose that the sand behaves elastically as the static stress increases 
from zero up to f = 5 psi. Suppose next that a regime of perfect plasticity is = 
maintained at f = 5 psi for strainsin the interval .000125 < ¢€ < .000512. Upon 
further straining suppose that f(¢€) increases linearly with one-fourth the 
elastic slope. Suppose further that the other parameters defining this idealized 
sand are given by the preceding values of k and S, and those which follow Eq. 1. ° 
The resulting responses of the impact and reaction gages to a 60 in. per sec 
impact are shown in Fig. 10. This figure also contains a plot which compares 
the smooth static stress-strain curve from Fig. 2 (dashed curve) with the 
curve which has the above region of slip (solid lines). Except for the differing 
static stress-strain relationships the theoretical conditions pertaining to Figs. 

9 and 10 are identical. The close agreement between the impact stress-time 
curves in the two cases is clear. Evidently the introduction of a small region 
of pure elasticity into the computations of Fig. 10 results only in a slight in- 
crease in the peak impact stress over that shown in Fig. 9. Of more impor- 
tance, Fig. 10 clearly shows a “stress plateau” at the reaction end of the 14 in. 
sample which bears a strong resemblance to that found experimentally. 

Of course, once the possibility of sudden slippage is admitted into the 
theofy the detailed manner in which this phenomenon is introduced in any 
particular case is subject to wide variations. If such slippage actually takes 
place in experimental situations, the manner of its occurrence might also be 
subject to considerable variations. However, it seems clear that the present 
theory containing the possibility of slippage can be made to reproduce certain 
of the experimental results with considerable fidelity. 

_ Finally, it must be noted that eventhough it has been found possible to cor- 
rect the MIT data so that excellent agreement between experiment and theory 
can be achieved, we shall make no further use of such correction procedures. 
The main reason for avoiding the further use of such procedures lies in the 
fact that the correction method allows for only one source of inconsistency in 
the data. Undoubtedly, there are other sources of uncertainty in these experi- 
ments. For this reason this paper will be restricted to the use of the MIT 
data as it is reported. By this means further discussions concerning the ap- 
propriateness of one or another correction scheme inany particular situation 
is avoided, and the ability to make useful qualitative comparisons is retained. 


FORT PECK SAND 


Some data on the propagation of compression waves in cylinders of Fort . 
Peck sand are given by Whitman.3 Details concerning static stress-strain 
curves for this sand are presented by Taylor and Whitman.* From these 
sources we find that Fort Peck sand, whenin a dense state with a void ratio of 
about 0.70, exhibits a peak compressive resistance of about 65 psi under a 
confining pressure of one atmosphere. This peak strength is reached at a 
strain of about 4%. Further straining decreases the resistance to about 40 psi. 
Clearly, the static peak strength of Fort Peck sand differs significantly from 
that of Ottawa sand. The primary reasons for this difference are attributed to 
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FIG. 10.—CALCULATED STRESS HISTORIES USING SMOOTH STATIC 
STRESS-STRAIN CURVE (DASHED CURVES) AND A STATIC 
STRESS-STRAIN CURVE WITH SUDDEN SLIPPAGE (SOLID 
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the facts that the grains for Fort Peck sand are more angular than are those 
of Ottawa sand and that the two sands have different grain size distributions. 
These two factors contribute to the phenomenon called “interlocking.” Thus, 
if one wishes to shear a sandy soil it is not enough to apply sufficient force 
to make the grains slide over one another. Additional force must be applied to 
the soil mass in order to displace the grains from their interlocked structure 
so that they can slide past one another. 
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FIG. 11.—STATIC STRESS-STRAIN CURVE FOR FORT PECK SAND 


For the present computations a static stress-strain curve for Fort Peck 
sand has been derived by the same methods as were used previously for 
Ottawa sand. Thus, data on wave speed and density were used to determine 
a representative value of Young’s modulus. Values which were found to be 
consistent for Fort Peck sand with a void ratio of 0.70 are 


E = 21,4000 psi, 
.00015 Ib sec2 per in.4 
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—_ c = 12,000 in. per sec 
For strains in the interval 0 < € < .009 it was found possible to approximate 
the static stress-strain curve by the parabola 


f= 6.902[1 + - 1]. ............. (24) 


This parabola was chosen to have the slope E at zero stress and it was found 
to give agreement with the experimental data inboth ordinate and slope of the 
stress-strain curve at € = 0.009. The resulting stress-strain curve which 
forms the basis for the present calculations is shown in Fig. 11. The unstable 
part of the static stress-strain curve, corresponding to strains which exceed 
4%, are excluded from the present considerations, The considerable difference 
in shape between this curve and the corresponding curve for Ottawa sand is 
noteworthy. 

As before, appropriate values of gage stiffness S and strain-rate parameter 
k must be added to the preceding data before detailed computations for wave 
propagation can be carried out. However, the value of S which describes the 
stiffness of the impact gage has already been determined from the Ottawa 
sand data. If this gage suffered no alteration between the tests on Ottawa and 
Fort Peck sands and if other factors not considered in the formulation of the 
theory do not come into play, it only remains to find that value of k which, for 
a specified impact velocity V,, reproduces the experimentally observed value 
of peak impact stress, This matching was carried out by successive numerical 
integrations of Eqs. 13, 16, and 17. The value of the strain-rate parameter 
which gives an adequate fit to the data is 


k = 6.5 x 103 per sec 
It will be recalled that the gage spring stiffness for this case is 
S = 4.6 x 103 lb per cu in. 


The value of k just given is based upon the condition that when V, = 60 in. per 
sec the peak pressure is 


pj = 77.8 psi 
The theoretical rise time required to attain this peak is 
t; = 0.65 milliseconds 


The theoretical rise time noted above must be compared with the rise time 
obtained in the experiments on Fort Peck sand. The experimental values of 
this quantity appeared never to be less than 0.8 millisecond and never to be 
more than 1.1 milliseconds. The theoretical rise time can be brought into 
better agreement with experiment if the value of the gage stiffness S is low- 
ered from the value which was found to be appropriate for Ottawa sand. 
Lowering the value of S also necessitates lowering the value of k in order to 
preserve the value of the peak stress. It was found that when these two para- 
meters were given the values of 


k = 3.0 x 103 per sec 
and 


S = 3.0 x 103 lb per cu in. 


SM 3 25 
. 
| 
| 


26 June, 1961 SM 3 
the experimental values of both peak stress and impact velocity were fitted - 
satisfactorily. To obtain these values the condition chosen for the fit was 
130 
Experiment (Whitman >) 
© 7 in. somples 
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100 
a 90 
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FIG. 12._PEAK STRESS VERSUS IMPACT VELOCITY FOR FORT PECK SAND - 
pj = 79.3 psi, V,, = 60 in. per sec Bs 
and 
t; = 0.97 millisecond 
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After obtaining the two pairs of values for Sandk described previously 
computations were made based upon them for the additional impact velocities 
of 15, 30 and 90 in. per second. Comparisons between these theoretical pre- 
dictions and the available data for Fort Peck sand are shown in Fig. 12 where 
the peak compressive stress is shown plotted against impact velocity. The 
dashed curve shows the theoretical relationship between the peak stress and 
Vo when the spring stiffness S is retained at that value determined from the 
experiments on Ottawa sand. The solid curve shows the result obtained when 
k and S are determined by matching both peak stress and rise time for one 
impact in Fort Peck sand. 

The fact that a new value of gage stiffness S is required to match the ex- 
perimental results each time the medium is changed is an important result 
which requires some consideration. It suggests that there is an important 
interaction between the gage and the sand whichis not explicitly accounted for 
in the formulation of the present theory. The nature of this interaction may be 
ascertained by recalling that a number of assumptions are implied in the 
definition of S. Most important among these is the implication that the gage 
stiffness somehow may be characterized by a uniform displacement and a 
uniform stress acting over the surface of the gage diaphragm. We may gain 
some idea of how imprecise these notions are by means of some simple cal- 
culations. 

In accordance with the present theory suppose for the moment that at the 
impact end of the cylinder the surface of the sand which is plane before impact 
remains plane during an impact. For animpact at 60 in. per sec the displace- 
ment of this surface between the time t = 0 and the time at which the peak 
stress is reached is about 0.015 in. inthe case of Ottawa sand and about 0.030 
in. in the case of Fort Peck sand. Next let us suppose, in conformance with the 
present theory, that this displacement gives rise toa peak pressure in the sand 
which is uniform over the impact end ofthe cylinder. At an impact velocity of 
60 in. per sec this peak pressure should have values of 60 psi for Ottawa sand 
and 78 psi for Fort Peck sand. If now the elastic plate theory8@ is used, it is 
found that as a result of these pressures the center of the impact pressure 
gage (gage Ne. 1 of Fig. 1) should suffer a maximum deflection of at least 
0.011 in. for Ottawa sand and a maximum deflection of at least 0.016 in. 
for Fort Peck sand. These values suppose that the gage diaphragm behaves 
like a rigidly-clamped circular plate. The actual gage deflections prob- 
ably exceed these values. However the important point, so far as the pres- 
ent examination is concerned, is that for both Ottawa and Fort Peck sands it is 
found that the displacement of the sand at the impact end of the column and the 
deflection of the impact-end pressure gage are approximately equal. This find- 
ing implies that the impact surface of the sand does not actually remain plane 
and that the pressure across it is not actually uniform as assumed. Evidently 
the experimental impact boundary condition corresponds only approximately to 
that required for one-dimensional wave propagation. Therefore, itis not sur- 
prising that the simplified gage stiffness of the present theory shows a depend- 
ence upon the detailed character of the medium under study. Perhaps one fact 
which should cause some surprise is that the present theory, although dras- 
tically simplified, is as successful as it has been found to be. However, this 
success should not blind the researcher to the fact that the lack of precise 
correspondence between the experimental boundary conditions and those re- 
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quired by the theory leads to additional undertainty in our determinations of 
values for the strain rate parameter k. 


OTHER MIT EXPERIMENTS 


All theoretical work discussed previously has sought to interpret the MIT 
wave propagation experiments. It is remarkable that the very simple theory 
formulated here has been so successful inthis task. However the wave propa- 
gation work does not constitute the bulk of the MIT experimental effort. There- 
fore it is of interest to examine some of the other experiments to see whether 
or not the present theory is consistent with results obtained from them. 

A major part of the MIT experimental program was carried out with an 
hydraulic apparatus and special instrumentation which was designed to measure 
the transient stress response of triaxial soil samples at predetermined rates 
of strain. From such tests, curves of compressive strength versus strain rate 
were determined for cohesive soils, dry sands, moist sands, and saturated 
sands. Ideally, the purpose of the tests was to study the relationship between 
shearing strength and rapidity of loading by subjecting one end of a short sam- 
ple of soil to a compressive displacement at constant velocity for various 
periods of time. It seems likely that the measured velocities were divided by 
the length of the unstrained sample inorder to determine a “strain rate.” For 
example, in the case of experiments with dry Ottawa sand the sample length 
was 3-1/2 in. so that data for a “strain rate” given as 3,000% per sec would 
correspond to a loading velocity of 105 in. per sec, and “strain rates” of 700% 
per sec would correspond to a loading velocity of 24.5 in. per sec, The time 
of loading reported for an average velocity of 40 in. per sec was 1 milli- 
second.2 Presumably, for lower velocities the loading time was proportionately 
longer. 

In order to measure the load on the soil, ring gages were mounted at the 
so-called piston end and at the reaction endof the sample. A schematic repre- 
sentation of the sample and the gages might look something like the diagram 
of Fig. 3 except that the ring gages, together with associated parts of the ap- 
paratus which transmits the force between gages and soil, require a somewhat 
different lumped parameter representation; and the impact velocity V(t) at the 
impact or piston end of the sample would now be replaced by a constant veloc- 
ity V5, which acts during the time of loading. At the time of the experiments 
on dry Ottawa sand it seems likely that the natural frequency of the gage at 
the piston end was about 2,000 cycles per sec and the natural frequency of the 
reacticn gage was about 1,500 cycles per sec. For reasons which are fully 
set forth by Taylor and Whitman‘ when the hydraulic rapid loading machine 
was put into operation it was found that these gages were subject to ringing. 

This behavior of the gages made the determination of strain rate effects 
on dry sand rather difficult. According to Whitman? no significant strain-rate 
effect could be detected for dry sand. However, owing to the difficulties pre- 
viously noted, to nonlinearities in the recording apparatus, and to unavoidable 
inconsistencies in sample preparation, these tests did not disprove the exist- 
ence of small strength increases. Hence, carefully performed creep and re- 
laxation tests were used for further study. These tests led to the conclusion 
that the compressive strength of sands increases about 10% to 15% between 
times of loading which range from several seconds to 50 milliseconds, It 
would appear that no further significant increase occurs at the time of loading 
decreases to 5 milliseconds. 
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Of course, it is not at all obvious that the theory of the present study is 
directly applicable to test results from the hydraulic apparatus primarily be- 
cause of the relatively short samples tested. Nevertheless it will be assumed 
that it is possible to compare the experimental results from this machine with 
predictions from the present theory in order to obtain qualitative information. 
Then if the frequencies of ringing are compared with the rise times computed 
at the impact end of the sample, it is concluded that there is a strong possibil- 
ity that the experimental records from the piston gage would be modified in 
such a way that the significant strain rate effects predicted by the theory 
would be entirely masked. Moreover, because of the strong dissipation in the 
idealized medium, the theory shows that by the time the initial spike has 
traveled the 3-in. sample length its magnitude will most likely have decayed 
to such an extent that little if any strain rate effect would be detected by the 
reaction gage. Therefore, in the case of the rapid loading experiments, the 
writers submit that if Fort Peck and Ottawa sands exhibit a strain rate effect 
of the sort contemplated by the present theory it is likely that evidence con- 
cerning its absence or its presence could not be obtained with the instrumen- 
tation installed on the rapid loading machine. 

The present theory is not strictly valid under the conditions previously 
noted. However, for creep and relaxation tests which are conducted at very 
slow rates of straining we may neglect all inertia effects and use the con- 
stitutive relationship, Eq. 7, to obtain a comparison between experiment and 
theory. Such a comparison was made for relaxation tests of Ottawa sand. In 
this case experiments indicated that such effects are present3 although, as 
noted previously, they are not very large. Whenthe appropriate value of E and 
the value of the strain rate parameter k from Eq. 23 were used in Eq. 7 and 
when computations were carried out for a 1% per sec strain rate it was found 
that there was no measurable strain rate effect. Comparison of this result 
with the experimental findings suggests that perhaps the present constitutive 
relationship does not hold for such low strain rates. On the other hand, it was 
noted previously in connection with wave propagation in Ottawa sand, that the 
present constitutive relationship seems to work quite well for strain rates 
which range from a lowof 500% per sec up to at least 15,000% per sec. There- 
fore, it may be concluded that it is not possible to use such slow relaxation 
tests in order to deduce the value of the strain rate parameter to be used in 
wave propagation calculations in accordance with the present theory. Of course 
one might seek to resolve this conflict by generalizing the constitutive relation- 
ship for the medium in order to obtain a theory which accounts for slow re- 
laxation effects as well as effects which are associated primarily with wave 
propagation. Such a possibility was considered by the MIT investigators in 
their analysis of the data.3 Their analysis was based upon linear viscoelastic- 
ity and it sought to provide for the slow and the rapid response of the material 
by means of a constitutive relationship which involved two relaxation times. 
A similar development of the present theory is certainly conceivable. However 
such a theoretical refinement will not be considered here. 


CONCLUDING REMARKS 


The results of the present study show that it is possible to formutate a 
theoretical model for impact wave propagation in columns of dry sand which, 
like the experimental observations, gives a sharp stress peak at the impact 
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16 “Stress-Strain Characteristics of Materials at High Strain Rates, Part VI. The 
Propagation of Plastic Waves in Finite Cylinders of Strain-Rate-Dependent Material,” 
by B. D. Tapley, Prepared for the Sandia Corp Contract AT(29-2)-621 by the Struct. 
Mechanics Research Lab., Univ. of Texas, Austin, August, 1960. See Also “Plastic 
Waves in Bars,” by E. A. Ripperger and H. J. Plass, Proceedings of the Symposium on 
Plasticity, Second Symposium on Naval Struct. Mechanics, Brown Univ., Providence, 
R.1., April, 1960 (in press). 
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end of a sand column immediately after it has been struck. The theory gives 
excellent agreement with experimental values of peak impact pressure for 
Ottawa and Fort Peck sands throughout the range of impact velocities from 
20 to 100 in. per sec. For Ottoawa sand this range in impact veiocity corre- 
sponds to computed impact strain rates which vary roughly from 500% to 
10,000% per sec. 

The present theory is unidimensional. It makes no appeal to any effects of 
lateral inertia. Of course the existence of such inertial effects is certainly 
not to be denied, although now their relative importance remains to be con- 
sidered. Such computations as have been carried out to show the influence of 
lateral inertia in sand24 are not particularly persuasive if agreement between 
computed and experimental results has any value. Obviously, in situations 
such as those studied here, the effect of lateral inertia constitutes a problem 
which is in neeed of further study. In this connection it is noted that B. D. 
Tapley16 has made such a study fora metal in which the constitutive relation- 
ship conforms in principle to that of Eq. 7. Although a more elaborate theory 
of this sort will eventually be most valuable for studies of wave propagation 
in sand, it would seem that the precision of the data now available tends to 
diminish its present usefulness. Nevertheless the very encouraging compari- 
sons obtained with the crude theory employed here lead to the hope that more 
precise data will become available. For then, development of a more elaborate 
phenomenological theory for the dynamical behavior of sand would have more 
than scholarly interest. 

Stress histories at the reaction end of Ottawa sand columns have been 
computed as a function of time for 14 in, and 20 in. sample lengths. In both 
cases the columns were subjected to an impact velocity of 60 in. per sec. 
Satisfactory overall agreement between experiment and theory has been found 
for both samples although the Ottawa sand experiments show a small stress 
“plateau” in the reaction stress-time curves. The possibility of obtaining 
such a plateau is also contained in the theory if a small amount of sudden 
slippage in the medium is provided for by the static stress-strain curve. 

Other experiments designed to examine strain-rate effects in soils were 
also reported by Whitman.? When these additional results were examined in 
the light of the present theory it was found that the correctness of the present 
theory was neither confirmed nor denied. 

An important difficulty in the MIT measurements was revealed by our 
comparison of the impact behavior of Fort Peck and Ottawa Sands. An inter- 
action between sand and impact stress gage was found. This interaction ap- 
pears to result from a bowing-in of the pressure sensing surface of the gage. 

A comparison of experiment with theory for histories of axial strain at 
two points in Ottawa sand columns showed rather poor agreement. In this 
comparison no allowance was made for the effect of nonzero strain-gage length 
or for slippage between the gages and the sand. 

In the present study the use of a phenomenological theory to predict the 
dynamical behavior of a dry granular medium has been proposed. This theory 
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due to Malvern, which permits the medium to undergo rather large amounts of 
irreversible straining and which specifically requires the consideration of a 
strain rate effect, shows considerable promise. However, itis certainly not a 
physical theory and no physical hypothesis concerning the origin of a strain- 
rate sensitivity in dry sand has been proposed. The absence of a suitable 

physical hypothesis regarding this strain-rate sensitivity is a serious lack. 
On the other hand, a most promising physical theory which seems to predict 
qualitatively some features of the behavior of dry sand is due to Mindlin and 
his co-workers.17 These investigators have used elastic contact theory to 
study the dynamic and static behavior of aggregates of elastic spheres having 
various packing geometries. Although the theory includes the effects of friction 
at the contact points, no specific mechanism which might lead to a strain-rate 
effect appears to have been postulated. Nevertheless, experiments on ag- 
gregates of spheres have given encouraging results. Those experiments which 
deal with stress wave propagation have been restrictedto the study of waves of 
small amplitude. This restriction is not required by our phenomenological 
theory and presumably this restriction is imposed upon the Duffy-Mindlin18,17 
theory only to secure the advantage of simplicity. In the present study there has 
been no attempt to bridge the gap between these two theoretical view points. 
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APPENDIX.—NOTATION 


The following symbols adopted for use inthis paper are listed her for con- 
venience of reference and for the aid of discussers: 


c = elastic wave speed 

E = Young’s modulus (slope of static stress-strain 
curve at zero stress) 

f, f(€) = static stress-strain curve 


17 “Mechanics of Granular Matter,” by H. Deresiewicz, “Advances in Applied Me- 
chanics,” Vol. V, Academic Press, Inc., New York, 1958, pp. 233-306. 
18 “Stress-Strain Relations and Vibrations of a Granular Medium,” by J. Duffy and 
R. D. Mindlin, Tech. Report No. 17 ONR Project NR 064-388, Contract Nonr-266(09), 
Dept. of Civ. Engrg. and Engrg. Mechanics, Columbia Univ., New York, October, 1956. 
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strain-rate parameter 

mass of gage per unit area 
compressive stress 

longitudinal particle displacement 
gage spring stiffness 

time 

particle velocity 

impact gauge velocity 

initial impact gauge velocity 
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Lagrange coordinate along sand column 


compressive longitudinal strain 


initial density of sand 
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UNCERTAINTIES IN EARTH DAM DESIGN 


By John M. Bird,! M. ASCE 


SYNOPSIS 


Some of the uncertainties that enter into embankment design and the meth- 
ods used in the designs by the Corpsof Engineers, U. S. Army, at the Termi- 
nus and New Hogan Projects to limit their effects are presented herein. Cer- 
tain areas in the design of high earth dams in which the engineer should take 
a conservative course are also examined. 


INTRODUCTION 


The increasing trend toward higher earth dams shows the need for a re- 
examining of the uncertainties that are inherent in earth-dam design. While 
it is always prudent for the engineer to take a critical look at questionable 
aspects of earth-dam design, the need has never been greater than today 
(1961) when dams in excess of 300 ft high are contemplated, designed, and 
built. The high dams, particularly those of unprecedented height, warrant 
extraordinary consideration since the designer is forced to work with many 
factors that either border on, or lie beyond, the range of his experience. It 


should be evident that extrapolation of behavior characteristics in this range - 


is fraught with danger.’ Every condition that the designer has previously con- 
sidered in establishing margins of safety for the lower dams, the exactness 
with which the loadings are known, the accuracy with which the resulting 
stresses and strains can be computed, the reliability of the materials, and the 


Note.—Discussion open until November 1, 1961. To extend the closing date one month, 
a written request must be filed with the Executive Secretary, ASCE. This paper is part 
of the copyrighted Journal of the Soil Mechanics Division, Proceedings of the American 
Society of Civil Engineers, Vol. 87, No. SM 3, June, 1961. 

1 Chf., Soils Sec., US Army Engr. Div., South Pacific, Corps of Engrs. 
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seriousness of failure if it occurs, must be examined with greater care in the 
case of the high dam, Not only will the safety of the structure depend on this 
careful study but also its economics. 

Examine some of the question marks in earth-dam design. Uncertainties 
common to all designs are found in the methods of analysis. To illustrate, 
proper evaluation of the stability of an earth structure under conditions of 
reservoir filling or drawdown demands that the ratesof saturation and drain- 
age of the embankment be determined. An exact solution to this problem has 
not been developed. An approach based on successive approximations —— 
progressive flow nets has been suggested by Harry Cedegren, M. ASCE. 
Further examples of inexactness in the design procedures are found in the 
slope stability analysis for earthquake loading and in the determination of 
settlement rates in soils which exhibit considerable secondary compression. 
Rational approaches to both of these problems are lacking. Inherent weak- 
nesses also exist in the soil sampling and testing procedures. Despite the 
advances which have been made in equipment development and techniques, 
sample disturbance and forces of unknown magnitude and effect continue to 
influence the testing results. Generally speaking, the uncertainties in the 
methods of analysis and in soil sampling and testing are beyond the control 
of the designer. There are, however, questionable elements peculiar to an 
individual design which are well within his control, 

This paper highlights the methods used in dealing with uncertainties of this 
type in the designs of two earth dams which form part of the Corps of Engi- 
neers’ extensive flood control program in California. One of these dams is the 
Terminus Dam, a 250-ft structure presently (1961) under construction on the 
Kaweah River approximately 50 miles southeast of the city of Fresno, Calif. 
The other is the New Hogan Dam, a200-ft structure to be located on the Cala- 
veras River approximately 30 miles northeast of Stockton, Calif. Construction 
of the New Hogan Dam is scheduled to commence in 1961. Although these 
structures are not to be classed among the high earth dams, the procedures 
used to remove certain questionable aspects in their designs are applicable 
to the higher dams. 

This paper presents background data on the Terminus and New Hogan 
Projects plus a description, in sequence, ofthe various stages in the develop- 
ment of the final embankment section. 


TERMINUS DAM 


Project Plan.—Terminus dam is a rolled earthfill structure which will have 
a maximum height of approximately 250 ft. An auxiliary dam is required to 
contain the reservoir at a saddle traversed by an existing state highway. The 
auxiliary dam is a rolled earthfill structure with a maximum height of ap- 
proximately 130 ft. The spillway, locatedinthe left abutment of the main dam, 
is unlined and consists of an approach channel, an uncontrolled broad-crested 
concrete weir sill, and an exit channel which will discharge spillway flows 
into an existing draw. Irrigation and flood control releases are provided by an 
outlet works consisting of a circular concrete-lined tunnel constructed in the 
left abutment. Releases are controlled by three hydraulic slide service gates 


2 Discussion by H. R. Cedegren of “Investigation of Drainage Rates Affecting Sta- 
bility of Earth Dams,” by F. H. Kellog, Transactions, ASCE, Vol. 13, 1948. 
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housed in a control tower located at the upstream portal of the tunnel. The 
layout of the main and auxiliary dams, the spillway, and outlet works is shown 
on Fig. 1. Pertinent data on the project are summarized in Table 1. 

The Terminus Reservoir will provide the maximum practicable amount of 
conservation storage without impairment of its primary flood control functions. 
As a result, high reservoir stages of long duration are expected. Based on 
hypothetical reservoir operation for observed hydrographs, it is reasonable 
to anticipate that heads of 150 ft would be sustained for periods of at least 2 


TABLE 2,—PERMEABILITY DATA FOR EMBANKMENT MATERIALS 


Mechanical Analysis Dry Weight,} Percent of] Permea- 
pe 


Average 
Grading 


Fine 
Grading 


Composite 
of Coarser 
Materials 


months, and heads of 180 ft, for periods of 2 weeks. High rates of drawdown 
are also to be expected for instances in which rapid releases must be made to 
provide space for high prospective snow-melt runoff and during flood flows 
in excess of the reservoir design flood. Such drawdown rates could be as rapid 
as 10 ft per day or more. 

Preliminary Embankment Design.—The central core type of section was 
selected for the preliminary studies. Recent streambed deposits, which cover 
the river valley both upstream and downstream ofthe site, were proposed for 


| Feature 
| percentage (ibs/cu.ft.)] AASHO |feet per day 
aan (1) (2) (3) (4) (5) (6) (7) 
| (a) Embankment Core 
i Fine 0 32 68 121 100 0.03 
| Grading 117 96 0,14 
i 115 92 0.37 : 
0 40 60 115 98 0.0001 
; 111 93 0,001 
| 105 89 0.02 
| 0 64 36 126 100 0.04 
122 96 0.40 
118 92 1.40 
i 0 57 43 126 100 0,03 
i 121 96 0,24 
117 92 0.82 
0 63 37 128 100 0.03 
125 96 0.30 
120 92 1,26 
Coarse 0 78 22 130 100 1,08 
Grading 126 96 1,62 
20 92 0.6 
(b) Embankment Shells 
i 17 14 120 100 0.13 “i 
108 85 
. 24 13 121 100 0.24 
44 3 132 Maximum 3.1 
117 Loose 61.0 
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use in the outer shells, and deeply weathered terrace deposits, located within 
short haul distance of the site, were selected for the core. Conventional 
laboratory studies indicated that these materials possessed suitable shear 
and consolidation characteristics for use inthe embankment. The permeability 
tests, on the other hand, did not yield as conclusive results. These tests in- 
dicated a wide range of values which reflect the fact that a considerable 
variation in soil grading exists in both the streambed and the terrace deposits. 
These grading ranges are shown in Fig. 2 and the permeability test results 
are given in Table 2. If it may be correctly assumed that the tests conducted 
on the average-graded materials from each of these borrow sources are ap- 
plicable to design, the streambed materials would be accepted as being suffi- 
ciently pervious to provide effective control of both drawdown and through 
seepage. Also, the terrace deposits would function adequately in the core to 
reduce the quantity of through seepage. 

The preliminary embankment section, shown on Fig. 3, was established 
on this basis. A narrow core has been used. This core, which has a width 
approximately equal to one-fourth the maximum head to which it would be 
subjected, was chosen to provide the best adjustment in the quantity and cost 
of the impervious soils available. The thickness of the slope protection is 
indicated as a minimum value. This provision was made to accommodate all 
usable rock from the spillway and outlet works excavations in these zones. 

Analysis of the preliminary embankment section by the critical arc and 
wedge methods® indicated adequate factors of safety for a variety of loading 
conditions. Included among these were the condition existing immediately 
following construction and prior to saturation, the condition of steady seepage 
under a sustained reservoir level, and conditions resulting from partial sub- 
mergence of the slope with a fluctuating pool and earthquake. Analysis of the 
rapid drawdown case was not made because, under the previously mentioned 
assumption, the material in the embankment shell would be capable of draining 
as rapidly as the pool would be drawn down and, hence, unbalanced hydro- 
static forces would not develop. 

Review of Preliminary Design.—The preliminary embankment section was 
examined objectively to determine the strong and weak points of the design. 
Each aspect of the design was reviewed inthe light of the criteria which must 
be met for the structure to be safe and stable. Consideration was given first 
to the embankment foundation. There appeared to be no uncertainty in this 
area, The embankment would rest directly on highly competent bedrock in 
both abutments and on approximately 20 ft of strong streambed material which 
overlie the bedrock in the river channel. These materials would be capable of 
resisting the stress imposed by the embankment. The possibility that detri- 
mental seepage flow through the bedrock would occur was dismissed on the 
basis that the proposed grout curtain, which would extend to a depth of 50 ft 
into the bedrock for the full length of the dam, would provide the necessary 
control. An examination of the embankment section from a functional stand- 
point revealed some uncertainty as to the stability of the slopes under various 
reservoir conditions. Despite the fact that the shell and core materials have 
high shearing strengths, the angles of internal friction as determined by con- 
‘solidated drained shear tests being 35° and 29°, respectively, a critical area 
in the design was found to center around the proper functioning of these mate- 
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rials in controlling seepage forces resulting from a drawdown, or a long re- 
tention period, of the reservoir. This uncertainty as to the operation of the 
shell and core materials stemmed from the recognition that nonuniformities 
existed within the borrow areas. The wide variations in gradation and the 
erratic distribution of the borrow soils are examined in the following sections. 

Heterogeneity of Pervious Borrow Soils.—The non-uniformity of the soils 
in the pervious borrow area is well illustrated by an examination of the grad- 
ings of soils obtained from six adjacent test holes selected at random from 
the more than ninety test holes used to explore the area. These gradings, and 
the grading band for the pervious borrow area, are shown in Table 3. A review 
of the data shows that a wide variation exists in the gradings of the soils both 
within individual test holes and in adjacent test holes. For example, the grad- 
ing of the soils from test hole 4B-156 practically envelop the entire range for 
the pervious borrow area. In contrast, for test hole 4B-157, the soils in the 
upper 16 ft are on the fine side of the grading band and those in the lower 8 ft 
approach the average grading for the borrow area. The soils in the upper 
layers of some of the test holes are relatively fine materials and, for other 
holes, they are of average or coarse grading. This heterogeneity exists 
throughout the entire pervious borrow area. It is to be noted from the per- 
meability data, shown in Table 2, that the finer fraction of the pervious bor- 
row materials would not be free-draining and therefore would not be suitable 
for use in the outer embankment zones. 

It is significant that the gradation of the finer fractions of the pervious 
borrow materials overlaps the gradation of the coarser fractions of the im- 
pervious borrow. To illustrate, Fig. 2 shows that the gradation curve for the 
10% finer material of the pervious borrow closely agrees with the 90% finer 
curve of the impervious borrow. The average gradation of the impervious 
borrow falls between the 10% finer and the finest gradation curves of the per- 
vious borrow. Considering these facts, and recognizing the erratic distribution 
in the borrow, there is a distinct possibility that sizable zones within the em- 
bankment shells could be constructed with materials of similar gradation and 
permeability to those placed in the core, The resulting discontinuities and 
stratification within the shells would tend to invalidate the assumptions of 
free-drainage of the shells, of adequate relative permeability between shell 
and core, and of adequate control of the line of seepage in the downstream 
shell, Although it must be emphasized that the probability of obtaining these 
discontinuities within the embankment shells may be somewhat remote, it is 
important to recognize that such a condition could result and that there is un- 
certainty that the structure would be built as designed. 

Heterogeneity of Impervious Borrow Soils.—The gradation curves for the 
core materials, shown on Fig. 2, indicate that there is a wide variation in the 
percentage of fines. As inthe instance of the pervious borrow, there is erratic 
distribution of the fine and coarse materials inthe terrace deposits. This non- 
uniformity in the borrow area could result in large portions of the narrow 
core being constructed with the coarser materials. The resulting stratifica- 
tion would tend to render the core ineffective for seepage control and could 
lead to a serious stability problem. 

Consideration of Field Methods to Minimize non-Uniformities.—It is rec- 
ognized that the problem of non-uniformities and discontinuities within zones 
of an embankment is common to many earth-dam projects. On some of these 
projects, various construction procedures have been used to minimize the 
variations in soil characteristics within a given zone. For example, many 
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specifications for earthwork contain a provision that the more impervious 
materials shall be placed toward the core and the more pervious materials 
toward the outer slopes. This procedure for controlling the placement of em- 
bankment materials contains a number of weaknesses inherent in the method 
itself. To illustrate, successful placement by this procedure depends on the 
ability of the field inspector to make the proper identification of the more 
pervious and the more impervious soils. In most instances this type of visual 
classification is undependable and impracticable. Another weakness of the 
procedure is its dependence on the favorable distribution of the materials in 
the borrow area. With the heterogeneous deposition of fine and coarse mate- 


TABLE 3,—PERVIOUS BORROW MATERIAL GRADATIONS 


Depth, in 
feet 


Gradation in Percentage Finer Than 


No. 200 


Test Hole 


NOON 


GRADING 


Finest 


Pervious 
Borrow 
Area No, 1 


10% Finer 
Average 


rials which exists in the borrow areas on this project, it is highly question- 
able that efforts to zone the materials by this method would prove effective. 
Another method, which has been used with some success on projects where a 
non-uniformity of the borrow materials has existed, is to blend the materials 
either at the source of excavation or on the embankment. This procedure was 
used most effectively during the construction of the Corps of Engineers’ 
Painted Rock Dam on the Gila River in Arizona. The success of this method 
requires development of a specific plan for the mixing operation. The de- 
velopment of such a plan is generally feasible only in instances where the 
boundaries of the various strata or zones of materials are clearly defined 
and uniform, 

Although conditions may be favorable for the use of the procedures out- 
lined herein, there are certain risks attendant to this dependence on con- 
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struction procedures to make up for deficiencies in the natural distribution 
of borrow materials. This does not imply that the construction operations 
would be faulty, but on the contrary, emphasizes that perhaps too much may 
be asked of the field men andthe construction procedure. The consequences of 
inability of these methods to produce the desired result should always be kept 
in mind. The risks of failure of such procedures may be warranted for the low 
embankment which would be exposed to high reservoir levels for only short 
periods, but would be indefensible for a dam such as the Terminus Dam. For 
the high dam which is to be subjected to severe reservoir conditions, the 
designer should take a conservative stand on the use of uncertain procedures. 
The design should not depend on such procedures but rather should be as 
nearly “fool-proof” as possible within the limits of good engineering practice. 

Design Modifications.—The modifications of the preliminary embankment 
section reflect tue prudent use of conservatism. The first design change in- 
volved replacement of the streambed materials in the outer portion of the up- 
stream shell with rock fill. This revision removed the uncertainty that cen- 
tered about construction of afree-draining shell with the streambed materials. 
The dimensions of this rock zone were limited to that required for stability of 
the slope against rapid drawdown. In establishing the zone thickness, it was 
assumed that the rock fill would drain as rapidly as the pool would be drawn 
down, and that porewater pressures would develop inthe remaining portion of 
the shell. It was recognized that some selectionof the rock to be placed in the 
outer zone might be required to assure free-draining characteristics. A review 
of the geologic data indicated that the major portion of the spillway excavation 
would be in fresh, hard, andonly moderately fractured rock. Except for a zone 
of fault gouge and breccia, most of the rock from the excavation would be suit- 
able as pervious fill. 

The next modification was to increase the thickness of the core to approxi- 
mately twice that prrviously used. The increased core thickness would reduce 
the possibility of developing stratification for the full core width. It was also 
felt that a more conservative ratio of core width to maximum pool depth would 
be desirable for a dam 250-ft high. It was realized, however, that increasing 
the core width would merely minimize, but not preclude, the possibility that 
seepage might outcrop at various points along the downstream slope of the 
core. 

On the basis that the random streambed materials would be used in the 
downstream shell, a number of alternatives were considered in the modifica- 
tion of the downstream shell to provide stability under steady seepage condi- 
tions. One method provided no positive control of the position of the saturation 
line in the shell and involved flattening the downstream slope as required for 
stability with a high seepage line outcropping into the random rock zone on the 
slope. A second method provided a rock toe section to control the position of 
the seepage line. Under this plan, the rock toe would have a minimum height 
of at least 100 ft above the foundation and would have an inside slope no steeper 
than 1 on 1.25 sloping upstream. 

Due to limitations on the quantity of available rock, adoption of this section 
would probably involve a reduction in the width of the rock zone in the upstream 
shell with a resulting decrease in the factor of safety of the upstream section. 
The final alternative provided positive control of the through seepage by means 
of a vertical drain adjacent to the downstream transition zone and a horizontal 
drain across the valley floor which would discharge into a nominal rock toe. 
A cost analysis of these alternative methods indicated that the embankment 


ai 
ad 


46 June, 1961 


Radius Cf Are 


Parallel To Redius 


Water Surface 
Center Of Gravity Of Moist 
And Sat. Soils 
Triel Feilure Are 


Point Of Application Of Tangential And 
Effective wW*t. (Moist + Buoyed) Normal Forces 


Tote! Wt. (Moist + Set.) 


Tengential Force 
Normal Force 


Resultant Earthquok 
Aad Tote! We. 


Earthquoke Force 


QOS (Moist + Set. wt.) 
SUBMERGED SLOPE WITH EARTHQUAKE Py! 
METHOD OF GRAPHICAL ANALYSIS 


CONDITIONS 


UPSTREAM -EARTHQUAKE WITH 
PARTIALLY SUBMERGED SLOPE 


MAIN 
MAXIMUM 


FIG. 4.—-STABILITY ANALYSIS 


SM 3 
| 
Critical Poo! E6200 4 ol 
1.609.07 
gon? Perviods Zone | 
Inactive Pool El. $70.0 
sf Core 
. 


SM 3 EARTH DAMS 


Rodius Of 


Center Of Gravity 


Triel Failure Arc 


Resultant Earthquoke 
And Total Wi. 


Eorthquoke Force 
Q05 (Moist wt.) Forces 


Point Of Application 


GRAVITY WITH EARTHQUAKE 


Parallel To Radius 


ofr 


Tongetial And Nermal 


crest E/,750.0 METHOD OF GRAPHICAL ANALYSIS 


3'Random Rock 
Zon 


8’ Vertical Drain 
8‘ Transition Zone a £1 565.0 
ous Rendom Hock 
— 8 Rock Toilwotes 
Ground Line Pervious Foundetion 
| Weathered Rock 
CONDITIONS 
— 
SECTION 


SUMMARY AT TERMINUS DAM 


47 


Q 
i 
Tongeatial Force 
Hormel 
Force 
= 
i 
Wi 


June, 1961 


section with the vertical and horizontal drain would be the most expensive; 
approximately 7% more than the rock toe section and 2.5% more than the 
flattened downstream slope section. However, it was decided that the merits 
of the positive drainage provisions more than outweighed the slight cost dif- 
ferential. As a result, that section was adopted as the final embankment sec- 
tion. The various alternative sections considered and the final design section 
are shown on Fig. 3. 

Stability of Final Embankment Section.—Stability analyses of the final em- 
bankment section indicated a safe, adequate structure. Part of a typical anal- 
ysis of the embankment slopes is shownin Fig. 4. A summary of the minimum 
safety factors obtained these analyses are givenin Table 4. The safety factors 
for the adopted embankment section are slightly lower than the recommended 
minimums for the drawdown, partial slope submergence, and steady seepage 
conditions. Acceptance of these lower safety factors is evidence that conser- 
vatism in the development of the section through the use of materials of cer- 
tain and known characteristics, and the use of positive seepage control meas- 
ures, permits some latitude in the margin of safety even in instances of a high 
dam. 

Construction Notes.—Construction of Terminus Dam was started in Feb- 
ruary 1959 and is scheduled for completion in May 1962. A view of the 
partially completed embankment is shown on Fig. 5. The photograph was taken 
from the left abutment in June 1960. As the effectiveness of the final design 
does not depend on positive control of gradation, the restrictions placed on the 
various components of the embankment are relatively minor. The specification 
provisions governing the embankment materials are outlined as follows: 


Material 


Random Rock To be obtained from required spillway excava- 
tion. Shall contain not more than 20%, by 
weight, of material passing No. 4 sieve. 


Restrictions 


Embankment Pervious Zone To be a sandy gravel obtained from specified 
borrow areas. Shall contain no material 
greater than 12 in. in size and not more than 
12%, by weight, shall pass No. 200 sieve. 


Embankment Transition Tobe obtained by selective borrowing in speci- 
fied borrow areas. Shall contain no material 
greater than 3 in. in size and not more than 


12% fines. 

Embankment Drain To be plant processed material meeting a 
specific gradation. 

Embankment Core To be silty or clayey material obtained from 
specified borrow areas. No grading limits 
specified. 


Construction of the embankment is proceeding (in 1961) according to plan. 
Although no unusual problems have been encountered, there have been some 
occurrences that show the wisdom of the conservatism used in the design. As 
anticipated, the pervious zones of the embankment are being constructed with 
materials having a wide range of grading. This grading range is shown on 
Fig. 6. 
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Fairly sizable sections of the pervious embankment zones are being formed 
of materials of considerably finer grading than the estimated average. This 
situation has developed through the use of materials from numerous sand bars 
which had been deposited in the upstream borrow area during high water 
stages of the preceding winter flood season. The embankment core material, 
on the other hand, has been found to be of more consistent gradation than had 
been anticipated. One development which has caused concernis the breakdown 
of the rock from the spillway excavation intofiner particles than had been ex- 
pected from observations of the in-situ rock and rock cores made during the 
design stage. As a result, some selection of the excavated rock for the rock 


TABLE 4,—MINIMUM SAFETY FACTORS 


Minimum Factor of Safety 


Upstream Sudden Drawdown from gross 
pool 

Upstream Partial submerged slope 

Upstream Earthquake 


Downstream | Steady seepage 
Downstream | Earthquake 


zones in the embankment has been necessary to insure a clean, free-draining 
material. An occurrence of this type helps to emphasize the need for conser- 
vatism in design, particularly in the instance of the high dams. 


NEW HOGAN DAM 


Project Plan.—The original plan for the New Hogan Project provided for 
raising the existing Hogan Dam, a variable radiusarch dam with gravity con- 
crete abutments, to create a reservoir of 237,000 acre-ft gross capacity. The 
plan was later modified to provide for the construction of an earth dam to 
create a reservoir of 325,000 acre-ft. The existing dam, which is located 
about 600 ft upstream of the proposed dam site, will have its gates removed 
following closure of the new dam and prior to the ensuing flood season, and 
will remain in place. 

New Hogan Dam will be a rolled-fill structure with a maximum height of 
200 ft. Four dikes, with a maximum height of 18 ft, will be required to con- 
tain the reservoir at three saddles north of the dam and one saddle to the 
south. Each of the dikes will have a homogeneous section, with a crest eleva- 
tion equal to that of the main dam. The spillway will be located in the first 
saddle south of the dam, and will consist of a gated ogee weir with a paved 
discharge chute terminated by a flip bucket. There will be three tainter gates 
at the weir. Flood control and irrigation releases will be provided by a four- 
barrel rectangular conduit constructed along the toe of the left abutment. 
Releases will be controlled by four hydraulically operated slide service gates 
located in a buried control chamber approximately 140 ft upstream of the 
dam axis. Access to the control chamber will be by means of a reinforced 
concrete gallery from the downstream toe. The layout of the main dam and 
dikes, the spillway, and the outlet works is shown on Fig. 7. Pertinent data 
on the project are summarized in Table 5. 
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The New Hogan Reservoir, in contrast to the Terminus Reservoir, will 
provide specific storage space for conservation of irrigation water in addi- 
tion to a defined reservation for flood control storage. As in the case of the 
Terminus Reservoir, however, encroachment on the flood control reservation 


FIG. 5.—CONSTRUCTION OF TERMINUS DAM 
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levels to the conservation pool are anticipated. These drawdowns could be as 
rapid as 5 ft per day to 10 ft per day. 

Preliminary Embankment Design.—The shallow deposition of recent al- 
luvium at the dam site precluded the use of a central core type of section 
with sand and gravel shells. Therefore, consideration was given to the use of 
a homogeneous type section composed entirely of the terrace gravels which 
cap the hills bordering the northwest portion of the reservoir area. Based 
on visual examination of these deposits, it appeared that the gravels were 
available in sufficient quantity for the dam and would be suitable for use in 
a homogeneous embankment section. The gradation range of the terrace 
gravels is shown on Fig. 8. Based on the high gravel content, it was assumed 
that the material would have high shearing strength. It was also felt that the 
percentage of fines would be sufficient to control seepage losses through the 
embankment to tolerable limits. Considering these factors, it was decided 
that an economical embankment section could be developed using only the 
terrace gravels. 

For purposes of the preliminary design, the terrace gravels were as- 
signed an angle of internal friction of 32°. The coefficient of permeability 
was considered to be less than 1 ft per day. The foundation materials were 
not assigned specific shear strengths for design as the weathered, but highly 
competent, bedrock that was considered to be much stronger than the em- 
bankment materials. The preliminary embankment section is shown on Fig. 9. 
Seepage control within the embankment is provided by vertical and horizontal 
drains. To reduce underseepage through the weathered rock foundation, a 
grout curtain is provided to a minimum depth of 75 ft. The critical arc and 
infinite slope methods were used in establishing the permissable slopes. As 
in the design of the Terminus embankment, the structure was analyzed for 
the end-of-construction condition, the earthquake condition, and for various 
reservoir conditions including the steady seepage and rapid drawdown case. 
The results of the stability analyses of the preliminary embankment section 
indicated that the safety factors obtained for each loading condition considered 
exceeded the accepted minimum values recommended for an adequate design. 

Investigation of Terrace Gravels.—To evaluate the adequacy of the pre- 
liminary embankment section, laboratory studies were conducted on the 
terrace gravels. These investigations consisted of routine classification, 
compaction studies, and consolidation, shear, and permeability tests. A 
critical review of the test data provided sound evidence that the initial evalua- 
tion of the terrace gravels had been incorrect. This review will be sum- 
marized in the following paragraphs. 

Soft Particles.—An examination of the soil grains revealed that throughout 
the entire range of particle sizes, approximately 30% of the grains consisted 
of soft particles. These particles were found to contain clay minerals with 
approximately one-third of the minerals being identified as bentonite. 

Particle Breakdown.—The compaction tests indicated that the soft particles 
would break down under dynamic loading. The results of one series of com- 
paction tests are shown on Fig. 8. The difference between the maximum 
densities obtained for the minus No. 4 and minus 1 1/2-in. materials is con- 
siderably less than that usually found with such a wide difference in grading. 
Normally, soils containing from 30% to 50% gravel will develop maximum 
densities from 10 pcf to 15 pef higher than a soil composed entirely of par- 
ticles finer than the No. 4 sieve. It was concluded that the failure of the 
coarser sample to attain a high density was due primarily to a breakdown of 
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TABLE 5,—PERTINENT DATA NEW HOGAN PROJECT 


Reservoir 


Elevation 
Spillway design flood pool $ 

Area 
Spillway design flood pool .........eeeeeeee .... 4,670 acres 

Storage capacity 
Standard project flood pool ° coe 
Spillway design flood pool . ° 


.. 15,060 acre-ft 
. . 325,000 acre-ft 
. - 325,000 acre-ft 
. . 357,000 acre-ft 


..., 
Spillway (gated, ogee) 


left abutment 
‘ 
ou yas 
12,000 cu yds 
3 42ft-by-36ft 
oe 


. 89,000 cu yds 
° ee . 22,100 cu yds 
. 44.5 ft 
Reinforced concrete conduits through. ee 
dam number and size ........ . . 4 - 5ft-8in,-by-10ft 
Gates, number and size..... ese 
Service, (hydraulic slide)..... . 4 - 5ft-8in.-by-10ft 
Emergency, (hydraulic slide) . . . . 4 - 5ft-8in,-by-10ft 
Discharge capacity (full gate opening) 
at inactive pool, elev 586 feet... . 7,450 cfs 
at gross pool elev 713 feet...... else dns «be . . 14,600 cfs 
Hydroelectric power facilities - None 


Control structure (tentative) - Reinforced concrete control chamber located approxi- 
mately 140 feet upstream of dam axis. Access by reinforced concrete gallery from down- 


stream portal of conduit. 

nee ; Main dam (rock & earth fill) 

Freeboard, above spillway 

design flood pool . . 

Length of crest .... . 1,960 ft 

Width of crest... 20 ft 

. 16.0 ft 

. 133,000 cu yds 
2,465,000 cu yds 
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. Eee Reservoir length at gross pool 
Type of Location ..... 
Crest length ........ 
Crest elevation ...... 
Total excavation...... 
Total concrete....... 
Sas Gates (Tainter) number & 
oy Head at gross pool.... 
Head at spillway...... ! 
Maximum capacity.... 
Outlet works i 
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the gravel-sized particles. It is also to be noted from the grain size distribu- 
tion curves that appreciable breakdown of the soil particles occurred in the 
compaction process, The percentage of fines is shownto have increased from 
3% to 12% after compaction, 

Low Shear Strength.—The triaxial shear tests revealed that the shearing 
characteristics of the terrace gravels were being controlled by the clay 
matrix rather than by the gravel fraction. The initial shear tests were of the 
consolidated-drained and consolidated undrained type and were performed on 
samples from which the plus 3-in. material hadbeen scalped, The test speci- 
mens were 12 in. in diameter and 24in. high and had been molded at optimum 
moisture to approximately 95% of the maximum density obtained for minus 
1 1/2-in material. The results of these tests indicated shearing characteristics 
which differed considerably from those usually associated with gravelly soils. 
The consolidated-drained tests yielded values of ¢ as lowas 18° and cohesions 
ranging from 0.16 to 0.35 tons per sqft. The consolidated-undrained tests gave 
@values of 11° and 12° and cohesions of 0.11 tons per sq ft and 0.22 tons per 
sq ft. Values of these magnitudes compare favorably to those which would be 
expected for fine-grained impervious soils. Fig. 10 indicates the results of a 
typical triaxial test. 

Consolidation During Shear.—Fig. 10, the specimen decreased in volume 
during shear despite a high initial density. This action is contrary to that 
expected of dense gravelly soils which normally expand on application of 
shearing stress, and, hence, build-up high shearing resistance due to the 
increased intergranular pressure. Other triaxial tests, conducted to determine 
the influence of particle size on the shearing strength, also showed that there 
was a definite decrease in volume during shear, particularly in the tests on 
the coarser material. A special series of unconsolidated-undrainedtests were 
conducted on minus 2-in. and onminus No, 4 materials, The 12-in triaxial ap- 
paratus was used for the coarser material and the 2.8-in. apparatus for the 
finer material. All of the specimens were tested at their molding moisture to 
simulate the condition of the embankment material immediately following con- 
struction and prior to any reservoir ponding. Although these tests indicated no 
unusual behavior as far as the effect of variations in molding moisture was 
concerned, one unusual relationship was observed. The strengths obtained 
from the tests on the minus No, 4 fraction were noted to be higher than those 
obtained on the minus 2-in. fraction. This provided further evidence that the 
coarse particles in the terrace gravels are breaking down during the shear- 
ing process. 

Revised Design Values.—On the basis of the laboratory studies, revised 
design values were selected to represent the various physical characteristics 
of the gravels after their placement in the dam: 


Unit Dry Weight, d 110 pef 

Coefficient of Permeability, K. 0.01 ft per day 

Angle of Internal Friction, ¢ 11°- 22°(Consolidated-Drained) 
Angle of Internal Friction, ¢ 11°- 15°(Consolidated-Undrained) 
Angle of Internal Friction, 11°- 15°(Unconsolidated-Undrained) 
Cohesion, C 0 (For all drainage conditions) 


The shearing properties have been expressed intermsof a range of values 
for the angle of internal friction and that cohesion has been neglected, The 
lower limit shown for the angle of internal friction represents the minimum 
value of the angle of internal friction that can reasonably be expected. The 
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upper limit is the average value obtained from the test data. The use of a 
range of design values permits the designer to judge whether adquate margins 
of safety are provided for the most adverse strength conditions as well as 
the average. This practice is particularly advisable in instances in which the 
embankment material may be altered by chemical or mechanical means to the 
extent that it will lose strength, decrease inpermeability, or consolidate less 
rapidly than at the time it was first placed in the structure. Cohesion has 
been neglected on the basis that it is not considered to be a unique property of 
the soil, but rather a function of the prestress history. Under this concept, 
cohesion is of transitory value and would diminish as the soil expands. The 
minimum value of % for the consolidated-drained and unconsolidated-undrained 
tests has been arbitrarily selected as 11° although the minimum values ob- 
tained from the laboratory tests were 18°and 15°, respectively. This conserva- 
tism in the selection of design values reflects the concern felt over the in- 
fluence of continued particle breakdown on the shearing strength and con- 
solidation characteristics of the terrace gravels. 

Development of the Final Section.—The evidence that the coarse fraction of 
the terrace gravels was not sufficiently durable to control the shearing 
strength, as had been assumed in the preliminary design, necessitated radical 
modification of the embankment section. Consideration was first given to main- 
taining the homogeneous type of section. However, the use of the revised 
design values required such drastic flattening of the embankment slopes, as 
shown in Fig. 9, that this type of section was immediately abandoned on the 
basis of economics alone. Furthermore, in view of the questionable behavior 
of the terrace gravels in the laboratory tests, it was decided that the basic 
prerequisite for the modified section would be to limit the use of these gravels 
to the central portion of the embankment where they would have little influence 
on the stability of the section. 

Considering the shallow deptns of the streambed alluvium at the site, it 
was concluded that the use of quarry rock was the only practicable solution. 
Rock suitable for inclusion in the dam as rolled rock fill was available in 
adequate quantity at the site selected as the source of riprap for protection 
of the outer slopes of the homogeneous section, The rock from this source is 
an interbedded meta-tuff and meta-agglomerate with infrequent beds of meta- 
sandstone. The grains of the meta-sandstone are flattened parallel to the 
plane of schistosity and have a leaf-like, micaceous appearance. The rock is 
moderately to slightly fractured below the weathered zone. Based on cores 
obtained from eleven exploratory holes, it was estimated that the probable 
average size of the quarried rock, as controlled by fracturing, would be 
somewhat less than 1 ft. It is anticipated that the rock will break into roughly 
cubic blocks in the less-than-1 ft size range and in tabular or slab-like frag- 
ments in the larger size ranges. Laboratory tests on the rock indicate that 
it is of high quality. For design purposes, it was assumed that the rock would 
have an angle of internal friction of at least 45°. 

Various sections utilizing rock fill in the outer zones were considered. 
These sections are shown in Fig. 9. As the quarry rock was estimated to be 
more expensive than the terrace gravels, the first of these rock-fill sections 
made maximum use of terrace gravels. This section had 1 on 2 outer slopes 
for the rock shells and a symmetrical core section with 1 on 0.75 side slopes. 
Although stability analyses of the section indicated adequate factors of safety 
for all loading conditions, the section was rejected on the basis that alteration 
of the core material with time could result in the development of excessive 
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pore-water pressures during future reservoir fillings and drawdowns, As there 
would be no way to predict this pore pressure build-up, it was concluded that 
the size of the core should be reducedto the minimum practicable width. This 
resulted in the development of the section with a 1 on 1.75 and 1 on 2.25 up- 
stream slope anda 1on1.5 downstream slope for the rock shells, and a central 
core with a 1 0n0.75 upstream side slope and a vertical downstream boundary. 
Stability analyses of this section, using the conventional wedge and critical 
circle methods of analysis, indicated adequate factors of safety for each load- 
ing condition considered. 

A summary of these analyses are shown on Figs. 11 and 12, Limiting fac- 
tors of safety are presented in Table 6. Additional soil characteristics are 
presented in Table 7 and 8. In these analyses, the unconsolidated-undrained 
shear strength, obtained from the tests conducted on specimens at molding 
moisture, was used to check the embankment slopes for the end-of-construction 
condition. The consolidated-undrained shear strength was used inthe analyses 
of the upstream slope for the sudden drawdown condition, and the consolidated- 
drained strength in the analysis of the drownstream slope for the steady seep- 
age condition. 

Final Refinements.—Although the conventional analyses indicated an ade- 
quate margin of safety, it was considered thatfurther conservatism was war- 
ranted in view of the height of the dam and the long duration of high reservoir 
stages. The first refinement involved flattening the outer slope of the down- 
stream rock shell to 1 on 1.75. This modification provides added stability in 
resisting movement of the downstream rock shell under full water load. Move- 
ments of this type, while not well understood, can be significant, and could 
cause cracking of the rock spall andtransition layers adjacent to the core and 
might ultimately lead to disruption of the core itself. Settlement records on 
existing rockfill dams, particularly on those with earth cores, are limited. 
Data presented by4,5 J. D. Galloway, J. D. Justin, Julien Hinds, Hon. M. 
ASCE, and William P. Creager, and compiledinan ASCE publication® provide 
some indication of the deformations that could be expected. To illustrate, 
vertical movement of a compacted rock shell, such as has been adopted in 
this design, could be of the order of 1% to 2% of the height. For the New 
Hogan dam, this would represent a movement of 2 ft to 4 ft. The horizontal 
movement under full water load could range from a fraction of 1 ft to several 
feet depending on such factors as the orientation of the earth core. It is ap- 
parent that these movements are more significant in the high dam since the 
water load which produces the horizontal movement varies with the square 
of the height and the vertical pressures which tend to crush the bearing points 
of the rocks and displace them are large. A second consideration involved in 
the modification was the possible influence of the flat, mica-like grains of the 
meta-sandstone on the shearing strength of the rock mass, The flatter slope 
provides an additional margin of safety against the development of potential 
sliding planes in the rock shells which could result through placement of thin, 
continuous layers of the mica-like par*icles. In this connection, the design 


4 “The Design of Rockfill Dams,” by J. D. Galloway, Transactions, ASCE, Vol. 104, 
1939. 

5 “Engineering for Dams,” by J. D. Justin, Julian Hinds, and William P. Creager, 
John Wiley and Sons, Inc., New York, 1945. : 

6 “Symposium on Rockfill Dams,” Transactions, ASCE, Vol. 125, Part 3, 1960. 


a 
m 


June, 1961 
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TABLE 7,—SOIL CHARACTERISTICS* 
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TABLE 8,—EMBANKMENT CHARACTERISTICS® 


Moist Bouyed | Saturated] Permea- 
Material |weight, in| weight, in| weight, in| weight, in} bility 
pef 


(2) (3) (4) (5) 
Rock 120 120 77 140 _ 40° to 45° 
Impervious} 110 128 68 131 0.01 |11° to 22°} 11° to 22° 


a Fig, 12 


Drained | Undrained 
Strength 


value of the angle of internal frictionfor the rock fill has been reduced to 40°. 

A second refinement involved revision of the downstream slope of the core. 
The vertical boundary results in an incomplete, or partial, wedge of rock to 
the downstream which is more susceptible to movement than a wedge with a 
sloping upstream boundary. Thus, the downstream boundary of the core was 
provided with a slight slope toward the upstream of 1 on 0.15. In accord with 
this change, the outer downstream slope was steepened to 1 on 1.6 to avoid the 
use of additional rock in the section. 

The final modification was the adoption of auniform slope of 1 on 2 for the 
outer upstream slope. This change results in an increase in the stability of 
the upper portion of the slope, principally during reservoir drawndown, and 
reduces the length of the outlet works conduit. The final embankment section 
is shown on Fig. 13. 
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Curvature of the Axis.—To provide an additional marginof safety against the 
development of cracks in the core and downstream shell, the alinement of the 
dam was modified to provide an arch upstream, of 3,000-ft radius in lieu of the 
8,000-ft radius used in the preliminary plan. The combination of the shorter 
radius of the arch and the upstream sloping of the core provides a more rigid 
support for the questionable core material to bear against. 


CONC LUSIONS 


The experiences in the development of the final sections for the Terminus 
and New Hogan Dams point out anumber of areas in which the engineer should 
be particularly alert to the possibility of introducing uncertainty into the 
design. 

Uniformity of Borrow Materials.—The designer should not place too much 
dependence on the uniformity of the borrow material. It is more often the 
exception, rather than the rule, that borrow areas for a particulr embankment 
zone will consistently yield soil of uniform gradation and plasticity character- 
istics. Lack of uniformity in these respects is almost certain to result in 
variations in the physical properties of the material. Where borrow areas 
are encountered in which no distinct pattern of distribution can be established, 
it is essential that the properties of eachof the individual soil types be deter- 
mined. No single design value would be appropriate for material from such a 
borrow area. In this case, a range of values is required to properly show the 
potential variations in the soil properties. 

In the instance of the high earth dam, where there is no room for uncer- 
tainty, it is essential that the value selectedfor use in the particular analysis 
be that which would yield the most conservative design. To illustrate, the 
design value to be used in determining the extent to which an embankment 
would become saturated during the rise of a reservoir pool should be based 
on the more permeable materials from the borrow area. In like fashion, for 
the determination of the rate and extent to which the embankment would drain 
following a rapid drawdown, the design value should be based on the least 
permeable material. Dependence on construction procedures to produce 
blending of the various soil types into ahomogeneous material is fraught with 
risks and, in general, should be avoided in high dam construction. In instances 
where the blending of various borrow components are used, it is essential 
that this be accomplished by routine construction procedures rather than by 
complex methods which are apt to break down. 

Questionable Embankment Materials.—The designer should be especially 
critical of the use of materials which appear to be in any way questionable. 
If there is indication that the soil may be altered by chemical or mechanical 
action to the extent that its shearing strength, consolidation, and permeability 
characteristics will be adversely affected, such material should be so posi- 
tioned within the embankment that it will not govern the stability of the section. 

Narrow Central Core.—The width of the central core is particularly critical 
in the high dams, particularly those with reservoir pools which would be 
maintained for long periods of time. While itis acceptable practice in earth 
dam design to provide a core width of 25% of the effective differential head 
acting against it, use of this criterion for the high earth dam may not be on 
the conservative side. It is considered better practice to use as wide a core 
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FIG. 13.—FINAL EMBANKMENT 


as is practicable commensurate with slope stability requirements. This pro- 
cedure provides added margin of safety against breaching of the core for its 
full width. Disruption of the core can result from a number of causes, As 
cited in the design of the rock fill section for the New Hogan Dam, vertical 
and horizontal movement of the rock shells can cause cracking of the pro- 
tective filter and bedding zones adjacent to the core. As a result, the de- 
formations of the core under full water load may not be adequately resisted. 
Breaching of the core can also be caused by differential settlements due to 
abrupt discontinuities in the abutment slopes or to widely different consolida- 
tion rates of contiguous foundation or embankment soils in either their natural 
state or on saturation. 

To illustrate this effect, an instance is cited in which an earth dam, 200 ft 
high with a narrow central core of fine-grained material, developed outcrops 
of cloudy seepage on the downstream slope. Seepage control is provided by a 
horizontal drainage blanket located downstream of the core and a rock toe 
drain. The seepage outcrop was noted during a period when the reservoir 
level was being maintained at approximately 15 ft below spillway level, equiva- 
lent to a differential head of approximately 150 ft. An examination disclosed 
that a crack had developed through the core in the abutment section of the dam. 
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The orientation of the crack was approximately at right angles to the dam 
axis. During the investigat'on, it was found that the embankment had been 
constructed directly on a haul road which had been benched into the abutment. 
The steep excavation for the haul road had provided an abrupt discontinuity 
in slope. The fact that the crack through the core closely followed the aline- 
ment of the haul road makes the cause of the failure apparent. This case em- 
phasizes the need for adequate scaling of abutment slopes to insure that the 
embankment will be constructed on uniform slopes. The case also points up 
the need for positive seepage control downstream of the core. Such positive 
control can be effectively obtained by use ofa vertical drain extending for the 
full height of the core. 

In final summation, the designers of high earth dams should strive to make 
the embankment section as simple and as “fool-proof” as possible. 
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~ GEOPHYSICS EFFICIENT IN EXPLORING THE SUBSURFACE 


By R. Woodward Moore! 


SYNOPSIS 


Two geophysical methods, refraction seismic and electrical resistivity, are 
briefly described, their countrywide use mentioned, and numerous examples 
given of their application to a variety of subsurface problems pertinent to high- 
way design and maintenance. Among these are studies of landslide conditions 
and tests made over water-covered areas. 


INTRODUCTION 


Much of the trouble and expense attendant to unforeseen and undesirable sub- 
surface conditions found after construction has begun could be avoidedor allevi- 
ated in some measure by a more thorough exploration of the subsurface prior 
to design and construction. Some of the reasons for failure to obtain more com- 
plete information on the geologic conditions existing in a proposed building 
area, ata prospective bridge site, or along the right-of-way of a projected road- 
way, are as follows: (1) insufficient funds, (2) lack of qualified personnel, 
and (3) insufficient time due to the need forearly construction of a given pro- 
ject. Occasionally, poor judgment with regard to the need for more complete 
subsurface information may be involved. 

Geophysics can effect much savings, both in time and cost, when used along 
with the usual direct exploration techniques. Ordinarily, when this is done, 


Note.—Discussion open until November, 1961. To extend the closing date one month, 
a written request must be filed with the Executive Secretary, ASCE. This paper is part 
of the copyrighted Journal of the Soil Mechanics Division, Proceedings of the American 
Society of Civil Engineers, Vol. 87, No. SM 3, June, 1961. 

1 Head, Geophysical Explorations Group, Soils, Foundations, and Flexible Pavement 
Branch, Div. of Physical Research, Bur. of Public Roads, Washington, D. C. 
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SUBSURFACE EXPLORING 


the geophysical tests have their greatest usefulness when used prior to any di- 
rect exploration programs utilizing borings, test pits, and core drilling. Much 
progress has been made during the past 10 yr in the application of two rela- 
tively simple but dependable geophysical tests, the refraction seismic test and 
the earth-resistivity test, to such problems as (a) determining foundation con- 
ditions for bridges, large buildings, dams and high earth fills, (b) obtaining 
information leading to the proper slope design on grading projects for highways, 
(c) studying existing and potential landslide conditions, and (d) locating 
sources of construction materials and determining the extent of the deposits. 
Fig. 1 shows the present status (1961) of geophysics in one civil engineering 
field, highway construction, in the United States. Hawaii also has geophysical 
equipment. Four states, Dlinois, Indiana, Kentucky, and Massachusetts, do not 
have geophysical equipment in their highway departments but have had a limited 
experience with the use of geophysical tests in cooperation with other state 
agencies suchas a university or the geological survey. About 1/5 of the states 
shown using geophysical methods make use of the seismic test. The Bureau 
of Public Roads uses earth-resistivity equipment in park and forest road sur- 
veys in 5 of its 11 regions and has used the equipment in the Philippines and 
Central America. Many other federal, state, and private agencies now make 
use of one or both of the simple geophysical tests, applying them toa wide vari- 
ety of field problems ranging from the location of uranium to finding adequate 
water supplies in areas of water scarcity. 


THEORETICAL ASPECTS 


With regard to the fundamental theory involved in the two test procedures, 
many of the theoretical assumptions made when evolving the formulas used in 
computing the geophysical data for both types of test become rather nebulous 
when making the actual field tests. The surface and near-surface conditions 
existing at many of the locations where geophysical tests are made vary rather 
widely from the smooth, level surfaces and interfaces between successive for- 
mation layers assumedintheory. Nor does the perfectly homogeneous charac- 
ter assumed for the various materials encountered in each layer materialize 
to the degree pictured by those responsible for the theoretical computations. 
For these reasons much of the geophysical data obtained have been analyzed 
by use of empirical rather than theoretical procedures. 

In the following only brief mention will be made of the theoretical aspects 
in the two test procedures. One should refer to the voluminous reference ma- 
terial available through the many published reports for detailed discussions of 
the theory involved. (A partial list of references is given in the Appendix.) 

Seismic Test.—Inthe refraction seismic test, as illustrated in Fig. 2, small 
man-made earthquakes are set up by exploding small charges of explosives in 
shallow shot holes and the time required for the vibration or sound wave to 
travel through the various materials is measured. The more dense layers 
such as solid rock will possess a high wave velocity of from 10,000 fps to 
18,000 fps, whereas the unconsolidated soil layers will have low wave velocities 
of about 1/10 of these values. For suchformations as hardpan, cemented sands 
and gravels, shales and weathered rock, intermediate wave velocities of 4,000 
fps to 8,000 fps will likely prevail. Little trouble is experienced in indicating 
the presence and depth at which the denser layers may be found. 
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As shown in Fig. 2, the wave travels through the low velocity soil from shot 
point to detector for the shorter distances. As the shots are placed farther 
away from the detector spread, the fastest wave, the one that must be used, 
passes down through the low velocity layer, rushes ahead in the high velocity 
layer, andis refracted to the surface through the soil layerto reach the detec- 
tors. Since a portion of the wave path is common for both detectors, the short 
additional interval of time required for the wave to reach the second detector 
is a measure of the velocity in the rock layer. The data are plotted in a time- 
distance graph such as that shown in Fig. 3 for Stevens Canyon Project, Mt. 
Rainier National Park, Washington. The velocities for each of the three layers 
involved were computed from the reciprocals of the slopes of each straight- 
line segment of the graph. 

Depths tothe layers involved canbe computed by a simple formula involving 
the intercept on the time axis and the wave velocity computed for each layer. 
Thus, for the surface soil layer in Fig. 3 


1.75 x 700 
Hy = —3x200 * 
the intercept on the time axis being 1.75 divided by 200 and the computed wave 
velocity being 700 fps. The factor 2 in the denominator is required since the 
wave travels through the surface soil twice, once beneath the shot point and 
again beneath the detector. Since the explosives producing the seismic waves 
are usually buried 3.0 ft below the surface, 1/2 of this depth, or 1.5 ft, must 
be addedas shown in the preceding equation and in Fig. 3, to give the total thick- 
ness of the surface layer. Accurate depth determinations to formation changes 
have been made many times in the seismic work of the Bureau, the Corps of 
Engineers, the United States Army, and by others (4),(6),(10),(11).2 

Earth-Resistivity Test.—The earth-resistivity test makes use of electric 
currents passedthrough the ground to produce a potential drop proportional to 
the resistance present inthe various subsurface materials from which the re- 
sistivity or apparent specific resistance of the materials is computed. The 
resistivity is usually recorded in ohms per centimeter cube (ohm-cms) or the 
resistance existing between the opposite faces of a centimeter cube of the ma- 
terial. The flow of current is electrolytic in nature, making use of the moisture 
in the soils and rocks. This moisture, together with dissolved impurities there- 
in, controls to some degree the relative values of the measured resistivity. 
Moist clays will possess low resistivities ranging from a few hundred to a few 
thousand ohm-cms. Clean sand and gravel, a much sought construction ma- 
terial, will generally possess much higher resistivities, in some instances as 
high as 1,000,000 ohm-cms. Dense, solid rock such as granite will have re- 
sistivities of from 1,000,000 to 4,000,000 ohm-cms, whereas some rock for- 
mations of volcanic origin may be quite hard but havea resistivity of only 5,000 
to6,000 ohm-cms. This emphasizes the fact that there is no fixed fundamental 
relation between the engineering properties of most soils and rocks and their 
measured resistivity values such that tabulated values may be used in analyz- 
ing the subsurface conditions from place to place. An intelligent and reliable 
interpretation of earth-resistivity data can be made only after adequate trial 
or calibration tests have been made over known rock and soil formations in a 
given area. Moving to another ridge in a mountain range may involve rock or 


2 Numerals in parentheses refer to corresponding references in the Appendix. 
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soilof a different type, making necessary new trial tests to prove their resis- 
tivity characteristics. This constitutes no problem, however, since the trial 
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tests are rapidly made, requiring only 10 min to 15 min time. When obtaining 
data for use in slope design in mountainous areas where deep cuts are needed, 
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the writer obtains many calibration curves as closeas possible to the right-of- 
way of the proposed roadway in order to make certain the materials involved 
in the trial tests are also present in the roadway prism. 

Fig. 4 shows a view of the milliammeter-potentiometer type of resistivity 
equipment used by the Bureau of Public Roads anda majority of the states. 
The Bureau prefers this type of apparatus making use of direct current, which 
is alternateda single time manually, in order to recognize the presence of stray 
potentials naturaltothe ground and perhaps make some use of them in special 
analyses of a particular field problem. The Bureau recommends the highest 
degree of sensitivity inthe potential measuring circuit that can be attained and 
maintained under the conditions which the field tests will involve. Perhaps this 
has been more necessary in the Bureau’s experience since work is done inall 
parts of the country and on a large variety of field conditions not necessarily 
present in the same degree in any one state or region. Much of the resistivity 
equipment now in use was built by state highway department personnel, or others 
retained by the state, making use of plans and specifications obtained from the 
Bureau. The galvanometer used in the potential circuit has a sensitivity of 
0.125 mu a per mm on its own scale but is sufficiently rugged to perform well 
in a portable appartus. 

Resistivity-Depth Test Described.—When making a resistivity-depth test 
the 4 electrodes used are placed in a line spaced a distance A apart, herein- 
after referred to as the electrode spacing. As illustrated in Fig. 5, it is as- 
sumed that an equipotential bowl or hemisphere of radius A is set up around 
each of the two outside current electrodes. Every point on the surface of a 
hemisphere is at the same potential due to the current flowing from Cj to C9. 
The potential drop or potential difference between the two hemispheres is meas- 
ured by a separate potential measuring circuit using P, and Pg placed on the 
ground where the equipotential hemispheres intersect the ground surface. The 
measured values of E, the potential drop, and I, the current flowing in the cir- 
cuit, are inserted in a simple formula, 


andthe resistivity per centimeter cube is obtained. This resistivity is assumed 
toapply to a zone of soilA centimeters in depth, thus giving rise to the empiri- 
cal relation used that the electrode spacing is equal to the depth involved. By 
expanding the electrode spacing A to successively greater values the resis- 
tivity obtained corresponds to increasingly greater depths. The resistivity is 
plotted against electrode spacing or depth as shown in Fig. 6. The field data 
as obtainedare plotted in the dashed-line relation andare replottedusing cumu- 
lative values of resistivity, keepingthe increment of electrode spacing or depth 
constant for the entire cumulative curve, that is, by using 3-ft increments of 
depth tothe full depth of 30 ft. Straight lines are drawn through as many points 
as possible on the cumulative curve. The first two lines intersect at a point 
which has been found to approximate closely the depth of the overburden to a 
second geologic layer such as the solid rock shown in Fig. 6 This empirical 
procedure was reported by the writer (7) in 1944, The trend in the dashed- 
line curve to much higher values of resistivity at a depth of 12 ft is signifi- 
cant in establishing the validity of the intersection obtained at 13.6 ft in the 
cumulative curve as representing the thickness of the surface layer. 

By keeping the 4 electrodes at a constant spacing, say 20 ft, a constant- 
depth resistivity traverse canbe made showing the effect of changes in the ma- 
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terial within 20 ft of the surface. The resistivity readings are made at frequent 
and regular intervals along the line of tests, usually at intervals equal to the 
electrode spacing chosen. The 4 electrodes are moved forward as a group 
after each reading is obtained. This tesi is used tolocate and outline the later- 
al limits of deposits of sand and gravel, ridges of solid rock, etc., that may be 
present within the depth involved in the test. 


APPLICATIONS OF GEOPHYSICAL TESTS 


Perhaps the most efficient way to emphasize the usefulness of these eco- 
nomical test procedures is to list and describe briefly a few of the many suc- 
cessful applications made of them by the Bureau of Public Roads in connec- 
tion withthe park and forest road construction program and while demonstrat- 
ing the test procedures in 43 states and in 4 Central American countries. 

Investigations Made with Seismograph.—The refraction seismic test has been 
applied to a study of the subsurface conditions involving a number of problems 
associated with highway construction and other civil engineering works. Usually 
acrew of 4 to 5men will perform efficiently, even in the less accessible areas, 
when using any one of several types of portable, multi-channel seismographs 
now available. The single-channel devices utilizing hammer blows in lieu of 
explosives for producing the seismic waves require only 2 or 3 men but may 
have some limitations for routine use on alltypes of seismic explorations. Ex- 
amples of the application of the seismic test to five problems are given sub- 
sequently. The equipment used was a small, portable, three-channel seismo- 
graph, developed by the Bureau in 1934, which makes use of simple carbon- 
button type detectors without employing electronic amplification in the circuit. 

1. Slope Design.—The seismic data of Fig. 3 were obtained some years ago 
while makinga survey for slope design purposes in Mt. Rainier National Park. 
The comparatively low velocity of 3,125 fps obtainedfor the second layer indi- 
cated at this test site is good evidence of the absence of solid rock within 44.0 
ft of the surface, the depth at which solid rock should be found, and suggests 
the needfor earth slopes inthe soil and boulder mixture comprising the layer. 

Fig. 7 shows data for a seismic test made on the Virginia Bank of the Potomac 
River in asurvey of a 2-1/2-mile long grading project only recently constructed. 
Here, again, the relatively low velocity of 2,855 fps is ample rroof of the ab- 
sence of hard rock within the first 21.5 ft below the surface and indicates the 
need for earth slopes to that depth. The material below this depth should be 
weathered to some appreciable degree to produce a velocity of only 6,890 fps. 
The planned grade depth at this location was only 19 ft. 

' 2, Tunnel Sites.—Seismic tests at proposed tunnel sites have utility in pro- 
viding the character of the materials at the portals and in determining the char- 
acter of materialabove the invert. Usually its usefulness ends at the first high 
velocity layer reached by the test, since the presence of less dense materials 
with their characteristically lower velocity cannot be indicated by this test be- 
neath the higher velocity layer. As a matter of fact, the seismic test loses 
much its effectiveness in areas where stratification results in alternate zones 
of nard and softer materials, that is, sandstones and shales, gravels and silts, 
etc., lying above the elevation to which the test must be made. Fig. 8 shows 
data obtainedin western North Carolina while investigating one of 5 tunnel sites 
on a 3-1/2-mile section ofthe Blue Ridge Parkway in 1950. The relatively low 
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velocity of 7,920 fps obtained for the final layer indicated is suggestive of weath- 
ered or badly fractured and jointed rock at a depth of 43 ft, neither of which 
would be particularly good for tunnel work. As a matter of fact, this particu- 
lar tunnel was eliminated and a relocation to permit an open cut of about 108 
ft in depth was made when initial excavation at one of the proposed portals 
failed to find material suitable for tunneling. 


Vez 14900 


SECOND 


200 


" 350 

" 350’ NORTH 


1.5 = 5 SHOT DEPTH 
7.3' SURFACE SOIL 


1:8. x1300 
400 


WW 
> 
a 
WwW 
> 
= 
= 
uw 
WwW 
= 


_ 65% 4950 _ 


400 80.5 CLAY OR GLACIAL TILL __ 


87.8' TO ROCK 


| 
200 250 300 
SHOOTING DISTANCE - FEET 


FIG. 10.—TIME-DISTANCE GRAPH FOR SEISMIC TESTS AT A PROPOSED DAM 
SITE IN NEW ENGLAND 


3. Bridge Foundations.—The seismic test is peculiarly adaptable to a study 
of bridge foundations, since any layer possessing high wave velocities can be 
expected to have excellent supporting power. Fig. 9 shows data obtained in 
northern California while investigating the foundations for a large bridge to 
cross an arm of Lake Britton on the Shasta-Lassen Highway. Based on other 
experience, the high velocity layer found at a depth of 53.6 ft should have ample 
density andstrengthto support the proposed structure. The overlying material 
witha wave velocity of only 2,380 fps would be in question pending careful test- 
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ing inthe laboratory to determine if it was sufficiently stable to support spread 
footings in lieu of piling driven to refusal on the deeper layer. 

4. Dam Sites.—Foundations for large dams can be investigated by the seis- 
mic testto advantage. This is particularly true when several possible locations 
are involvedand when preliminary or reconnaissance survey using the rapidly 
made seismic test are possible to eliminate all save the best location for more 
detailed investigation by supplemental core drilling. Fig. 10 contains data for 
one of many seismic tests made in studies of proposed dam sites in various 
parts of the country. This test, made by Bureau personnel when demonstra- 
ting the usefulness of the test procedure to the United States Army, Corps of 
Engineers, in 1937, indicated the presence of a dense, hard rock at a depth of 
about 88 ft. Subsequently, the Corps of Engineers made good use of the seis- 
mic test in investigations of several hundred prospective dam sites in connec- 
tion with their flood control program. 

5. Quarry Sites.—The seismic test is useful in the location and evaluation 
ofa potential quarry site. In addition to obtaining a good estimate of the depth 
ofthe overburden andamount of stripping required, this test can produce some 
indication of the character of the rock involved. The higher wave velocities 
will indicate the more dense nature of the rock. However, when relating ve- 
locity todensity of the rock care must be taken to obtain a true representative 
velocity by making reverse profiles and, perhaps, right angle tests or profiles 
over the same test areasto prevent obtaining erroneous values of velocity due 
toa sloping interface between the rock and the surface materials. Due consid- 
eration also must be given to known geologic conditions characteristic of the 
local rock, since jointing and seams in relatively dense rock might produce 
an intermediate wave velocity likely to be attributed to a much less desirable 
rock formation. Fig. 11 shows data for a seismic test made in southwest Iowa 
in connection with a search for quarry material. 

Single-Channel versus Multi-Channeled Seismic Equipment.—Trouble may 
be encountered when attempting to analyze data obtained when using recently 
developed single-channel seismic equipment, which does not provide perma- 
nently recorded time data, inlieu of the more conventional multi-channel seis- 
mographs that do make use of photographically recorded time data. There are 
conditions when the relative velocities of surface soils, intermediate layers, 
and the underlying rock formations are such that the arrival of sound waves 
from one or more layers can be masked by wave arrivals from the deeper ly- 
ing but higher velocity rock layer. Fig. 11 shows data which would have pro- 
duced a very normal-appearing, time-distance graph, based on first wave ar- 
rivals, with no hint of the presence of the 2,650 fps layer shown immediately 
above the limestone rock if a permanent time record had not been available 
for detailed study to disclose the existence of the secondary time events used 
to produce this very important segment of the time-distance graph. Actually, 
if this higher velocity layer had not been identified and used in the computations, 
the depth tothe limestone formation would have been in error by approximately 
31%. Use was made of such secondary wave-arrival times when plotting the 
time-distance graphs of Figs. 7 and 8. Fig. 12 shows a seismic record obtained 
with the Bureau’s three-channel seismograph showing first and second arrivals 
such as those referred to previously. 

Investigations Made with Earth-Resistivity Test.—In the writer’s experience 
it has been demonstrated that the earth-resistivity test possesses even greater 
utility for studies of shallow subsurface conditions over a wider range of geo- 
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logic conditions than is possible when using the seismograph. For this reason 
alarge percentage of the geophysical exploration work carried on by the Bureau 
of Public Roads in recent years has made use of the earth-resistivity test. 
Four men are normally used for efficient operation of the resistivity test, one 
to operate the instrument, one to record and plot the data obtained and two men 
(usually laborers) to move the four electrodes. Considerable work has been 
done by the writer with only two helpers and at times with a single helper: Ex- 
amples of the application of the resistivity test to 8 field problems are pre- 
sented. 

1. Slope Design.—With the trend toward use of deep cuts and high fills in 
building our modern parkways and expressways, there is even greater need for 
acareful consideration of the design of cut bank slopes. Often it is possible to use 
composite slopes or benching to assure the stability of high cvt banks. The 
earth-resistivity test has been used successfully bythe Bureau: Public Roads 
and others in connection with slope design problems in many p’ ts of the coun- 
try. Two graphs have been prepared from data obtained for tests made locally 
in the Washington, D. C., area in connection with slope design on a section of 
the George Washington Memorial Parkway betweer the Spout Run Parkway and 
Chain Bridge on the Virginia bank of the Potomac River. Fig. 13 shows a re- 
sistivity-depth curve obtainedat Station 129+50, westbound roadway, where the 
planned grade was 25 ft belowthe surface. The resistivity data indicate changes 
at 9.2 ft and at 33.0 ft, the latter change being to solid rock. The material to 
be found between depths of 9.2 ft and 33.0 ft was interpreted as weathered rock 
in the report furnished by the design engineer prior to construction. Fig. 14(a) 
shows several calibration curves obtained over weathered material exposed in 
the river bank along the right-of-way (lower part of graph) and over solid rock 
(upper three curves). In Fig. 14(b), 11 resistivity curves are shown that were 
rather typical of about 90%of the 82 tests made on the 2-1/2-mile long project 
while furnishing training in the use of the resistivity test to several groups of 
engineers. Obviously, the materials involvedonthis project must be weathered 
for the most part since the shape and trends in 90%of the curves obtained ap- 
proximate those appearing in the curves plotted for the calibration test data 
obtained over such subsurface materials. One ortwotest pits dugin cuts where 
the resistivity tests indicated 100,000 chm-cms and 200,000 ohm-cms resist- 
ivity, respectively, would provide adequate proof of the existence of the weath- 
ered material. This section of the parkway has been constructed. The ma- 
terials above grade were weathered as indicated by the resistivity test, requir- 
ing slopes of 1 to 1 as designed to assure stable cut slopes. 

2. Bridge Foundations.—The earth-resistivity test has proved itself as a 
very fast, dependable means of obtaining the preliminary information on sub- 
surface foundation conditions required by those concerned with bridge design. 
A1/2-hr test may provide adequate evidence that footings placed on solid rock 
for a bridge pier will not be practical by indicating a depth to solid foundation 
material in excess of that considered economical. In some instances pilings 
will be used in design and their approximate length determined from the results 
of the geophysical tests. Fig. 15 shows data for subsurface conditions found 
at 1 of the 4 typical bridge sites selected on the Virginia side of the Potomac 
River at Washington, D. C., where many field tests have been made in recent 
years while demonstrating this test procedure to hundreds of engineers from 
all over this country as well as from 20 foreig: countries. This test, made 
just upstream from the old 14th Street Bridge, indicates a continuous, high re- 
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sistivity formation at a depth of 102.2 ft. The average of 15 tests made here 
was 100.8 ft to rock, the maximum depth 107.5 ft, and the minimum 95 ft. These 
tests, made at various times during the seasons and by engineers from many 
parts of the world, have indicated the subsurface conditions as accurately as 
would be requiredfor design. Piling with lengths averaging 100 ft would be re- 
quired, driven to refusal on the rock formation. 

Trends inthe measured values of resistivity, such as that appearing at 99 ft 
in Fig. 15 (middle curve), assume their true significance as evidence of the 
presence of ledge rock when one considers that the higher resistivity associated 
with each 9-ft increase in depth, as the test penetrates deeper into the layer, 
has overcome the effect of 99 ft of overburden to cause the curve to trend up- 
wardas shown. The middle curve of Fig. 15 is a replot of a portion of the lower 
curve using a scale large enough to emphasize the uptrend in the curve. Even 
less pronounced trends become significant in the analysis of resistivity data 
under such conditions. Solid rock was found by the drill at a depth of 110 ft at 
the location of the abutment for the new 14th Street Bridge some 800 ft to 900 
ft downstream from the location where the data for Fig. 15 were obtained. Re- 
sistivity tests made at the other3 test sites have shown evidence of solid rock 
at depths of 33 ft, 51 ft, and 75 ft at locations near the Key Bridge in George- 
town, just below the Memorial Bridge and about halfway between the Memorial 
Bridge and the 14th Street Bridge, respectively. 

The resistivity testis useful in limiting the number of drill holes at a par- 
ticular structure site, when design features dictate some direct examination 
of materials above a bedrock formation and within some stated depth below the 
rock surface. Itis also useful in locating areas of significant change, both lat- 
erally and at depth, thus permitting the placing of drill holes at strategic loca- 
tions and assuring that they be carried to sufficient depth to obtain the maxi- 
mum informationforthe investment made. Fig. 16 shows data obtained in west- 
ern North Carolina while testing the site for a large structure to carry the 
Blue Ridge Parkway over a highway and railroad passing through Balsam Gap 
some 8 miles west of Waynesville. Two resistivity crews spent about one work- 
ing day here while making 10 separate depth tests to depths averaging 120 ft 
to 130 ft. The nearly level rock surface inferred from these tests could be 
checked by use of a sounding rod, driven to refusal at a depth approximating 
the elevation of the rock surface as shown. One drill hole to study the charac- 
ter of the overlying material and obtain a core in the bedrock would prove this 
location rather well. 

3. Tests Made over Water-Covered Areas.—The fact that the earth-resist- 
ivity test can be made on water almost as easily and effectively as on land 
makes it even more valuable for investigating bridge foundations. During the 
course of demonstration work, tests have been made directly on water surfaces 
in Alabama, Louisiana, South Carolina, Alaska, and Nicaragua in Central 
America. Fig. 17 shows a view of atest in progress in the middle of the Ashley 
River at Charleston, S. C. Two large oil drums were anchored 400 ft to 500 ft 
apart and connected by a small rope floated on the water surface. Using five 
small wooden boats, the test was begun in the center and expanded in either 
direction to encompass a total depth of 120 ft. In such tests, the 4electrodes 
are floated on the water surface and spaced by marked, nonmetallic tapes in 
the same manner as tests madeon land. Fig. 18 shows the data obtained at 
this location together with drilling results shown inthe lower part of the graph. 
An 8-ft tide, together with adverse winds, kept the small boats off the water 
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for much of the 2-day period available for the tests and prevented 1 or 2 ad- 
ditional tests being made to the depth indicated in Fig. 18. The salt water pos- 
sesseda rather uniform resistivity of about 50 ohm-cms. The several changes 
indicated by the drill at this location apparently have been shown equally well 
by the rapidly made resistivity test. The increased effort required to drive 
a 4-in. diameter casing under the action of a 300-lb hammer falling 30 in. is 
shown along the lower part of the graph. The rather definite increase in the 


@CORE DATA 


ROCK Ev Ev REMARAS 

ROCK tat 
22.2 
ano Roce 
' e 


FIG. 19.—PORTION OF ROCK CONTOUR MAP, PREPARED FROM DATA OBTAINED 
BY EARTH RESISTIVITY TESTS 


measured resistivity at depths of 41 ft, 65 ft, and 99 ft appears to be closely 
identified withthe bottom of the water and the underlying stiffer layers requir- 
ing more blows per foot penetration of the casing. Makingtests in this manner 
at each pier location over a 1,200-ft water crossing such as that involved in 
this instance could disclose similar trends in the resistivity curves obtained 
ateach pier which could be provedat one or more locations, perhaps on either 
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bank, with the elimination of costly drilling at some of the pier sites. The 2 
holes drilled at the locations of the piers for the bascule spans crossing the 
ship channel cost about $10 per ft for depths of approximately 100 ft. The defi- 
nite rising trend in the resistivity curve beyond the 99-ft depth is ample evi- 
dence of a continuous layer composed of the muchstiffer material found by the 
drill at about that depth. The test results indicate this material is at least 21 
ft thick. The layer thickness could have been proved for much greater depth 
if the wind and tide had permitted. 

Some provision for control of the applied current will be required when test- 
ing on salt water surfaces. In the test under discussion the writer resorted to 
a suitable rheostat in the current circuit to keep the test current within the 
range of the apparatus used (300 ma). Various applied voltages ranging from 
1.5v to 9.0v were used to permit full use of this maximum range in current. 
This was essential because even with the fuli-scale current readings the po- 
tentials measured ranged from a value of 15.4 mv for the first 3-ft depth of 
the salt water to an average reading of only 1.79 mv for 37 separate readings 
made for depths greater than 24 ft below the water surface. This emphasizes 
the needfor using a very sensitive and dependable apparatus capable of making 
such small potential measurements whenever the necessity arises. 

When testing over water-covered areas, one of the potential hazards will 
be the local game warden who is always on the watch for those attempting to 
obtain fish illegally and will likely question those found using mysterious appa- 
ratus floated on the water surface. 

4. Resistivity Tests Applied to Building Sites.—In 1942 the Navy Depart- 
ment’s Bureau of Ships requested thata contour mapbe made of the rock sur- 
face underlying a 150-acre tract near Washington, D. C., where a facility had 
been constructed for testing ship models. Following the completion of approxi- 
mately 500 depth tests and about 10-1/2 miles of constant-depth resistivity 
traverse, a rockcontour map was made,a portion of whichis shown in Fig. 19. 
Two years later the existing building was extended to a point 1,800 ft farther 
tothe east, the area involved being about 120 ft by 1,800 ft in plan. The amount 
of soil stripped from the area to expose the rock surface, comprising some 
100,000 cu yd, differed by less than 6%from the amount of stripping indicated 
by the rock contour map and surface contour map made 2 yr earlier. 

Fig. 20 shows a curve typical of the constant-depth resistivity traverse data 
obtained in this study and illustrates the possible use of such data to eliminate 
many unnecessary direct borings or drill holes in an area where the resistivity 
data suggest a uniform depth of overburden as indicated by the flat lying portion 
of the traverse between stations 0 + 00 and 8 + 00. Any obviously anomalous 
localized conditions which should receive further and detailed exploration are 
alsolocatedin a striking manner as shown by the two resistivity highs appear- 
ing inthe right-hand portion of Fig. 20(c). Thetraverse involved a 20-ft depth 
along a 2,000-ft line. The numbers shown underlinedare results of resistivity 
depth tests for depth of overburden. Curves for two such tests are shown in 
Figs. 20(a) and (b). 

5. Locating Deposits of Construction Materials.—Large quantities of granu- 
lar material are sought and used in subbase and base courses, pavements, and 
other highway structures. The earth-resistivity test has a particularly good 
application to the location and delineation of the boundaries of buried deposits 
of such material having no surface indication of their presence. The constant- 
depth traverse is used to great advantage to locate the higher resistivity zones 
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usually associated with sands and gravels wherever they are found. Fig. 21 
shows data for 2 constant-depth resistivity traverses, of the many that were 
made on sections of the Baltimore-Washington Parkway when locating large 
amounts of granular material required for use in subsequent construction. Both 
tests involved a depth of 30 ft. The lower curve, by reason of the very low val- 
ues of resistivity, failed to show evidence ofany useful granular material in the 
cut area tested. The upper curve, however, showed strong evidence of an ex- 
cellent source of granular borrow material suitable for the purpose for which 
it was required. 

Data from curves of this type are plotted in a resistivity contour map such 
as that shown in Fig. 22. The zones of higher resistivity suggest that the area 
involved contains large amounts of granular material such as that needed for 
alayer of topping material tobe placed over clay soils present in the subgrade 
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FIG. 21.—CONSTANT DEPTH RESISTIVITY TRAVERSE OVER SAND AND GRAVEL 
DEPOSIT AND OVER CLAY FORMATION 


elsewhere on the project. Only a few borings are required to prove the pre- 
dominately granular character of this high resistivity material. A boring to 
sample the materials at one of the lesser resistivity contours, say the 250,000 
ohm-cms contour, disclosing material suitable for topping material, would indi- 
cate material of equal or better quality for such use present in the zones of 
higher resistivity without further borings being made. This is true because 
the higher resistivity normally results from a preponderance of material of a 
nonplastic nature devoid of deleterious amounts of clays and silts. 

Fig. 23 shows data obtained in Mississippi while demonstrating the test pro- 
cedure in the Jackson area in the fall of 1958. Due to the extreme scarcity of 
granular material inthat area, as well as other parts of the state, borings were 
being made to depths of 80 ft to 90 ft to prospect possible sources of granular 
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material associated with the Citronelle Formation from which much of the 
granular material obtained in the area is taken. Overburden depths of 30 ft to 
40 ft were not uncommon in the area tested. In this instance the resistivity 
test showed strong evidence ofthe granular materialata depth of 40.0 ft, about 
3.0 ft deeper than that shown by the borings. The test also showed definite 
trends towards higher resistivity at depths of 60.5 ft and 90.8 ft, which are 
probably associated with layers containing higher percentages of granular ma- 
terial. The change at 60.5 ft corresponds quite well with the depth of 57.0 ft 
where the boring found A-3 soil having about 50% retained on a No. 40 sieve. 
The layer immediately overlying this material was of the same type soil, A-3, 
but with only 35%retained on the No. 40 sieve. The borings were terminated 
ata depth of 87 ft, just 4 ft short of the 90.8-ft depth where the resistivity test 
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FIG. 23.—RESISTIVITY DEPTH TEST SHOWING CORRELATION WITH BORINGS 
OBTAINED OVER DEEP-LYING SAND AND GRAVEL FORMATION 


snowed a further strong increase in resistivity likely associated with a layer 
of even better granular qualities. Since the resistivity test carried to such 
depths embraces huge volumes of material, a series of tests made at regular 
intervals, perhaps comparable to those arbitrarily chosen for the borings, 
woutd provide overlapping zones of influence and produce depth data much more 
representative of the total volume than could be gotten by small diameter, iso- 
lated bore holes. An entire prospective pit could be investigated by the rapidly 
made resistivity test, evaluated, option obtained and finally proved by a few 
well placed borings at locations dictated by the evaluation of the resistivity 
survey. Samples taken from the materials found by the borings could be used 
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in a laboratory analysis of the character of the material comprising the various 
layers involved. 

Use of the resistivity test to locate quarry material is a practical applica- 
tion. Fig. 24 shows data obtained ata quarry site near Fenn, Idaho. Figs. 24(a) 
and (b) show data for two constant-depth resistivity traverses, one over an 
area underlain by good quarry material (Fig. 24(a)), the other indicating the 
presence of unsuitable material for the first 150 ft and outlining the lateral con- 
tact of this material with the lava rock being sought (Fig. 24(b)). 

The smal! resistivity contour map shown in Fig. 24(c) shows this contact 
to better advantage. The 10,000-ohm-cm contour should delineate the boundary 
conditions rather closely. Note that the early part of traverse No. 1 parallels 
an exposed face ofan existing quarry where suitable material exists to a depth 
of about 22.0 ft, or 2.0 ft deeper than the depth involved in the traverse. 

The depth curves shown in Figs. 24(d) and (e) show the way in which the 
resistivity data disclose good material (Fig. 24(d)) or unsuitable material 
(Fig. 24(e)). Two or 3 hr work was sufficient to get the information shown in 
Fig. 24. A third resistivity traverse running at right angles to traverse No. 1 
in the vicinity of station 2 + 00 would have served to prove the quarry area 
further, together with extensions of the two traverses shown. 

6. Investigating Landslides.—The growing importance of landslide studies 
warrants emphasis on the use of the resistivity test on such problems. The 
resistivity test has producedsome good results in landslide studies, by reason 
of the effect upon the measured resistivity of the excess water at the slip sur- 
face or shear zone. Althoughthetest procedure has been used in brief demon- 
stration tests on landslide problems in about 15 states, a complete and rather 
thorough survey ofa particular landslide has been made in only a few instances. 
Much work along these lines remains to be done. Fig. 25 shows one of several 
cross sections obtained by earth-resistivity tests through a slide area originat- 
ing on Ohio Route 7 between East Liverpool and Wellsville in the fall of 1951. 
The close check on the predictedcontour of the firm material based on the re- 
sistivity data is apparent from the depthstofirm material obtained by the core 
drill at the three locations shown. Much useful information of this type can be 
obtained prior to making borings, the resistivity data being used to locate the 
borings to obtain the greatest amount of useful information. 

7. Investigating Tunnel Sites.—The rapid investigation of proposed tunnel 
sites is possible with the resistivity test. Fig. 26 shows data obtained several 
years ago when testing for subsurface conditions for 5 tunnel locations and 
about 3-1/2 miles of grading onthe Blue Ridge Parkway in western North Caro- 
lina. This investigation required about 4 days. The solid-line curves are for 
calibration data obtained on a nearby project. Thetwo lower calibration curves 
were obtainedat the location of a tunnel that had collapsed during early con- 
struction operations following which a redesignforan open cut was made which 
involved a depth of approximately 100 ft. The upper calibration curve was ob- 
tained over a hard-rock cut which was selected as typical of good tunnel ma- 
terial, whereas the dotted-line curve is one of the several curves obtained at 
the 5 tunnel sites. Obviously, the location of the proposed tunnel No. 3 is ina 
ridge where poor tunnel rock will be found, equally as poor as that existing on 
the other project where the tunnel had failed. This is the same tunnel where 
the seismic data of Fig. 8 was obtained. It is of interest to relate that, follow- 
ing failure to find suitable material for tunnel construction at the proposed 
portals when construction was begun in the fall of 1957,a relocation was made 
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with plans for an open cut more than 100 ft deep. The writer was asked to make 
further resistivity tests over the relocated line to determine if the final 60 ft 
of the cut could be made using 1/4 to 1 rock slope in a composite slope design. 
The test results indicated the need for earth slopes of 1-1/4 to 1 for the entire 
depth of cut. This was based on a comparison ofthe currently obtained resist- 
ivity curves with those obtained in 1950 for tunnel No. 4. Tunnel No. 4 was 
built without use of explosives inthe soft weathered rock involved. Some blast- 
ing was done in the last 6 ft or 8 ft in the cut on the relocated line. 

8. Tests Made In Swampy Areas.—The resistivity test can be very useful 
in an investigation of swamps. The bottom formations can be located rather 
easily and rapidly in areas where it is possible to walk on the surface materi- 
als. Fig. 27 shows a resistivity party working in a cypress swamp some 14 
miles north of Mobile, Ala. Tests were made at 58 locations across a 5-mile 
section of the swamp along the centerline of a proposed 4-lane divided highway. 
These tests found firm sand as bottom materialin some areas and marine clay 
in other areas. The depth of the muck ranged from 17.0 ft to 60.0 ft in 40 of 
the tests made where muck was definitely present at the surface. The average 
depth to firm sand or clay was 34.0 ft. A few soundings made here and there 
should suffice to prove the entire 5-mile long survey. Due to commitments 
elsewhere in the state, time was available for only 1 or 2 such checks to be 
attempted during the geophysical survey. At one location a rod equipped with 
a 1-1/2-in. to 2-in. diameter sampling barrel was pushed to a depth of 32.0 ft 
where increasing frictional resistance and stiffening of the muck, perhaps, 
made further penetration slow and laborious. The resistivity test had indi- 
cated firm material at a depth of 40.0 ft. At another location sand was found 
by the sounding ata depth of 19.0 ft where the electrical test had indicated 22.5 
ft. This constitutes a fairly good check since the resistivity test has a reason- 
ably large area of influence and produces an average depth which would likely 
be checked more closely by an average of several probings made throughout 
this area. Further comments of resistivity tests made over swampy ground is 
presented elsewhere (7) ,(13),(15),(18). 

Other successful applications of the electrical test include the study of ground 
water conditions, location of suitable water supplies, explorations for mining 
operations, and studies of undergroundcorrosion problems. Nofurther consid- 
eration will be given these rather important uses, however, since the writer 
has had little first-hand experience in such explorations. 


CONCLUSIONS 


The economical andefficient use of the refraction seismic test andthe earth- 
resistivity test for obtaining information on subsurface materials and condi- 
tions has been reported on heretofore in the many published papers dealing 
with such work, only a few of which have been included in the Appendix. Al- 
though the materiai presented in the current paper shows successful results 
also, it must not be construed that these rapid test procedures are infallible. 
There will be geologic formations and combinations thereof which present only 
slight degrees of change in the velocity characteristics or resistance to flow 
ofan electric current. In instances where suchconditions exista very thorough 
and careful geophysical survey may require more frequent check tests by di- 
rect means and, in some instances, perhaps a rather complete investigation by 
borings or test pits will be necessary. Fortunately, it has been the writer’s 
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experience that these unusually difficult field conditions occur infrequently, 
the geophysical tests being applicable in a large majority of cases when they 
are applied to problems that are more or less amenable to their particular 
advantages or limitations. 

The seismic test will normally have good application to foundation studies 
for bridges, buildings, and dams, and can be used to good advantage in deter- 
mining the character of the materials present in roadway cuts and at the por- 
tals of proposed tunnel sites. This test has little value for use in locating ma- 
terials of construction such as sand and gravel and it may prove rather inef- 
fectual in a detailed study of landslide areas. Also, in areas where thin layers 
of relatively hard material, such as sandstone, limestone, or comparatively 
freshanddense volcanic rock are underlain by weathered materials or forma- 
tions of lesser density, the seismic test will be ineffective since the higher 
velocity wave inthe dense layer will reach the detectors first and tend to mask 
the effect of the arrival of the waves through the lower-velocity layers lying 
below. Further, seismic work in water-covered areas and swampy ireas will 
be more complicated when using most of the conventional multi-channeled 
seismic equipment and may become impossible when using some of the recently 
developed single-channel equipment now being offered for use in shallow sub- 
surface exploration. At the present time (1961), the Bureau makes use of the 
seismic test to supplement subsurface data obtained by use of the resistivity 
test and to prove the results obtained at 1 or 2 test locations where resistivity 
tests were made when some question exists regarding the analysis of the re- 
sistivity data. 

The resistivity test is equally as good as the seismic test for foundation 
studies when proper calibration of the test is made over known materials and 
it has obtained excellent results on many slope design problems. This test is 
not greatly hampered by thin layers of hard material overlying less dense for- 
mations, it being possible to locate a shale bed or otherless resistant material 
beneath sandstone or other materials possessing high resistivity. Its use for 
locating construction materials has been well demonstrated in the field and 
reported on inthe literature. In landslide studies, since the water associated with 
many landslide areas will likely affect the values of the measured resistivity, 
the resistivity test has promise for use in obtaining more accurate and detailed 
information on the formations or conditions associated with the slip surface. 
Also, the possibility of the water causing the sliding being responsible for 
measurable electric potentials that can be associated with a particular slip 
surface or contact between soils of differing moisture content or plasticity 
must not be ignored. 

Use of the resistivity test over water-covered areas constitutes one of its 
more valuable applications, since it is then in direct competition with relatively 
expensive direct procedures for exploring the subsurface conditions. 

The need for complete and reliable information onthe subsurface conditions 
connected with the larger and more important construction and maintenance 
problems will dictate the use of allavailable means to meet this need in order 
thata stronger, safer and more dependable structure will result ata more rea- 
sonable final cost devoid of the added expense of extra work orders, etc., at- 
tendant to unforeseen conditions found after construction is begun. 
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GLOSSARY OF TERMS AND DEFINITIONS IN SOIL MECHANICS 


Closure by the Committee on Glossary of Terms and Definitions in 
Soil Mechanics of the Soil Mechanics and Foundations Division 


Introduction.—The seven discussions have contributed to the value of the 
report, and have provided information, related opinions, and corrections that 
are of importance and benefit to the committee in its continued review and re- 
vision of the report. 

The number of discussions the report received, together with comments re- 
ceived personally by individual members, convince the committee that prepara- 
tion of a revised glossary such as this is desirable now, and on a continuing 
basis. 

Since publication of the original report, the committee has been reorganized, 
and a joint ASCE-ASTM Committee on Glossary of Terms and Definitions in 
Soil Mechanics and Foundations has been formed. In addition to commenting 
on the several discussions, the committee takes this opportunity to report on 
its subsequent activities related to terms and definitions. 

This closure has been endorsed by action of the joint ASCE-ASTM Commit- 
tee on Glossary of Terms and Definitions in Soil Mechanics and Foundations. 

Comments on Discussions.—The discussions are concerned with both gen- 
eral considerations and details. Because the glossary is a compilation, rather 
than a creative effort, it is essential to obtain such information from users of 
the nomenclature and notations. Assembly of a record of what is being used 
is a continuing process; the establishment of general preferences is evolution- 
ary. 

Barber and Focht make a number of pertinent suggestions. These, and all 
comments and recommendations made by the discussors, will be considered 
carefully inthe further work of the committee. However, the following specific 
observations may be of value at this time. Mr. Barber suggests the use of the 
No. 10 sieve size as the upper limiting size for sand. ASTM Committee D-18 
on Soils for Engineering Purposes has adopted the No. 4 sieve size for this 
upper limit. The matter will be reconsidered by the joint ASCE-ASTM com- 
mittee. Term 324, Silt, should be changed to agree with Term 325 as noted by 
Barber and Focht. 

The committee appreciates the inclusion by Jakobson of the “ Proposed List 
of Symbols” prepared by the Swedish Geotechnical Society. Close correlation 
with work being done in other countries is essential. This is being achieved 
through a committee of the International Society. Participation in the work of 
this committee will be discussed later. 

Ahlvin’s comments are timely. The joint ASCE-ASTM committee will, no 
doubt, wish to consider this problem at an early date. 


4 October 1958, Report of the Committee on Glossary of Terms and Definitions in 
Soil Mechanics of the Soil Mechanics and Foundations Division. 
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While sympathetic to Parisi’s desire for simplicity inthe glossary, the com- 
mittee is of the opinion that it is preferable to continue the use of symbols that 
have been widely accepted. The notations, in context with appropriate sub- 
scripts, are not as generally confusing, or duplicating, as is indicated by an 
alphabetized listing. 

The committee is most greatful to Legget for his careful and discerning re- 
view. It is in agreement with his “minor” comments, and with many of his 
“major” ones. They will all be carefully considered in the future deliberations 
of the committee. 

Heidema notes the absence of some terms that are associated with recent 
developments in soil mechanics. As stated in the preface to the report, the 
preliminary compilation was completed in 1951. Few additional terms were 
subsequently added during the lengthy review. The committee is cognizant of 
the fact that its compilation is incomplete and has assigned itself the task of 
corretting these omissions. 

ASTM Action.—The glossary presented in Paper No. 1826, with those cor- 
rections that were brought to the ASTM Committee’s attention by the discus- 
sions, or throughother channels, was recommended for adoptionas a standard, 
ASTM Designation D653-58T (See ASTM Standards, Part 4, Pg 1211). This 
recommendation was submitted toletter ballot by the Administrative Commit- 
tee on Standards on September 30, 1959, and passed. The following statement 
was included: “These definitions were prepared jointly by the American Soci- 
ety of Civil Engineers and the American Society for Testing Materials.” 

International Definitions and Symbols .—R. E. Fadum isthe U.S. representa- 
tive to the committee on Glossary of Terms of the International Society of Soil 
Mechanics and Foundation Engineering. A study of the proposed list of sym- 
bols appearing in the International Standards Association’s Draft Recommen- 
dation No. 276 has been undertaken by the International committee. The ASCE 
Committee on Standards and Definitions contributed a review of this listing. 

A Russian translation of Proceedings Paper 1826 was issued by the Central 
Bureau of Technical Information, Moscow in 1960. 

Committee Reorganization.—At the October, 1959, meeting of the Executive 
Committee of the Soil Mechanics and Foundations Division, the committee was 
renamed the COMMITTEE ON STANDARDS AND DEFINITIONS, whose stated 
purpose is: 


To give continuing study to definition of terms and symbols in soil me- 
chanics, and to report periodically on matters relative thereto; to rec- 
ommend actions related tothe establishment of endorsement of standards 
in the field of interest of this division. 


The following are the current members of this committee: E. D’ Appolonia, 
R. E. Fadum, C. R. Foster, J. K. Mitchell, S. R. Stearns, and G. A. Leonards, 
chairman. 

In view of the enlarged scope of the committee’s activities, and in order to 
avoid duplication of effort, it was deemed desirable to organize a joint ASCE- 
ASTM Committee on Glossary of Terms and Definitions in Soil Mechanics and 
Foundations. Suchaction was consummated by the two Societies by June, 1960. 
ASCE’s representatives on the joint committee are the members of the Com- 
mittee on Standards and Definitions, as noted previously; ASTM is represented 
by Sub-committee G-3 of Committee D-18, Soils for Engineering Purposes, 
whose members are as follows: E. A. Abden-Nur, D. M. Burmister, A. I. John- 
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son, M. D. Morris, A. A. Wagner, and R. G. Abnlvin, chairman. The joint com- 
mittee is in the process of organizing its activities, that will be directed along 
the following general lines: 


1. The correction of errors and omissions in the report. 

2. Continued re-evaluation of the report, and its discussions, leading to the 
improvement and addition of definitions wherever possible, and to the selection 
of preferred symbols wherever feasible. 

3. Preparation of periodic reports to inform the profession of progress and 
to stimulate general and specific discussion. 

4. Cooperation with the International Committee in the formulation of an 
international glossary. 


In conclusion, the committee expresses its sincere appreciation to those 
who accepted the invitation to comment on the glossary. It is hoped that dis- 
cussion will continue both before and after subsequent reports so that an up- 
to-date and complete glossary can be assembled. 
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STRUCTURE AND STRENGTH CHARACTERISTICS OF 
COMPACTED CLAYS* 


Closure by H. B. Seed and C. K. Chan 


H. B. SEED,! M. ASCE, and C. K. CHAN,2 A. M. ASCE.—The writers ap- 
preciate the additional information previded in the discussions of this paper. 
Da Silveira has suggested an interesting alternative method of presenting the 
results of compaction tests and following the changes in composition of com- 
pacted samples. He also points out the significance of volume change charac- 
teristics in determining pore pressures in soils and the importance of consid- 
ering these effects in developing a better understanding of strength character- 
istics. The writers would certainly agree with this, but believe it to be rela- 
tively incidental to the purpose of the paper that is to illustrate the influence 
of soil structure on strength characteristics. Because structure will also af- 
fect compressibility, consideration of this factor is indirectly included in the 
results. The writers doubt that Da Silveira’s computation of the volume change 
required for complete saturation in connection with the test data presented in 
Fig. 7 has any real significance in regard to the shape of the stress-strain 
curves. Stress-strain curves having an essentially similar form were obtained 
in the tests described in Fig. 8, but the samples in this case had a quite dif- 
ferent initial composition. Because the two samples compared in each of these 
figures were essentially alike in all respects except structure, it is believed 
that the differences in their characteristics can reasonably be attributed tothe 
influence of this factor. 

The writers welcome the additional data on the influence of structure on 
swelling characteristics and pore water tensions provided by Gilbert. Gilbert 
also raises the interesting question of the uniqueness of effective stress soil 
strength parameters for compacted samples in the partially saturated and sat- 
urated conditions. Unfortunately, it does not seem possible at this stage to 
resolve this question, primarily due to the great difficulty in making reliable 
measurements of negative pore water pressures in samples compacted dry of 
optimum. There is also considerable disagreement on appropriate methods of 
determining effective stresses in partially saturated soils. However, the writ- 
ers tests were conducted on samples that are either saturated or close to sat- 
uration and Gilbert is quite correct in pointing out that the data presented show 
only for samples in this condition that structure has little or no effect on soil 
strength characteristics expressed in terms of effective stresses. It.is also 
gratifying to learn that he agrees with the concept that the initial structure of 
a compacted soil is progressively destroyed as samples undergo significant 


4 October 1959, by H. B. Seed and C. K. Chan (Proc. Paper 2216). 

1 Asst. Prof. of Civ. Engrg., Engrg. Materials Lab., Univ. of Calif., Berkeley, Calif. 

2 Junior Research Engr., Inst. of Transp. and Traffic Engrg., Univ. of California, 
Berkeley, Calif. 
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deformation so that large differences in initial structure may be eliminated by 
the time failure occurs. 

The evidence provided by McRae in support of the flocculated structure of 
samples compacted dry of optimum and the dispersed structure of samples 
compacted wet of optimum is also a welcome addition to available information. 
Like McRae the writers have also noted the relatively spongy or springy prop- 
erties of samples compacted wet of optimum. However, it should be noted that 
the degree of springiness is considerably different for samples prepared wet 
of optimum by kneading and static compaction methods. Fig. 1 presents data 
on the elastic compression and rebound observed in repeated loading tests on 
samples of a silty clay prepared to identical densities and water contents by 
these methods of compaction. Both samples appeared spongy or springy in 


Sample preporedby 
Aneading compaction. 


Sample prepared by 
Static compaction . 


Resilient percent 


9 
% 


Silty clay; water content = /4.3% 
Ory density = 1/9.5/b percuft 
Contining pressure =0.25 kg per sqcm 
Repeated deviator stress = 0.7 kg per sq cm 
Duration of stress application =0./5 second 
Frequency of stress application = 30 per minute 


/ 10 100 1000 10,000 100,000 
Number of Stress Repetitions 


O./ 


FIG. 1.—EFFECT OF METHOD OF COMPACTION ON RESILIENCE CHARACTERISTICS 


comparison with samples compacted dry of optimum. It wili be seen however 
that the recoverable deformations are almost 100% greater for the samples 
prepared by kneading compaction indicating that relative springiness cannot 
readily be assessed from qualitative types of tests and that quantative deter- 
minations are required to accurately asses differences in compacted samples. 

The groupings of soils adopted by the Corps of Engineers provide a logical 
basis for differentiating between different soil types and the writers are not 
surprised at McRae’s statement that “soils in group (c) (highly swelling soils) 
are not extremely sensitive, by the CBR test, to molding water content in the 
as-molded condition.” The writers suspect that high swelling soils are more 
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difficult to disperse by compaction wet of optimum than are clay soils exhib- 
iting low expansion characteristics, thereby reducing the possibility of produc- 
ing the extremely low strength at low strains characteristic of samples with 
dispersed structures. 

The structures described by McRae are not greatly at variance with those 
prepared by the writers, and it would seem that McRae is in agreement with 
the suggested effects of shear strain by his statement “This picture would be 
modified as to degree of dispersion to account for variations in the amount of 
shear strain when using various methods of compaction.” Because this is one 
of the more important points made in the paper it is gratifying to see that in 
spite of minor differences, McRae is in agreement with the writers inthis im- 
portant concept. 

Finally it should be noted that the data in McRae’s Fig. 4 were apparently 
obtained from samples prepared by static compaction and tested in the as- 
compacted condition. The results are therefore comparable to the data pre- 
sented by the writers in Figs. 21(a) and 21(b). Furthermore, they have the 
same form as the data in these figures. It is difficult therefore to understand 
McRae’s statement that the form of the results is directly opposite to the trend 
shown for statically compacted specimens. It is believed that there must be 
some confusion between the text and the figure with reference to this point. 

The explanation of the shape of the compaction curve and the drying shrink- 
age curve given by Olsen and Scott forms as interesting addition tothe original 
paper. However, it would appear that they have unfortunately misconstrued 
some of the writers statements. For example, no attempt was made by the 
writers to explain the shape of the compaction curve; rather the paper was 
simply concerned with clarifying the type of structure associated with different 
positions of compacted samples on a typical compaction curve - a somewhat 
different problem. Nor did the writers attribute the engineering properties of 
compacted clays to double-layer phenomena. In fact, a deliberate attempt was 
made to show that factors other than double-layer effects are responsible for 
the final structure and properties of compacted soils. It would be difficult to 
point this out more forcibly than by the statement, 


“This leads to a modification of the original hypothesis for the develop- 
ment of soil structure during compaction as follows: at low water con- 
tents the high electrolyte concentration prevents the double layer from 
developing fully resulting in low inter-particle repulsion and consequent- 
ly, for most soils, a tendency for flocculation of the clay particles. An 
increase in water content causes a decrease inelectrolyte concentration, 
expansion of the double layer, and therefore increased repulsion between 
clay particles together with a tendency for higher pore-water pressures 
to develop when the tamping pressure is applied; consequently it leads 
to adecrease in shear strength. The somewhat reduced tendency to floc- 
culate but more particularly the greater shear deformations under the 
tamping foot which result from the reduction in shear strength can lead 
to an increased degree of dispersion or an increased degree of orienta- 
tion of the clay particles in the compacted soil.” 


In other words, dispersed structures at higher water contents tend to result 
primarily from the greater shear deformations resulting from the reduction 
in strength due to increased pore pressures - a quite different concept than 
the simple application of double layer effects. 
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In order to avoid confusion, it should also be pointed out that double layer 
effects cannot be completely excluded from consideration and that the data in 
the discussers’ Fig. 1 may be somewhat misleading to readers unfamiliar with 
other work on this problem. It is clearly apparent that for the tests reported 
by Olsen and Scott, changes in pore water salt concentration were ineffective 
in altering the position of the compaction curve. However, the tests were con- 
ducted on a clay made homo-ionic to calcium and, therefore, one that would 
have strong flocculation tendencies regardless of the pore water salt concen- 
tration. Under these conditions, significant structural changes would not be 
anticipated, and the results presented are simply what might have been ex- 
pected for this particular soil. However, there is abundant data from other 
investigations? to show that the addition of small amounts of dispersants can 
cause appreciable changes in the position of compaction curves for soils. It 
would indeed be dangerous to ignore these data completely and draw conclu- 
sions only from the data in the discussers’ Fig. 1. Furthermore, for reasons 
pointed out by Lambe, even when electrolyte concentration does cause changes 
in compacted density, there is no reason to expect corresponding changes in 
‘as compacted strength’ at any given water content. Thus the data in Fig. 2 
could be anticipated regardless of any variation in the position of the compac- 
tion curves in Fig. 1. 

The writers believe that the shear strength properties of clays can be sat- 
isfactorily explained either in terms of soil structure, that involves consider- 
ations of shear strain, double layer effects and other interparticle forces, or 
in terms of effective stresses and pore water pressures, but that consideration 
of all these phenomena is required for a complete analysis of the properties 
of compacted clays. 


3 “The Improvement of Soil Properties with Dispersants,” by T. W. Lambe, Journal, 
Boston Soc. of Civ. Engrs., April, 1954. (Fig. 5) 
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SOIL STRUCTURE AND THE STEP-STRAIN PHENOMENON@ 


Closure by D. H. Trollope and C. K. Chan 


D. H. TROLLOPE,37 M. ASCE, and C. K. CHAN,38 A. M. ASCE.—The writ- 
ers are indebted to D. F. Coates for his interesting contribution concerning 
experience with the Ottawa clay and look forward to reading further published 
details of this work. Coates questions the use of the term “colloidal friction” 
to describe the component of shear strength contributed by the clay fraction. 
In support of his objection he quotes the well known experimental fact that 
“witha constant void ratio it would not be expected that this ‘colloidal friction’ 
would vary with the magnitude of the normal stress.” This statement does not 
contradict the concept of colloidal friction. The writers interpretation is that, 
because the void ratio is constant, the component of effective normal stress 
that is derived from inter-particle forces is a constant and thus, the related 
part of the total shear resistance is also a constant. 

In this context the writers may well be accused of “splitting hairs.” Recent 
work at the University of Melbourne39 has confirmed, however, that the con- 
cept of colloidal friction is of considerable interest in understanding the be- 
haviour of ‘pure’ clays as distinct from composite scils. As a result of this 
work it has been suggested that the shear resistance (S) of svils may be ex- 
pressed in the form. 


in which $y is the angle of internal friction; o is the applied normal stress; U 
is the pore pressure; and p is the component of effective normal stress due to 
inter-particle forces. For sands p: is effectively zero, and the @r is then di- 
rectly related to the coefficient of micro-friction between individual particles. 
Similarly for a “pure” clay, ¢r is also directly related to the coefficient of 
micro-friction. For composite soils, however, in which both sand, silt, and 
clay fractions contribute to the strength, the measured value of $y is an em- 
pirically determined statistical mean dependent not only on the physical char- 
acteristics of the sand and clay particles but also ontheir relative proportions. 

Coates also inquires whether there is any substantiation of the hypothesis 
concerning movement of the coarse grains under shear. Fig. 34 shows a cut 
section through a cylindrical sample composed of glass marbles embedded in 


@ April 1960, by D. H. Trollope and C. K. Chan (Proc. Paper 2431). 

37 Senior Lecturer in Civ. Engrg., Univ. of Melbourne, Melbourne, Australia. 

38 Junior Research Engr., Inst. of Transp. and Traffic Engrg., Univ. of Calif., Berke- 
ley, Calif. 

39 “The Fabric of Clay in Relation to Shear Strength,” by D. H. Trollope, Proceed- 
ings, 3rd Australia-New Zealand Conf. on Soil Mechanics and Foundation Engrg. (In 
press) 
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a kaolin matrix. The specimen was prepared by placing symmetrical layers 
of marbles at equal vertical intervals in the compacted kaolin and then tested 
in atriaxial compression machine. This investigation40 showed that the move- 
ment of the marbles depended mainly on the characteristics of the kaolin. The 
molding water content of the sample in Fig. 34 was 35% compared with the 
Proctor standard optimum moisture content of 24% for the kaolin. For sam- 
ples prepared at or close to optimum, although some relative movement of the 
marbles was observed, failure occurred primarily through the kaolin matrix. 
A svu>sequent investigation?! was carried out in which the marbles were re- 
placed by a standard (Leighton Buzzard) sand (98% passing BS No. 14 sieve 


FIG. 34._MOVEMENT OF MARBLES IN KAOLIN 
AT FAILURE 


and retained on BS No. 25 sieve) and thekaolin matrix was prepared in a stand- 
ard manner.42 Triaxial samples 1 1/2 in. diameter x 3 in. long were prepared 
from slurries containing 40%, 60%, and 80% sand content on adry weight basis. 
Each sample was consolidated in the triaxial cell under 60 psi ambient stress 


40 “A Step-Strain Model Incorporating Kaolin and Single-sized Spheres,” by A. G. 
Johnson and J. D. Kaye, Student Thesis, Univ. of Melbourne, Melbourne, Australia. 

41 “Stress-Strain Behaviour of Sand-Clay Mixtures,” by G. D. Handley and R. A. G. 
Vines, Unpublished Student Thesis, Univ. of Melbourne, Melbourne, Australia. 

42 “The Measurement of Soil Properties in the Triaxial Test,” by A. W. Bishop and 
D. J. Henkel, London, Arnold, 1957. 
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and then subjected toan undrained test with pore pressure measurement. Only 
those samples with 80% sand content showed some indication of step-strain 
behaviour as determined from the incremental strain plot and even this evi- 
dence was inconclusive. Similarly, no significant conclusions can be drawn 
from the measured pore pressures in these tests. The writers recognise that 
in the absence of adequate pore pressure information the detailed evaluation of 
the step-strain mechanism is deficient. It does appear however, that to resolve 
this deficiency will entail measurement of pore pressures within the samples 
and this presents some experimental difficulties. 

Coates has criticized the criterion adopted by the writers in referring step- 
strain behaviour to the erratic nature of the incremental strain plot following 
matrix yield. As the writers see it this erratic deformation ina controlled 
rate of loading test is the direct consequence of what has been termed step- 
strain behaviour, and as such serves to identify soils that are subject to this 
phenomenon. Within our present knowledge it is not possible to specify the 
degree to which this erratic behaviour is developed. Primarily it is a function 
of internal strain velocities. 

This has been pointed out in the paper where it has been observed that rapid 
testing procedures may obscure any steps in the stress versus strain graph, 
and the only evidence that the soil is a step-strain type is obtained from the S 
shape of the curve. It is because of the high internal strain velocities developed 
in the early stages of the constant load creep tests that no steps are observed 
in the curves of Figs. 22 below 60% of the normal unconfined compression 
strength. The steps only show up in those samples that cannot reach an equi- 
librium structure and progressively deform until failure develops. 

Coates’ final question concerns the presence of bound water films and as- 
sociated hydrogen bonding of these films as a contributing factor to cohesion 
in clays. It has been noted previously in this discussion that the first named 
writer is of the opinion that what has been called cohesion in clays is best de- 
scribed as an effective stress term in a basically frictional mechanism. The 
effective stress in this term being derived from interparticle forces. This 
raises the question of the nature of “contact” between colloidal particles, or 
alternatively the structure of the “link bond” as it has been described in the 
paper. At present there are two schools of thought on the characteristics of 
the water films adjacent to these contacts. One school holds that the water 
films are virtually “solid,” the other that the water retains its fluid charac- 
teristics. Little is known of the nature of these “contacts” however, andit does 
not seem possible at present to resolve these arguments. It would appear to 
the writers, however, that a likely explanation is that the water loses its iden- 
tity, as water, at the contacts, the hydrogenions and the hydroxyls become dis- 
sociated under the high stresses and recombine in some way with the constit- 
uents of the colloid particles in a manner similar to the welding of asperity 
contacts in metal friction. Thus, hydrogen bonding may be one of a number of 
mechanisms contributing to the contact strength. Such considerations are 
merely conjectural at the present time however and further research is needed 
to clarify the position. 

The writers are grateful to B. K. Hough for his resume of some of the back- 
ground work that has contributed so much to improved understanding of soil 
structure. Owing to space limitations it was not possible to do justice to this 
background work in the original paper and therefore the engineering approach 
of dealing only qualitatively with the forces of attraction and repulsion between 
colloid particles was adopted. The writers were aware of the contributions of 
Bolt and other workers from Cornell University, Ithaca, N. Y., but it must be 
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admitted that at the time of writing they were more concerned with the inade- 
quacies of the parallel plate model in describing the nature of soil structure 
than they were in considering the exact nature of the forces in this model. It 
is in this context that the writers felt that attention had not been paid to the in- 
fluence of external stresses in the formation of soil structure. Hough is right 
in drawing attention to the fact that Bolt had considered external stresses in 
the behaviour of the parallel plate model. Subsequent developments have indi- 
cated that Bolt’s work, and that of others who have contributed to the parallel 
peak model, is of significance in domain formation,49 but this is only one as- 
pect of the complex nature of soil structure. 

The writers would join Hough in emphasizing the mutual advantages of co- 
operation between soil scientists and engineers. In particular the first named 
writer would acknowledge the great benefit he has derived from association 
with soil scientists both in his own University and in the Waite Agricultural 
Research Institute, University of Adelaide, and the Division of Soils, CSIRO 
Australia. 

In his contribution Eremin infers that step-strain characteristics may be 
observed in a colloidal matrix. It may well bethat there are certain conditions 
in which a sensitive clay soii may progressively develop an altered structure 
analogous to that described in the paper, but the writers have no evidence of 
this. 

McRae has emphasized the need for further researchintothe question of the 
influence of soil structure on the general stress-strain properties of soils, 
particularly in relation to compacted materials. The writers heartily agree 
with this viewpoint. Despite significant advances in recent years inour under- 
standing of the internal fabric of soils, it appears that we are little advanced 
in the development of a satisfactory mathematical model for describing their 
general stress-strain properties. 

The complexity of this problem is obvious in that we are forced to deal with 
relatively large volume changes of the material under stress. In other engi- 
neering materials, for example, steel and concrete, this volume change is of a 
much lower order and is generally neglected. Hence, we cannot expect to ex- 
trapolate directly to soils experience with these other materials. Studies of 
soil structure in relation to stress-strain behaviour are therefore essential 
preliminaries to the development of a satisfactory general theory. 

The writers were most interested to read that their hypothesis has a great 
deal in common witha similar hypothesis developed at the Waterways Experi- 
ment Station, Vicksburg. McRae rightly infers that the step-strain behaviour 
was observed in samples compacted on the wet side of the optimum water con- 
tent. In the paper the writers express the view that because of the high mois- 
ture tensions associated with samples compacted on the dry side of the opti- 
mum water content, the clay matrix will develop stiff-brittle characteristics 
and, thus, under shear it is unlikely that the step-strain mechanism will de- 
velop. It should not be concluded, however, that step-strain behaviour is nec- 
essarily restricted to soils compacted on the wet side of the optimum water 
content. Similar behaviour may well occur in natural soils, as for example in 
the collapsible soils studied in South Africa.43,44 


43 “The Additional Settlement of Foundations Due to a Collapse by Structure of Sandy 
Subsoils on Wetting,” by J. E. Jennings and K. Knight, Proceedings, 4th Internatl. Conf. 
of Soil Mechanics and Foundation Engrg., 1957. 
44 “Studies of Shear Strength and Bearing Capacity of Some Partially Saturated 
Sands,” by A. A. B. Williams, Proceedings, 4th Internatl. Conf. of Soil Mechanics and 
Foundation Engrg., 1957. 
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FUNDAMENTAL ASPECTS OF THIXOTROPY IN SOILS# 


Closure by James K. Mitchell 


JAMES K. MITCHELL,® A. M. ASCE.—Newland and Allely raise several 
questions regarding the writers’ interpretation of thixotropic behavior and sug- 
gest alternate explanations for several observations. In the study of any phe- 
nomenon as complex as thixotropy, in which the investigator is forced to pos- 
tulate micro-scale mechanisms on the basis of macro-scale observations, dif- 
ferences of opinion are to be expected. The discussers’ interpretations form 
a useful addition to the paper because alternate points of view are helpful in 
gaining a true understanding of phenomena. Several of the questions raised by 
Newland and Allely are considered in the following paragraphs. 

Newland and Allely take exception to the writers’ suggestion that the re- 
sults of their experiments (24) that showed that a sensitivity was developed 
during consolidation of 2 remolded clay were possibly due to thixotropy. While 
they list six points in connection with their results, that, considered together, 
would indeed suggest that thixotropy alone could not account entirely for the 
behavior, there is no proof that thixotropy could not be responsible in part for 
the behavior. There is no reason to believe that thixotropic strengthening pro- 
cesses are not occurring simultaneously with the consolidation. Thi, tropy is 
usually defined as a phenomenon dependent solely on a lapse of time and gen- 
erally assumes conditions of constant composition and yolume. However, the 
mechanisms active in the system at rest may still be active during consolida- 
tion and shear. 

Newland and Allely’s second interpretation of the method used by Day (26) 
for the determination of pore water tension in clay pastes as a function of time 
after stirring isthe correct one. This is, in order to maintain equilibrium and 
prevent water from being drawn into the sample, it was necessary to increase 
the pore water tension with increasing time of rest. The results of similar 
measurements on clay pastes have been reported by J. H. Kolaian and P. F. 
Low (35) and these are in agreement with those of Day. Bishop, Alpan, Blight, 
and Donald (36) presents data which show the same type of behavior in com- 
pacted clays, and it is suggested that the decrease in pore pressure (increase 
in tension) with time after compaction is a possible cause of thixotropy. 

It has been suggested by Newland and Allely that this increase in tension 
may be caused by increased osmotic pressure differences and thus increased 
repulsive forces between particles. If this is the case, they believe that there 
is then no evidence for a net increase in the effective stress that could contri- 
bute to any gain in strength. If, as is usual, effective stress is regarded as the 
difference between total pressure and pore water pressure then there must be 


a June 1960, by James K. Mitchell (Proc. Paper 2522). 
5 Asst. Prof. of Civ. Engrg. and Asst. Research Engr., Univ. of California, Berkeley, 
Calif. 
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an increase in effective stress with time if the pore pressure is continually 
decreasing because total stress remains unchanged. This definition of effective 
stress does not include directly the effects of interparticle repulsive forces; 
however, it should be noted that, theoretically at least, an increase in the os- 
motic repulsive forces must be reflected by a decrease in the measured pore 
pressure (37). It should alsove noted that an effective stress increase at con- 
stant volume is only possible if there is an accompanying change in the clay 
and/or water structure such that the compressibility of the system is reduced 
sufficiently to prevent volume change. 

On the basis of recent analyses of effective and intergranular stresses, 
somewhat different conclusions are possible. Equations have been proposed 
(38, 39, 40) that include interparticle attractive and repulsive forces in the 
formulation of intergranular stress in saturated soils. These equations are of 
the general form 


in which 0; is the intergranular pressure; o represents total pressure; A is 
the interparticle attractive pressure; u denotes hydrostatic pore water pres- 
sure; and R is the interparticle repulsive pressure. 

It may be seen from Eq. 3 that if the repulsive forces tend to increase with 
time, equilibrium may be maintained by a corresponding decrease in the hy- 
drostatic pore water pressure without change in the intergranular pressure 0j;. 
Even if the intergranular stress remains constant, however, it does not nec- 
essarily meanthat the effective stress remains constant, if by effective stress 
is meant that stress that controls clay strength. In fine-grained soils where 
interparticle attractive and repulsive forces are active the terms intergranular 
stress and effective stress are not necessarily synonomous. It may be con- 
cluded, therefore, that if oj is truly the effective stress, if the decrease in pore 
water pressure with time in a thixotropic soil at rest is the reaction to an in- 
crease in repulsive forces, and if the attractive forces remain constant with 
time, then Newland and Allely’s contention that there is no increase in effective 
stress that could contribute to a gain in strength is correct. 

Newland and Allely present arguments to show that it is not likely that the 
pore water tension arises from a tendency of the pore phase to shrink. The 
writer is in agreement with these points. A point not considered in detail by 
either the discussers or the writer, however, is that the increase in water ten- 
sion may arise from a time-dependent change in the pore water structure it- 
self at constant density. As pointed out by Martin (41) present evidence con- 
cerning the nature of adsorbed water structure is conflicting, and much re- 
mains to be learned. 

The writer would be quick to agree with Newland and Allely that the con- 
clusion that thixotropic hardening results from a spontaneous change from a 
dispersed to a flocculated state is not based on direct observation. Even if 
methods forthe direct observation of the structure were available, it is doubt- 
ful that the change in structure would be readily detectable. Particle transla- 
tions and rotations of only a very small magnitude would be necessary to effect 
considerable changes in strength, because at the low water contents of com- 
pacted clays particles approach within a very few angstroms of each other. It 
can be stated, however, that the test results indicate that the clays behaved in 
all respects as if they underwent a change from a dispersed to a more floccu- 
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lent structure. The results of a recent study on the effect of aging on the per- 


* meability of a compacted thixotropic clay lend weight to the argument for an 


actual change from a dispersed to a more flocculent structure. It was found 
that the permeability of samples aged for 21 days prior to test was up to seven 
times greater than that for samples of the same density tested immediately 
after compaction. The very pronounced effect of structure on permeability has 
been discussed by Michaels and Lin (42) and it is shown that flocculent struc- 
tures may be many times more permeable than dispersed structures. Finally, 
in the writers view a change from a dispersed to a more flocculent structure 
presupposes an increase in number and intensity of interparticle bonds. 

Eremin notes that clarification of the conclusion -- “while thixotropic ef- 
fects may cause appreciable increases in strengthin terms of total stress, they 
may cause little or no effect in terms of effective stresses” -- is needed. Un- 
fortunately, insufficient data are available to permit drawing a more positive 
conclusion. The results in Figs. 14 and 15 show that the curves of effective 
principal stress ratio 01/03, versus axial strain are not appreciably different 
for the aged and remolded specimens, even though the curves for deviator 
stress versus axial strain show a significant difference. These results would 
tend to suggest that the thixotropic strength gain can be attributed to the lower 
pore water pressures developed in the aged samples. Such results are in a- 
greement with data presented by Seed and Chan (23) that show that while the 
strengths of flocculated samples are considerably greater than those of dis- 
persed samples if expressed in terms of total stresses, the strengths in terms 
of effective stresses are the same. In the writers’ opinion, however, confir- 
mation of this tentative conclusion with regard to thixotropic behavior should 
be based on the results of a determination of failure envelopes in terms of ef- 
fective stresses for both aged and remolded samples. Determination of the 
effective cohesion and the effective angle of friction according to the procedure 
described by Seed, Mitchell, and Chan (40) would also provide useful informa- 
tion. 
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EXPERIENCE WITH A PIER-SUPPORTED BUILDING OVER PERMAFROST 


Discussion by Edward F. Lobacz 


EDWARD F. LOBACZ.4—This paper is of particular interest to the writer 
because he had the opportunity to examine the pier-supported structure at 
Churchill, Manitoba and other structures in the townsite while at Churchill in 
July of 1955 on foundation investigations for the ICY Rocketry Program. The 
paper illustrates a satisfactory structure and foundation design that utilizes a 
gravel fill and an air space, the need for proper control of ground water and 
drainage, and good construction practices in permafrost areas. It is unfortu- 
nate that the taking of ground temperatures was discontinued so early after 
completion of the foundation, and that other interrelated factors suchas ground 
surface temperatures, degradation and/or aggradation of permafrost, and ver- 
tical movements were not programmed for correlation with structure behavior. 

As the authors indicate, very little attention was given to suitable foundation 
designs on permafrost until after World War II. Because of this, the Arctic 
Construction and Frost Effects Laboratory initiateda building foundation study 
in 1946 at the Arctic Construction Investigations Area in Fairbanks, Alaska, 
(known as the Fairbanks Permafrost Research Area prior to 1959). The study 
was undertaken to aid in abetter understanding of building performance and to 
assist inthe development of improved foundation designs and construction tech- 
niques for arctic and subarctic areas. 

This study was undertaken, in part, by the construction of eleven test build- 
ings. Various ground exposures and insulators were tried to determine the ef- 
fect that different types of foundation construction had on ground temperatures, 
especially in the permafrost, and on the vertical movements of the building 
foundations. 5 

The program of observations covered 5 yr on eight of the buildings and 10 
yr on the other three. Many of the foundation types investigated were found to 
be satisfactory for temporary buildings only. Twotypes were found acceptable 
for permanent-type buildings: one combining an air space and an underlying 
non-frost-susceptible fill and one with piles bonded in permafrost. The per- 
formance of a building with the former type foundation is given subsequently. 

The building is a 32 ft by 32 ft frame residence that was erected in 1947. 
Foundation details are shown in Fig. 10 together with the permafrost table lo- 
cation. As indicated, the permafrost table remained very stable at the northern 
end of the building, whereas a steady degradation took place at the southern end 
between 1948 and 1957. This lowering was attributed, in part, to a concentrated 


4 October 1960, by H. B. Dickens and D. M. Gray (Proc. Paper 2618). 

4 Coordinator and Head, Engrg. Proj. Sect., Arctic Constr. and Frost Effects Lab., 
U. S. Army Engr. Div., Waltham, Mass. 

5 “A Study of Building Foundations on Permafrost,” by M. S. Kersten and E. F. Lo- 
bacz. Presented at the October 15, 1958 Convention of the ASCE at New York City. 
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flow of ground water in a natural drainage ditch along the access road just 
south of the structure. Until the winter of 1950, the natural flow inthe drainage 
ditch was augmented by an estimated 200,000-gal per day discharge from an 
uncontrolled artesian well. Such thawing of the permafrost illustrates the need 
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FIG. 10.—FOUNDATION DETAILS AND LOCATION OF PERMAFRO®? TABLE 
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FIG. 11.—DISPLACEMENT OF THE CORNERS 


for locating buildings away fiom drainage ditches, small streams, or channels 
of ground water flow. The permafrost levels maintained under the north por- 
tion of the building are considered normal for this type building. 
Vertical-movement observations of the four corners of the building shown 
in Fig. 11, indicated that level variations followed a sine-wave pattern with an 
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amplitude of about 0.15 ft and a net average downward displacement of 0.2 ft. 
This building has shown only minor indications of structural strain, that is, 
some cracking in the plastered walls. 

It was concluded from the 10-yr program of observations that (1) an air 
space under a building can effectively dissipate heat and prevent serious per- 
mafrost degradation; and (2) that, if subgrade conditions are reasonably uni- 
form, a non-frost-susceptible fill can reduce the magnitude of seasonal heave 
and distribute its force sufficiently to prevent the vertical movements from 
damaging the structure. 

A 3-ft-high air space in conjunction with pipe piers was used with apparent 
success in the design and construction of a Communications Center and Quar- 
ters Building at Bethel, Alaska in 1957. This foundation design was, in prin- 
ciple, that used for the pier-supported building at Churchill. Based on personal 
observations, a brief resume of the foundation treatment follows. 

This wooden structure, 30 ft by 68 ft in plan, was designed and constructed 
under the supervision of the U.S. Army Engineer District in Alaska. The sub- 
grade materials at the site consisted of a 2-ft-thick silt mantle underlain with 
sandy silts and silty sands to an explored depth of 13 ft. The subgrade mate- 


_ Vials were first excavated to a depth of 11 ft and then a layer of classified fill 


was placed and compacted to a depth of 5 ft. Concrete footings were founded 
on this fill to support 6-in.-diameter, standard-pipe piers (Fig. 12). Finally, 
unclassified fill was placed and compacted around the foundation footings and 
pipe piers to a height of 6 ft. 

It is of interest to note that, because of thaw, ground water difficulties were 
encountered on the job site and were controlled by ditching until the classified 
fill was in place. 

To check the temperature of the fill materials and the material beneath the 
footings during construction of the structure, the resident engineer placed a 
thermocouple string that extended horizontally beneath 5 piers. 

In July, 1958 a close visual examination of the completed building by the 
writer did not indicate any structural distress (Fig. 13). Temperature ob- 
servations in September, 1958 indicated the temperature beneath the piers to 
be 32° F. Periodic observations are now being made (1961) to determine the 
success of the structure. 
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GENERALIZED SOLUTIONS FOR LATERALLY LOADED PILES* 


Discussion by Aleksandar B. Vesic 


ALEKSANDAR B. VESIC, 19 M. ASCE.—A generalized procedure for ap- 
proximate analysis of a laterally loaded pile together with additional analyti- 
cally compatible patterns of bending moments and deflections are presented in 
this interesting paper. The writer would like to point out to some limitations 
of the approach used and to present a few other comments of general interest. 

Approach to the Problem.—The problem of a laterally loaded pile is, toa 
certain extent, complicated by the fact that, sometimes, a part of the soil ad- 
jacent to the top of the pile may undergo considerably higher shearing strains 
than the rest of the mass. Still, setting aside the possibility of a basically dif- 
ferent treatment by a generalized theory of earth pressure, 20 the most plau- 
sible approach to the problem would be to consider the pile as an elastic bar 
surrounded by an elastic (or visco-elastic) continuum and loaded at its end. 

The solution, by the methods of Theory of Elasticity, of problems of this 
kind, even in a most elemental form, treating, for instance, the soil as a homo- 
geneous, isotropic, linearly deformable solid, presents serious difficulties, 21 
A considerable mathematical simplification of the problem is achieved by in- 
troducing the hypothesis that the magnitude and variation of the ratio of soil 
pressures and deflections along the pile are known. 

However, if such a simplification is used, as in the present paper, the de- 
termination of the ratio of pressures and deflections, called “soil modulus” or 
“modulus of soil reaction,” becomes a problem in itself, by no means less im- 
portant than the derivation of solutions based on the mentioned hypothesis. 

Unfortunately, there has been in the past, in the case of the related problem 
of beams on elastic foundation, a strong tendency among the mathematically 
minded investigators to concentrate their efforts only on the first part of the 
problem: derivation of solutions. As vividly described elsewhere, 2 the sec- 
ond part of the problem, determination of the modulus of soil reaction, has been 
so largely neglected that a widespread erroneous conception has been created 
among engineers that this modulus is a kind of soil constant. 

The general approach to the problem in this paper displays a similar tend- 
ency. Efforts appear tobe concentrated on computations of various patterns of 
influences along the pile, corresponding to different assumed variations of the 


& October 1960, by Hudson Matlock and Lymon C. Reese (Proc. Paper 2626). 

19 Assoc. Prof. of Civ. Engrg., Georgia Inst. of Tech., Atlanta, Ga. 

20 “The Single Pile Subject to Horizontal Force,” by P. W. Rowe, Géotechnique 6, 
1956, pp. 70-85. 

21 Thesis presented to the Technical Univ. of Belgrade in 1956 in partial fulfillment 
of the requirements for the degree of Doctor of Science. 

22 “Evaluation of Coefficients of Subgrade Reaction,” by Karl Terzaghi, Géotech- 
nique 5, 1955, p. 300. 
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soil reaction modulus with depth. There are no statements regarding the fac- 
tors that influence the magnitude and variation of that modulus, nor any exper- 
imentally observed data. 

This, in the writer’s opinion, greatly restricts the general value of the pa- 
per, although many of the comparisons presented may be helpful as aids to 
judgment. 

Modulus of Soil Reaction.—It is probably not necessary to re-emphasize 
here the multitude of physical and geometrical factors influencing the magni- 
tude and variation with depth of the soil reaction modulus. Most of these fac- 
tors have been noticed already by the earlier investigators inthis field. 22, 23,24 
Many of them have been mentioned or discussed in more recent papers, 3,5 
(The sometimes complex influences of adjacent piles and of the pile caps, as 
well as the axial stress conditions and method of placing of the piles are per- 
haps the least frequently considered. ) 

Nevertheless, it is necessary, thinking of a very simplified mechanistic pic- 
ture of trying to set up a mathematical model of the problem, to indicate the 
following parameters as being of utmost importance: 


1. Stress-strain characteristics of the soil; and 
2. Dimensions, shape and flexural rigidity of the pile. 


To illustrate this argument, it may be interesting to note that, in the ideal- 
ized case of a sufficiently long pile of diameter D and flexural rigidity Ep I in 
a homogeneous isotropic medium defined by a Young’s modulus E and a Pois- 
son’s ratio v, the following expression for the modulus of subgrade reaction K 
can be found:25 


(88) 


Eq. 88, originally derived and verified experimentally for beams on elastic 
subgrade,26 proved to be very useful in interpretation of pile test results in 
some cases in which the soil conditions were such that an assumption of con- 
stant E was justified. 

Some analogous relationships should be found betweenthe soil reaction mod- 
ulus and the soil and pile parameters in other cases. As pointed out in the pa- 
per, the case of a soil profile displaying a linear increase of the modulus of 
deformation with depth deserves, in this sense, a particular attention. 

The advantages of using, in interpretation of test results, as well as in de- 
sign problems, raticna! formulas for the quantity in question, instead of em- 
pirically observed values, cannot be overemphasized. 27 


23 “Lateral Pile Loading Tests,” by L. B. Feagin, Transactions ASCE, Vol. 102, 
1937, p. 236. 

24 “Symposium on Lateral Load Tests on Piles,” papers by G. A. McCammon and 
J. C. Ascherman, L. B. Feagin, L. T. Evans, G. P. Tschebotarioff, J. O’Halloran, A. A. 
Wagner and Sol M. Gleser, ASTM Spec. Technical Publication, No. 154, 1954. x 

25 “Bending of Beams Resting on an Isotropic Elastic Solid,” by Aleksandar B. Vesic, 
Proceedings ASCE, Vol. 87, No. EM 2, April, 1961. 

26 “Beam on Elastic Subgrade and the Winkler’s Hypothesis,” by Aleksandar B. 
Vesié, Proceedings, Fifth Internatl. Conf. of Soil Mechanics and Foundation Engrg., 
Vol. I, Paris, 1961. 

27 Discussion on “Mechanics of Fluid Turbulence,” by Th. v. Karman, Transactions 
ASCE, Vol. 102, 1937, pp. 522-523. 
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Dimensional Analysis and Conditions of Similarity.—An attempt was made 
in the paper to use the principles of dimensional analysis in order to establish 
non-dimensional relations for a laterally loaded piie. In this way Eqs. 9 to 13 
were found and served the immediate purpose. 

It should be noted, however, that due to the nature of the approach used, only 
a minor part of the problem of similarity of laterally !oaded piles was actually 
considered. No attention was paid to the problem of similarity of the soil re- 
action moduli (the magnitude of whichis certainly influenced by the size of the 
loaded piles), 28,22 

It is, thus, hard to imagine physical conditions under which the authors’ Eqs. 
9 to 13 can be satisfied. The solutions presented in the paper, are, in a sense, 
to be considered as only formally non-dimensional. Therefore, unless more 
complete analytical and experimental arguments are presented, the procedure 
proposed at the end of the paper should be applied only to analyses based on 
fullscale test results. 

However, when such results are availiable there are alternative, simpler 
procedures to obtain data necessary for design. Some ideas on this have been 
published elsewhere. 29 

Terminology and Notations.—Finally, a minor remark will be made con- 
cerning the term “soil modulus” and the notation Es. Although frequently used 
in the past, principally in this country, this term and notation do not conform 
with the recommended ASCE standards?9 and, in writer’s opinion, may lead to 
erroneous interpretations (Es isthe generally adopted notation for the modulus 
of deformation of the soil). The term “modulus of soil reaction” or “soil re- 
action modulus” and a notation like K would seem to be more appropriate. 


_28 Discussion on “Evaluation of Coefficients of Subgrade Reaction,” by P. W. Rowe, 
Geotechnique 6, 1956, p. 94. . 
29 “Contribution a l’étude des fondations sur pieux verticaux et inclinés,” by Alek- 


sandar B. Vesic, Annales des Travaux Publiques de Belgique, 1956, No. 6. 
30 “Glossary of Terms and Definitions in Soil Mechanics, ” Report of the ASCE Com- 


mittee, Proceedings, ASCE, Vol. 84, No. 4, October, 1958, pp. 1826-9 and 1826-23. 
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NEW METHOD OF CONSOLIDATION- COEFFICIENT EVALUATION@ 
Discussion by K. Y. Lo and M. Arnold 


K. Y. LO.8—An interesting and convenient method for the determination of 

the coefficient of consolidation has been given by the author. As the basis of 
this method rests on Terzaghi’s equation of consolidation, it is necessary to 
examine the assumptions made in the derivation of the equation sothat the im- 
plications involved in the method may be more fully explored. 
It is well known that both the permeability and compressibility of clays change 
with time during the process of consolidation. The nature of these variations 
has been studied by an approximate method.9 Thus, consolidation is essentially 
a nonlinear phenomenon. A. McNabb10 has derived the one-dimensional con- 
solidation equation in a very general form in which the only assumptions made 
are the validity of Darey’s Law, and the decrease in volume of the soil is equal 
to the quantity of water displaced at the boundary. To facilitate discussion, it 
is sufficient to write this nonlinear equation in the form 


8 a 


in which u is the excess pore pressure, z and t are the space and time varia- 
ble, and k, my and yw are the permeability, compressibility of soil, and density 
of water, respectively. Substituting cy = k/( Yw My), Eq. 18 may be written in 
the explicit form 


In Terzaghi’s equation, the right hand member of Eq. 19 is assumed tobe zero. 
In effect, the nonlinear Eq. 18 is linearized by making it piecewise linear so 
that both the permeability and compressibility are assumed constant during the 
consolidation process. 

Thus, it may be seen that the variation of the coefficient of consolidation 
during the process of consolidation may well be a physical phenomenon, and 
different values obtained at various times are not necessary due to errors in 


a February 1961, by R. F. Scott (Proc. Paper 2746). 

8 Civ. Engr., Norwegian Geotech. Inst., Blindern, Oslo, Norway. 

9 Discussion by K. Y. Lo of “Measurement of the coefficient of consolidation of la- 
custrine clay,” by Rowe, Géotechnique Vol. 10, No. 1, 1960, p. 36. 

10 “A mathematical treatment of one-dimensional soil consolidation,” by A. McNabb, 


Quarterly of Applied Mathematics, Vol. 17, No. 4, 1960, p. 337. 
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the initial dial reading as the author attributed. In fact, experience shows that 
the initial compression is only important for the first few increments, and cor- 
rections are unnecessary for the subsequent loadings if the tests are carefully 
performed and samples saturated. On the contrary, consolidation tests invar- 
iably show variation of coefficient of consolidation with effective stresses for 
successive increments. This fact implies that the coefficient of consolidation 
varies with time during the same increment asthe effective stresses increase. 

In comparing the different methods for the determination of coefficient of 
consolidation, it is interesting to note that Taylor’s method using the settle- 
ment versus square root of time plot has a sound theoretical basis. In the 
same paper McNabb has shown by dimensional analysis that even in the general 
case of nonlinearity, the settlement at the early stage proceeds in a linear re- 
lationship with t1/2 for boundary conditions as those in the oedometer test. 

It appears, however, that the author’s method is quite general and in case a 
tractable solution of Eq. 18 is obtained, the method may be adopted tothe more 
general case without difficulty. 


M. ARNOLD,!! M. ASCE.—As the author states, this method of computing 
the coefficient of consolidation would appear to have some value in providing 
spot-checks on either of the two methods commonly used, but the writer feels 
that it cannot supplant the older methods and that the author’s statement that 
continuous compression time readings are not required in somewhat mislead- 


ing. 

Using the author’s figures for the example on the one-dimensional consoli- 
dation test in which Taylor computed the coefficient of consolidation with the 
use of 10 dial readings up to a time of 20 min inhis square-root of time meth- 
od, the following observations may be made. 

As the author’s method is dependent on the value of the corrected zero dial 
reading, it is necessary to determine this withsome care and thus at least two 
evaluations are necessary, the minimum times for these being 0.25 min and 1 
min, and 1 min and 4 min. The corrected zero dial reading differs for each of 
these two combinations and it is, thus, unlikely that either of the two values, 
or their mean, would be completely correct; hence, it would be necessary to 
compute the coefficient of consolidation at least twice to obtain a reasonably 
correct value. Because the larger time has to be beyond the point of 50% con- 
solidation, the use of two further compression dial readings involves those at 
6 min and 9 min. 

Thus, the absolute minimum time required in this example is 9 min with a 
total of 6 dial readings (that at 2.4 min only would be unnecessary). 

The method is also sensitive to inaccuracies in dial readings as will be 
shown. If the compression dial reading for any particular time was to be in 
error by 5 units (anot improbable error where dial units represent 0.0001 in.), 
it is possible that the errors might be cumulative, see Table 3 for example. 

These values represent 74% and 123% of the author’s quoted value of 23.8 x 
10-4 sq cm per sec, and the discrepancies would also make it difficult to as- 
sess the probable inaccuracy in the corrected zero dial reading. It is felt, 
therefore, that more compression dial readings and further computations would 
be necessary to establish a reasonably accurate value for the coefficient and 
the method then loses the advantages that the author claims for it. 


11 Senior Lecturer in Civ. Engrg., Univ. of Adelaide, Adelaide, South Australia. 
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It should be noted alsothat the slope of the compression ratio curves at val- 
ues immediately beyond the point of 50% consolidation is very flat, particularly 
at low values of N, and it would appear that for any satisfactory degree of ac- 
curacy tobe obtained, it would be necessary to use values of N of not less than 
2.5. 

The direct use of numerical values of compression dial readings inthe com- 
putation of the coefficient of consolidation seems to be unjustified because of 


Compression 


TABLE 3,—CORRECTED ZERO DIAL READING = 2093 
Compression 
dial reading ratio 


Time factor 
(2) (3) (4) 


1820 0.685 0,186 
1695 


Cy, in sq cm 
per sec x 16-4 
(5) 


their divergence in practice from the theoretically exact values. In 1952 the 
writer developed a numerical method,12 but this was also highly sensitive to 
compression dial reading inaccuracies, and the only satisfactory way in which 
it could be used was to plot the dial readings and draw the best-fitting sraooth 
curve through the plotted points after which numerical values were taken from 
the curve. 


12 “The Calculation of the Coefficient of Consolidation of a Soil from Experimental 
Results,” by M. Arnold, Transactions, South African ICE. Vol. 2, No. 2, February, 
1952. 
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PROCEEDINGS PAPERS 


The technical papers published in the past year are identified by number below. Technical-division 
sponsorship is indicated by an abbreviation at the end of each Paper Number, the symbols referring to: Air 
Transport (AT), City Planning (CP), Construction (CO), Engineering Mechanics (EM), Highway (HW), Hy- 
draulics (HY), Irrigation and Drainage (IR), Pipeline (PL), Power (PO), Sanitary Engineering (SA), Soil 
Mechanics and Foundations (SM), Structural (ST), Surveying and Mapping (SU), and Waterways and Harbors 
(WW), divisions. Papers sponsored by the Department of Conditions of Practice are identified by the sym- 
bols (PP). For titles and order coupons, refer to the appropriate issue of “Civil Engineering.” Beginning 
with Volume 82 (January 1956) papers were published in Journals of the various Technical Divisions. To 
locate papers in the Journals, the symbols after the paper number are followed by a numeral designating 
the issue of a particular Journal in which the paper appeared. For example, Paper 2703 is identified us 
2703(ST1) which indicates that the paper is contained in the first issue of the Journal of the Structural Di- 
vision during 1961. 


VOLUME 86 (1960) 


JUNE: 2494(IR2), 2495(IR2), 2496(ST6), 2497(EM3), 2498(EM3), 2499(EM3), 2500(EM3), 2501(SM3), 2502 
(EM3), 2503(P93), 2504(WW2), 2505(EM3), 2506(HY6), 2507(WW2), 2508(PO3), 2509(ST6), 2510(EM3), 2511 
(EM3), 2512(ST6), 2513(HW2), 2514(HY6), 2515(PO3), 2516(EM3), 2517(WW2), 2518(WW2), 2519(EM3), 2520 
(PO3), 2521(HY6), 2522(SM3), 2523(ST6), 2524(HY6), 2525(HY6), 2526(HY6), 2527(IR2), 2528(ST6), 2529 
(HW2), 2530(IR2), 2531(HY6), 2532(EM3)°, 2533(HW2)°, 2534(WW2), 2535(HY6)¢, 2536(IR2)°, 2537(PO3)°, 
2538(SM3°°, 2539(STA)C, 2540(Ww2)°. 

JULY: 2541(ST7), 2542(ST7), 2543(SA4), 2544(ST7), 2545(ST7), 2546(HY7), 2547(ST7), 2548(SU2), 2549(SA4), 
2550(SU2), 2551(HY7), 2552(ST7), 2553(SU2), 2354(SA4), 2555(ST7), 2556(SA4), 2557(SA4), 2558(SA4), 
2559(ST7), 2560(SU2)°, 2561(SA4)°, 2562(HY7)°, 2563(ST7)C. 

AUGUST: 2564(SM4), 2565(EM4), 2566(ST8), 2567(M4), 2568(PO4), 2569(P04), 2570(HY8), 2571(EM4), 
2572(EM4), 2573(EM4), 2574(sM4), 2575(EM4), 2576(EM4), 2577(HY8), 2578(EM4), 2579(PO4), 2580 
(M4), 2581(ST8), 2582(ST8), 2583(EM4)¢, 2584(PO4)¢, 2585(SsTS)°, 2586(SM4)°, 2587(HY8)¢, 

SEPTEMBER: 2588(IR3), 2589(IR3), 2590(WW3), 2591(IR3), 2592(HW3), 2593(IR3), 2594(IR3), 2595(IR3), 2596 
(HW +), 2597(WW3), 2£98(IR3), 2599(WW3), 2600(WW3), 2601(WW3), 2602(WW3), 2603(WW3), 2604(HW3), 
2605(SA5), 2606(WW3), 2607(SA5), 2608(ST9), 2609(SA5)°, 2610(IR3), 2611(WW3)°, 2612(ST9)°, 2613(IR3)°, 
2614(HW3)°. 

OCTOBER: 2615(EM5), 2616(EM5), 2617(ST10), 2618(SM5), 2619(EM5), 2620(EM5), 2621(SsT10), 
2622(EM5), 2623(SM5), 2624(EM5), 2625(SM5), 2626(SM5), 2627(EMS), 2628(EM5), 2629(ST10), 
2630(ST10}, 2631( PO5)¢, 2632(%M5)°, %673(ST10,, 2634(ST10), 2635(ST10)°, 2636(SM5) ¢. 

NOVEMBER: 2637(ST11), 2638(ST11), 2639(C:03), 2640(ST11), 2641(SA6), 2642(WW4), 2643(ST11), 2644(HY9), 
2645(ST11), 2646(HY9), 2647(WW4), 2648(WW4), 2649(WW4), 2650(ST11), 2651(CO3), 2652(HY9), 2653(HY9), 
2654(ST11), 2655(HY9), 2656(HY9), 2657(SA6), 2658(WW4), 2659(WW4)°, 2660(SA6), 2661(CO3), 2662(CO3), 
2663(SA6), 2664(CO3)°, 2665(HY9)°, 2666(SA6)°, 2667(ST11)°. 

DECEMBER: 2668(ST12), 2669(IR4), 2670(SM6), 2671(iR4), 2672(IR4), 2673(IR4), 2674(ST12), 2675 
(EM6), 2676 (IR4), 2677(HW4), 2678(ST12), 2679(EM6), 2680(ST12), 2681(SM6), 2682(IR4), 2683(SM6). 
2684(SM6), 2685(IR4), 2686(EMG), 2( ./{EM6), 2688(EM6), 2689(EMS), 2690(EM6), 2691(EM6)°, 2692 
(ST12), 2693(ST12), 2694(HW4)°, 2695(IR4)°, 2696(sM6)°, 2697(ST12)°¢, 


VOLUME 87 (1961) 

JANUARY: 2698(PP1), 2699(PP1), 2700(HY1), 2701(SA1), 2702(SU1), 2703(ST1), 2704(ST1), 2705(SU1), 2706 
(HY1), 270%(HY1), 2708(EY1), 2709(PO1), 2710(HY1), 2711(HY1), 2712(ST1), 2713(HY1), 2714(PO1), 2715 
(ST1), 2716(HY1), 2717(SA1), 2718(SA1), 2719(SU1)°, 2720(SA1)°, 2721(ST1), 2722(PP1)°, 2733(PO1)°, 2724 
(HY1)°, 2725(ST1)°. 

FEBRUARY: 2726(WW1), 2727(EM1), 2728(EM1), 2729(WW1), 2730(WW1), 2731 (EM1), 2732(SM1), 2733 (WW1), 
2734(SM1), 2735(EM1), 2736(EM1), 2737(PL1), 2738(PL1), 2739(PL1), 2740(PL1), 2741 (EM1), 2742(8T2), 
2743(EM1), 2744(WW1), 2745(WW1), 2746(5M1), 2747(WW1), 2748(EM1), 2749(WW1), 2750(WW1)¢, 2751 
(EM1)°, 2752(6M1)°, 2753(PL1)°, 2754(8T2)°, 2755(PL1). 

MARCH: 2756(HY2), 2757(IR1), 2758(AT1), 2759(CO1), 2760(HY2), 2761(IR1), 2762(IR1), 2763(HY2), 2764 
(ST3), 2765(HY¥2), 2766(HW1), 2767(SA2), 2768(CO1), 2769(IR1), 2770(HY2), 2771(SA2), 2772(HY2), 2773 
(Coli), 2774(AT1), 2775(IR1), 2776(HY2), 2777(HY2), 2778(SA2), 2779(ST3), 2780(HY2), 2781(HY2)°, 2782 
(HW1)°, 2783(SA2)°, 2784(CO1), 2785(CO1)°, 2786(IR1)°, 2787(ST3)°, 2788(AT1)°, 2789(HW1). 

APRIL: 2790(EM2), 2791(SM2), 2792(SM2), 2793(SM2), 2794(SM2), 2795(SM2), 2796(SM2), 2797(SM2), 2798 
(EM2), 2799(EM2), 2800(EM2), 2801(EM2), 2802(ST4), 2803(EM2)°, 2804(SM2)°, 2805(ST4)°. 

MAY: 2806(SA3), 2807(WW2), 2808(HY3), 2809(WWw2), 2810(HY3), 2811(Ww2), 2812(HY3), 2813(ww2), 
2814(HY3), 2815(WW2), 2816(HY3), 2817(HY3), 2818(SA3), 2819(WW2), 2820(SA3), 2821(Ww2), 2822 
(ww2)°, 2823(HY3), 2824(SA3), 2825(HY3), 2826(SA3)°, 2627(HY3)°. 

JUNE: 2828(SM3), 2829(EM3), 2830(EM3), 2831(IR2), 2832(SM3), 2833(HW2), 2834(IR2), 2835(EM3), 2836 
(IR2), 2837(IR2), 2838(SM3), 2839(SM3)°, 2840(IR2)°, 2841(HW2)°, 2842(EM3)°, 2843(ST5), 2844(STS), 2845 
(STS), 2846(ST5)°. 


c. Discussion of severa! papers, grouped by divisions. 
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